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Multi-year sets of ground-based sun-photometer measurements conducted at 1 2 Arctic sites and 9 Antarctic sites 
were examined to determine daily mean values of aerosol optical thickness r(A) at visible and near-infrared 
wavelengths, from which best-fit values of Angstrom's exponent a were calculated. Analysing these data, the 
monthly mean values of t(0.50 pm) and a and the relative frequency histograms of the daily mean values of 
both parameters were determined for winter-spring and summer-autumn in the Arctic and for austral summer 
in Antarctica. The Arctic and Antarctic covariance plots of the seasonal median values of a versus t(0.50 pm) 
showed: (i) a considerable increase in t(0.50 pm) for the Arctic aerosol from summer to winter-spring, without 
marked changes in a\ and (ii) a marked increase in t(0.50 pm) passing from the Antarctic Plateau to coastal sites, 
whereas a decreased considerably due to the larger fraction of sea-salt aerosol. Good agreement was found when 
comparing ground-based sun-photometer measurements of r( A) and a. at Arctic and Antarctic coastal sites with 
Microtops measurements conducted during numerous AERONET/MAN cruises from 2006 to 2013 in three Arctic 
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Ocean sectors and in coastal and off-shore regions of the Southern Atlantic, Pacific, and Indian Oceans, and the 
Antarctic Peninsula. 

Lidar measurements were also examined to characterise vertical profiles of the aerosol backscattering coefficient 
measured throughout the year at Ny-Alesund. Satellite-based MODIS, MISR, and AATSR retrievals of r(X) over 
large parts of the oceanic polar regions during spring and summer were in close agreement with ship-borne 
and coastal ground-based sun-photometer measurements. An overview of the chemical composition of mode 
particles is also presented, based on in-situ measurements at Arctic and Antarctic sites. Fourteen log-normal 
aerosol number size-distributions were defined to represent the average features of nuclei, accumulation and 
coarse mode particles for Arctic haze, summer background aerosol, Asian dust and boreal forest fire smoke, 
and for various background austral summer aerosol types at coastal and high-altitude Antarctic sites. The main 
columnar aerosol optical characteristics were determined for all 14 particle modes, based on in-situ measure- 
ments of the scattering and absorption coefficients. Diurnally averaged direct aerosol-induced radiative forcing 
and efficiency were calculated for a set of multimodal aerosol extinction models, using various Bidirectional Re- 
flectance Distribution Function models over vegetation-covered, oceanic and snow-covered surfaces. These gave 
a reliable measure of the pronounced effects of aerosols on the radiation balance of the surface-atmosphere sys- 
tem over polar regions. 

© 2014 Elsevier B.V. All rights reserved. 
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1. Introduction 

Aerosols are one of the greatest sources of uncertainty in climate 
modelling, as their microphysical, chemical and optical characteristics 
as well as their concentration vaiy in time and in space, inducing signif- 
icant direct radiative forcing effects on the surface-atmosphere system. 
In addition, they can alter cloud optical features and indirectly impact 
climate. The aim of this paper is to present an overview of the optical 
characteristics of atmospheric aerosol observed in polar regions during 
the past two decades, including recent measurements conducted with 
ground-based and ship-borne sun-photometers, or retrieved from re- 
mote sensing data recorded with visible and infrared sensors mounted 
onboard various satellite platforms. Optical instruments (e.g., lidars, 
sun-photometers) measure the characteristics of the atmospheric light 
field (internal, reflected, or transmitted). Specific procedures therefore 
need to be applied to convert optical signals to aerosol characteristics, 
such as particle size and shape distributions, or chemical composition. 
Similar procedures are also needed to derive the vertical concentra- 
tion distribution from columnar measurements. They are based on the 


solution of the inverse problem of radiative transfer theory accounting 
for multiple light scattering, molecular and aerosol scattering and 
absorption, and surface reflectance effects. 

The presence of a visibility-reducing haze in the Arctic was already 
noted by early explorers in the 19th century (see Garrett and Verzella, 
2008, for a historical overview). The explorers also documented that 
haze particles were deposited on snow in remote parts of the Arctic 
(e.g., Nordenskidld, 1883) and haze layers were also observed later by 
pilots in the 1950s (Mitchell, 1957). The source of the haze was debated 
for almost a century but poorly understood until the 1970s when it was 
suggested that this “Arctic Haze” originated from emissions in northern 
mid-latitudes and was transported into the Arctic over thousands 
of kilometres (Rahn et al., 1977; Barrie et al., 1981). The seasonality of 
the haze, which peaks in winter and early spring, was explained by 
the fact that removal processes are inefficient in the Arctic during that 
time of the year (Shaw, 1995). 

Polar aerosols originate from both natural and anthropogenic sources 
(Shaw, 1988, 1995). In the Arctic regions, natural aerosols have been 
found to contain an oceanic sea-salt mass fraction that frequently 
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exceeds 50% on summer days, and a mass fraction of 30-35% due to 
mineral dust, with lower percentages of non-sea-salt (nss) sulphate, 
methane sulphonic acid (MSA), and biomass burning combustion prod- 
ucts. In contrast, anthropogenic particles have higher concentrations of 
sulphates, organic matter (OM) and black carbon (BC) with respect 
to natural aerosol (Quinn et al., 2002, 2007; Sharma et al., 2006). In 
fact, boreal forest fire (hereinafter referred to as BFF) smoke transported 
from North America and Siberia often contributes to enhance soot con- 
centration in summer (Damoah et al., 2004; Stohl et al., 2006). Rather 
high aerosol mass concentrations of anthropogenic origin are frequently 
transported from North America and especially Eurasia in the winter and 
spring months, leading to intense Arctic haze episodes (Shaw, 1995). For 
instance, Polissar et al. (2001) conducted studies on the BC source 
regions in Alaska from 1991 to 1999 finding that predominant contribu- 
tions have been given by large-scale mining and industrial activities in 
South and Eastern Siberia. In the North European sector of the Arctic, 
the dominant sources of sulphates and nitrates (and to a lesser extent 
of water-soluble OM and BC) are located in Europe and Siberia, due to 
both urban pollution and industrial activities (Hirdman et al., 2010). 
Episodes of Asian dust transport have also been observed over the 
past years in the North American sector of the Arctic, especially in spring 
(Stone et al., 2007), together with local transport of soil particles 
mobilised by strong winds, which provisionally enhance the mass 
concentrations of elemental components, such as Al, Si, Mg and Ca 
(Polissar et al., 1998). The Arctic atmosphere's stratification is highly sta- 
ble, with frequent and strong inversions near the surface, which limit 
turbulence and reduce the dry deposition of aerosols to the surface 
(Strunin et al., 1997). They also decouple the sea ice inversion layer 
from the Arctic free troposphere, leading to very different chemical and 
physical properties of aerosols in the sea ice inversion layer where 
aerosols are depleted, and higher up where a sulphate-rich background 
aerosol typically of anthropogenic origin is often found (Brock et al., 
2011). In addition, organic-rich biomass burning layers occur in the 
free troposphere but rarely reach the surface (Warneke et al., 2010; 
Brock et al., 2011). The low-altitude high-latitude atmosphere in the 
southern hemisphere is similarly stably stratified as in the Arctic's but 
is also influenced strongly by katabatic winds bringing down air from 
the high altitudes of interior Antarctica (Stohl and Sodemann, 2010). 

On larger scales, the stable Arctic stratification leads to the so-called 
polar dome, where isentropes form shells above the Arctic. As atmo- 
spheric transport tends to follow the isentropes, direct transport of air 
masses from mid-latitude pollution source regions into the Arctic 
lower troposphere is very inefficient. According to Stohl (2006), pollut- 
ed air masses from lower latitudes typically follow one of five major 
transport pathways (see Fig. 1): (1) lifting at the Arctic front, where 
wet scavenging is efficient; (2) lifting at lower latitudes (at the polar 
front or convection), where wet scavenging is even more efficient; 


(3) most importantly for Arctic surface aerosol concentrations, low- 
level transport over land in winter where strong radiative cooling 
allows air masses to enter the polar dome; (4) slow descent by radiative 
cooling of upper-tropospheric air masses into the polar dome; (5) slow 
mixing across the lateral boundaries of the dome. Forest or agricultural 
fires are important, as they produce strong aerosol plumes in the mid- 
to high-latitude free troposphere, which can subsequently enter the 
Arctic by one of the previously mentioned processes. 

In addition to long-range pollution transport, local emission sources 
can be important. For instance, emissions from cruise ships can lead 
to measurable enhancements of BC and other aerosols in the Svalbard 
archipelago (Eckhardt et al., 2013). Diesel generators can also locally 
pollute the environment (Hagler et al., 2008). Aircraft emissions north 
of the Arctic circle are primarily injected into the stratosphere where 
removal is inefficient and these emissions can slowly descend (Whitt 
et al., 2011). All these local sources can also enhance the Arctic aerosol 
background; however, quantification of their contribution relative to 
long-range transport of pollution from sources outside the Arctic 
remains uncertain. In addition, sulphate and volcanic ash from high- 
latitude eruptions can occasionally influence the Arctic troposphere 
(e.g., Hoffmann et al., 2010). The stratospheric background aerosol in 
the Arctic as elsewhere can be perturbed by explosive volcanic erup- 
tions (especially in the tropics) for several years. Some volcanic aerosol 
emission episodes have been observed by Bourassa et al. (2010) and 
O'Neill et al. (2012) over the last few years, involving the low strato- 
sphere over short periods of a few months. 

Transport processes in the high-latitude southern hemisphere are 
similar to those sketched for the Arctic in Fig. 1 . As in the northern hemi- 
sphere, polluted air masses from the lower-latitude continents are 
quasi-isentropically lifted to higher altitudes and, furthermore, there is 
no low-altitude transport pathway over land in winter (i.e., the ana- 
logue to transport pathway number 3 in Fig. 1 is missing in the southern 
hemisphere). Consequently, as in the Arctic, the lowermost troposphere 
in the Antarctic is very isolated and, thus, contains little anthropogenic 
pollution transported from lower-latitude continents (Stohl and 
Sodemann, 2010). A major difference to the Arctic, however, is the 
high topography of the Antarctic continent. This means that the most 
isolated air masses (as measured by the time since last exposure to pol- 
lution sources at lower-latitude continents) are not found close to the 
pole, as in the Arctic, but in the coastal areas surrounding Antarctica 
(Stohl and Sodemann, 2010). Descent over the Antarctic continent is 
stronger than over the Arctic and can also bring down air from the 
stratosphere, and air from the Antarctic interior is transported down 
to coastal areas by strong katabatic winds. 

In Antarctica, aerosols sampled at coastal sites originate almost 
totally from natural processes, with a prevailing oceanic sea-salt mass 
content of 55-60%, and lower percentages of nss sulphate (20-30%) 
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Fig. 1. Schematic of the main transport pathways of aerosols into the Arctic, as described by Law et al. (in press) [used with permission of the American Meteorological Society]. 
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and mineral dust (10-20%) (Tomasi et al., 2012). Only very low mass 
fractions of nitrates, water-soluble OM and BC have been monitored in 
Antarctica, mainly associated with transport from remote anthropogen- 
ic sources (Wolff and Cachier, 1998) or biomass burning (Fiebig et al., 
2009) in mid-latitude areas. More than 60-80% of particulate matter 
suspended over the Antarctic Plateau has been estimated to consist of 
nss sulphates formed from biogenic sulphur compounds and/or MSA, 
due to long-range transport in the free troposphere and subsequent 
subsidence processes. Therefore, aerosols sampled at these high- 
altitude sites contain only moderate mass fractions of nitrates, and 
minor or totally negligible mass percentages of mineral dust, water- 
soluble OM and BC (Tomasi et al., 2007, 2012). 

The paper is organised as follows. In the next section the ground- 
based remote sensing measurements of atmospheric aerosol are 
reviewed (sun-photometers, lidars). Section 3 gives a description of 
ship-borne aerosol remote sensing instruments. Section 4 discusses 
aerosol backscattering coefficient profiles from lidar measurements, 
while Section 5 is dedicated to airborne and satellite observations of 
polar aerosols. The last section presents the most important optical 
characteristics and size-distribution features of polar aerosols, which 
are appropriate for calculations of direct radiative forcing effects 
induced by aerosols on the climate system. 

2. Ground-based remote sensing measurements 

2.1. Spectral measurements of aerosol optical thickness 

Ground-based remote sensing of the optical characteristics of aero- 
sols in the atmospheric column is usually conducted with multi- 
wavelength sun-photometers. A sun-photometer is oriented towards 
the Sun to detect the solar radiation attenuated along the slant path 
from the top-of-atmosphere (TOA) to the ground. The atmospheric 
aerosol load leads to a decrease in the solar radiation transmitted 
through the atmosphere. This decrease depends on the aerosol optical 
thickness (hereinafter referred to as AOT and/or using symbol r(A)), 
which is given by the integral of the volume aerosol extinction coeffi- 
cient along the vertical path of the atmosphere. 

The networks, sites and sun-photometers whose data were em- 
ployed in this paper are defined and characterised in Tables 1 and 2 
for the Arctic and Antarctic regions, respectively. The largest networks 
of sun-photometers in the world are AERONET and SKYNET. Spectral 
measurements of r(A) are performed with AERONET sun-photometers 
(Holben et al., 1998) at 8 wavelengths ranging from 0.340 to 
1.600 pm, and with SKYNET instruments (Nakajima et al., 2007) at 
10 wavelengths from 0.315 to 2.200 pm. The Cimel CE-318 sun- 
photometers of the AERONET network are currently used at several 
Arctic sites: Barrow (since March 2002), Thule (since March 2007), 
Hornsund (since April 2005), Sodankyla (since February 2007), Tiksi 
(since June 2010), Resolute Bay (since July 2004), Eureka-OPAL (since 
April 2007) and Eureka-PEARL (since May 2007). The last three sites, 
located in the Nunavut region of Canada, are part of the AEROCAN/ 
AERONET network. In addition, an AERONET sun-photometer has 
been intermittently used since 2002 by the Atmospheric Optics Group 
(GOA) (University of Valladolid, Spain) at the Arctic Lidar Observatory 
for Middle Atmosphere Research (ALOMAR), located on the Andoya 
Rocket Range, near Andenes (Northern Norway) (Toledano et al., 
2012). AERONET sun-photometers have been used to obtain Level 2.0 
r(A) measurements at the South Pole (since November 2007), and 
occasionally at Dome Concordia (in January and December 2003) in 
Antarctica. AERONET Level 1.5 measurements of r(A) are available for 
McMurdo on the Ross Sea (from February to December 1997, in the aus- 
tral summer 2001/2002, and in January-February 2011), at Marambio 
(Antarctic Peninsula) since October 2007, at Vechernaya Hill (Thala 
Hills, Enderby Land) since December 2008, and at Utsteinen (Dronning 
Maud Land) since February 2009. However, because these data were 
not promoted to Level 2.0, they were not considered in the present 


study. A PREDE POM-OIL sun/sky radiometer was used in Antarctica 
during the 2001-2002 austral summer by Di Carmine et al. (2005) at 
the Mario Zucchelli station. PREDE instruments have been used since 
2001 in Antarctica at Syowa (East Ongul Island, Liitzow-Holm Bay) by 
the National Institute of Polar Research (NIPR, Tokyo, Japan) since 
2001, at Rothera by the British Antarctic Survey (BAS) since January 
2008, and at Halley by BAS since February 2009. In addition, Precision 
Filter Radiometer (PFR) sun-photometers (Wehrli, 2000; Nyeki et al., 
2012) from the Global Atmosphere Watch (GAW) PFR network are 
currently used in the Arctic at Summit by PMOD/WRC (Switzerland), 
at Ny-Alesund by NILU (Norway), and at Sodankyla by FMI (Finland), 
and in Antarctica at Marambio by FMI and at Troll by NILU. The instru- 
mental and geometrical characteristics of the PFR sun-photometer are 
described by Wehrli (2000). 

The monochromatic total optical thickness r T0T W of the atmo- 
sphere is commonly calculated in terms of the well known Lambert- 
Beer law for a certain sun-photometer output voltage J(A) taken within 
a spectral channel centred at wavelength A and for a certain apparent 
solar zenith angle 6 0 . The monochromatic value of r TO r(^) is given by 
(Shaw, 1976): 

r ror (A) = (l/m)ln[RJ 0 (A)/J(A)], (1) 

where:(i) m is the relative optical air mass calculated as a function of 6 0 
using a realistic model of the atmosphere, in which wet-air refraction 
and Earth/atmosphere curvature effects on the direct solar radiation 
passing through the atmosphere are properly taken into account 
(Thomason et al., 1983; Tomasi and Petkov, 2014); (ii)J(A) is the output 
signal (proportional to solar irradiance) measured by the ground-based 
solar pointing sun-photometer; (iii) JoW is the output signal that 
would be measured by the sun-photometer outside the terrestrial 
atmosphere, at the mean Earth-Sun distance; and (iv) R accounts for 
J 0 (A) variations as a function of the daily Earth-Sun distance (Iqbal, 
1983). 

The solar radiation reaching the surface for cloud-free sky conditions 
is attenuated not only by aerosol extinction but also by Rayleigh scatter- 
ing (Tomasi et al., 2010) as well as absorption by minor gases (mainly 
water vapour (H 2 0), ozone (0 3 ), nitrogen dioxide (N0 2 ) and its dimer 
(N 2 0 4 ), and oxygen dimer (0 4 )). The spectral values of r(\) within the 
main windows of the atmospheric transmission spectrum are according- 
ly calculated by subtracting the Rayleigh scattering and absorption optical 
thicknesses from r m AA). 

AOT is usually a smooth function of wavelength A (measured in pm), 
which can be approximated by the following simple formula: 

r(A) = r(A 0 )(A/A 0 )-“ (2) 

where a is the so-called Angstrom (1964) wavelength exponent, and A 0 
is usually assumed to be equal to 1 pm. In reality, the analytical form 
defined in Eq. (2) can be convex or concave depending on the relative 
contents of fine and coarse particles in the atmospheric column. 
O'Neill et al. (2001a) demonstrated that the variation in r( A) and 
its first and second spectral derivatives (named here a and a ' , re- 
spectively) can be realistically described in terms of the spectral in- 
teraction between the individual optical components of a bimodal 
size-distribution. O'Neill et al. (2001a) then showed that one can ex- 
ploit the spectral curvature information in the measured t(A) to permit 
a direct estimate of a fine mode Angstrom exponent [otf] as well as the 
optical fraction of fine mode particles. However, an analysis of a and a' 
determined in real cases and taking into account that both a(0.44- 
0.87 pm) and a' are closely related to the spectral features of r(A) 
showed that propagation of errors leads to an error Aa/a ~ 2 Ar(A)/ 
r( A) and an error Act'/ct' ~ 5Ar(A)/r(A), respectively (Gobbi et al., 
2007). These estimates yield values of Aa/a and Act'/ct' that are 
>~20% and >~50%, respectively, for r(A) < 0.10 and a typical sun- 
photometry error equal to -0.01. To avoid relative errors > -30% in 
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Table 1 

List of the Arctic stations, where regular ground-based sun-photometer measurements have been conducted over the past decades, using different instrument models equipped with a variable number of narrow-band interference filters to determine 
the spectral values of aerosol optical thickness r(A) and Angstrom wavelength exponent a in the visible and near-infrared wavelength range. 


Sun-photometer stations 

Managing institutions 

Geographical 

Overall number of 

Measurement period 

Sun-photometer 

Peak wavelengths (nm) 

Spectral 

References 



coordinates and 

measurement days 


model 

of the spectral channels 

interval 




altitude 





(nm) of a: 


Barrow, Alaska (USA) 

GMD/NOAA, Boulder, Colorado, USA 

71° 19' N, 156° 36' 

832 (cloud-screened 

March 2000- 

Carter Scott SP02 

412, 500, 675, 862 

412-862 

Stone (2002) 



W, 8 m a.m.s.l. 

by GMD/NOAA) 

September 2012 

Carter Scott SP01-A 

367,610, 778, 1050 

367-778 







Carter Scott SP022 

368,610, 778, 1050 

368-778 


Barrow, Alaska (USA) 

AERONET, NASA/GSFC, (Greenbelt, 

71° 19' N, 156° 40' 

579 (Level 2.0) 

March 2002- 

Cimel CE-318 

340, 380, 440, 500, 675, 870, 

440-870 

Holben et al. 


Maryland, USA); U.S. DoE Atmospheric 
Radiation Measurement Program, USA 

W, 0 m a.m.s.l. 


September 2013 


1020, 1640 


(1998) 

Resolute Bay, Nunavut 

AERONET/AEROCAN, Environment 

74° 44 ' n, 94° 54' W, 

361 (Level 2.0) 

July 2004- 

Cimel CE-318 

340, 380, 440, 500, 675, 870, 

440-870 


(Canada) 

Canada, Downsview, Ontario, Canada 

40 m a.m.s.l. 


October 2012 


1020, 1640 



Eureka-OPAL, Nunavut 

AERONET/AEROCAN, CARTEL, University 

79° 59' N, 85° 56' W, 

360 (Level 2.0) 

April 2007- 

Cimel CE-318 

340, 380, 440, 500, 675, 870, 

440-870 


(Canada) 

of Sherbrooke, Canada 

0 m a.m.s.l. 


September 201 1 


1020, 1640 



Alert, Ellesmere Island, 

GMD/NOAA, Boulder, Colorado, USA 

82° 28' N, 62° 30' W, 

810 (cloud-screened 

August 2004- 

Carter Scott SP02 

412, 500, 675, 862, or 368, 

412-862 or 

Stone (2002) 

Nunavut (Canada) 


210 m a.m.s.l. 

by GMD/NOAA) 

September 2012 


610, 778, 1050 

368-778 


Thule, NW Greenland 

AERONET, NASA/GSFC, Greenbelt, 

76° 31' N, 68° 46' W, 

605 (Level 2.0) 

March 2007- 

Cimel CE-318 

380, 440, 500, 675, 870, 1020 

440-870 

Holben et al. 


Maryland, USA 

225 m a.m.s.l. 


September 2012 




(1998) 

Summit, Central Greenland 

PMOD/WRC, Davos, Switzerland 

72° 35' N, 38° 28' W, 

391 (cloud-screened 

January 2001- 

PFR#34 

368,412, 500, 862 

412-862 

Wehrli (2000) 



3250 m a.m.s.l. 

by PMOD/WRC) 

October 201 1 





Ittoqqortoormiit, 

AERONET, NASA/GSFC, Greenbelt, 

70° 29' N, 21° 57' W, 

307 (Level 2.0) 

May 2010- 

Cimel CE-318 

340, 380, 440, 500, 675, 870, 1020 

440-870 

Holben et al. 

Eastern Greenland 

Maryland, USA 

68 m a.m.s.l. 


October 2013 




(1998) 

Ny-Alesund, Spitsbergen 

AWI, Bremerhaven, Germany 

78° 54' N, 11° 53' E, 

749 (cloud-screened 

April 2000- 

SP1A 

371, 380, 416, 443, 500, 532, 609, 

443-864 

Herber et al. 

(Svalbard, Norway) 


5 m a.m.s.l. 

by AWI) 

September 2013 

SP2H 

675, 778, 864, 1025, 1046, 1062 
367, 380, 413, 441, 501, 531, 605, 
673, 776, 862, 1023, 1045 

441-862 

(2002) 






STAR01 

390,441,501,531,605, 673, 
776, 862 

441-862 


Ny-Alesund, Spitsbergen 

N1LU, Kjeller, Norway 

78° 54' N, 11° 53' E, 

693 (cloud-screened 

March 2002- 

PFR#18 

367.7, 411.9, 500.6, 862.5 

411.9-862.5 

Wehrli (2000) 

(Svalbard, Norway) 


5 m a.m.s.l. 

by N1LU) 

September 

2004, and March 2006- 
September 2013 





Barentsburg, Spitsbergen 

IAO-SB-RAS, Tomsk, Russia 

78° 04' N, 14° 13' E, 

56 (cloud-screened 

April-August of 201 1 

New portable SPM 

339, 380, 442, 500, 547, 675, 

442-871 

Sakerin et al. 

(Svalbard, Norway) 


20 m a.m.s.l. 

by I AO) 

and 2012 


871, 1020, 1240, 1553,2134 


(2009, 2012, 
2014) 

Hornsund, Spitsbergen 

AERONET, NASA/GSFC, Greenbelt, 

77° 00' N, 15° 34' E, 

514 (Level 2.0) 

April 2005- 

Cimel CE-318 

380, 440, 500, 675, 870, 1020 

440-870 

Holben et al. 

(Svalbard, Norway) 

Maryland, USA; Institute of Geophysics, 
PAS, Warsaw, Poland 

10 m a.m.s.l. 


August 2013 




(1998) 

Sodankyla, 

FM1, Helsinki, Finland 

67° 22' N, 26° 38' E, 

312 (cloud-screened 

May 2004- 

PFR#32 

367.6, 411.4, 500.5, 861.6 

411.4-861.6 

Wehrli (2000) 

Northern Finland 


184 m a.m.s.l. 

by FM1) 

September 2013 





Sodankyla, 

AERONET, NASA/GSFC, Greenbelt, 

67° 22' N, 26° 38' E, 

119 (Level 2.0) 

February 2007- 

Cimel CE-318 

340, 380, 440, 500, 675, 870, 

440-870 

Holben et al. 

Northern Finland 

Maryland, USA; FM1, Helsinki, Finland 

184 m a.m.s.l. 


November 2013 


1020, 1640 


(1998) 

Tiksi, Northern-Central 

AERONET, NASA/GSFC, Greenbelt, 

71° 35' N, 128° 55' E, 

162 (Level 2.0) 

June-October of 2010, 

Cimel CE-318 

340, 380, 440, 500, 675, 870, 1020 

440-870 

Holben et al. 

Siberia (Russia) 

Maryland, USA 

0 m a.m.s.l. 


2011 and 2012 




(1998) 
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Table 2 

List of the coastal and high-altitude Antarctic stations, where regular ground-based sun-photometer measurements have been conducted over the two past decades, using different instrument models equipped with a variable number of narrow- 
band interference filters to determine the spectral values of background aerosol optical thickness r(A) and Angstrom wavelength exponent a in the visible and near-infrared wavelength range. 


Sun-photometer stations 

Managing institutions 

Geographical 

Overall number of 

Measurement 

Sun-photometer 

Peak wavelengths (nm) 

Spectral 

References 



coordinates and 

measurement days 

period 

model 

of the aerosol spectral 

interval 




altitude 




channels 

(nm) of a 


Marambio, Seymour- 

FMI, Heksinki, Finland 

64° 14' S, 56° 37' 

139 

August 2011- 

PFR #29 

367.6,411.4, 500.5, 861.6 

367.6-861.6 

Wehrli (2000) 

Marambio Island 


W, 205 m a.m.s.l. 

(cloud-screened by FMI) 

March 2013 





Neumayer, 

AW1, Bremerhaven, Germany 

70° 39' S, 8° 15' 

234 

September 2000- 

SP1A 

371, 380, 416, 443, 500, 532, 609, 

443-864 

Herber et al. (2002) 

Weddell Sea coast 


W, 40 m a.m.s.l. 

(cloud-screened by AWI) 

April 2007 

SP2H 

675, 778, 864, 1025, 1046, 1062 
367, 380, 413, 441, 501, 531, 605, 
673, 776, 862, 1023, 1045 

441-862 







STAR01 

390, 441, 501, 531, 605, 673, 776, 
862 

368.7,411.9, 500.6, 862.5 

441-862 


Troll, Queen Maud Land 

NILU, Kjeller, Norway 

72° 01' S, 2° 32' E, 

547 

January 2007- 

PFR#40 

411.9-862.5 

Wehrli (2000) 



1309 m a.m.s.l. 

(cloud-screened by NILU) 

April 2013 

PFR#42 

368.9,412.1,499.7, 862.2 

412.1-862.2 


Novolazarevskaya, 

AARI, San Petersburg, Russia; 

70° 46' S, 11° 50' 

83 (Level 1.5, cloud- 

December 2008- 

Hand-held 

440, 500, 675, 870 

440-870 

Smirnov et al. (2009) 

Queen Maud Land 

AERONET, NASA/GSFC, Greenbelt, 

E, 119 m a.m.s.l. 

screened by NASA/GSFC 

February 2009; 

Microtops, 





Maryland, USA. 



November 2009- 
February 2010 

calibrated at GSFC 




Mirny, Davis Sea coast 

AARI, St. Petersburg, Russia 

66° 33' S, 93° 01' 

725 (ABAS and SPM data 

March 2000- 

ABAS 

395, 408, 479, 581, 651, 789, 873, 

408-873 

Radionov et al. (2002), 



E, 40 m a.m.s.l. 

cloud-screened by AARI; 

October 2013 


1041 


Radionov (2005) 




Microtops data cloud- 


SPM 

340, 379, 443, 499, 548, 676, 871, 

443-871 

Sakerin et al. (2009, 




screened using the 



1019, 1244, 1555,2134 


2012), Tomasi et al. 




Smirnov et al. (2009) 





(2012) 




procedure) 


Microtops 

440, 500, 675, 870 

440-870 

Smirnov et al. (2009) 

Syowa, East Ongul Island, 

Japan Meteorological Agency 

69° 00' S, 39° 35' 

987 

January 2000- 

EKO MS-110 model 

368, 500, 675, 778, 862 

368-862 

Ohno (2005) 

Liitzow-Holm Bay 

(JMA), Tokyo, Japan 

E, 21 m a.m.s.l. 

(cloud-screened by JMA) 

December 201 1 

(with 2.5° field-of- 
view diameter) 




Mario Zucchelli, Terra 

ISAC-CNR, Bologna, Italy 

74° 42' S, 164° 07' 

87 (cloud-screened by 

November 2001- 

Prede POM-OIL 

400, 500, 675, 870, 1020 

400-870 

Di Carmine et al. (2005) 

Nova Bay, Ross Sea 


E, 15 m a.m.s.l. 

ISAC-CNR) 

February 2002 





coast, Victoria Land 




December 2005- 

ASP-1 5WL 

381, 412, 451, 500, 551, 610, 673, 

412-861 

Tomasi et al. (2007) 





February 2006 


775, 861, 1026 



Dome Concordia, East 

AERONET, NASA/GSFC, Greenbelt, 

75° 05' S, 123° 18' 

44 (Level 1.5 data, cloud- 

January and 

Cimel CE-318 

440, 675, 870, 1020 

440-870 

Holben et al. (1998), 

Antarctic Plateau 

Maryland, USA 

E, 3260 m a.m.s.l. 

screened by NASA/GSFC) 

December 2003; 
January 2004 




Six et al. (2005) 


GMD/NOAA, Boulder, Colorado, 

75° 06' S, 123° 21' 

65 (cloud-screened by 

January 2006- 

Carter Scott SP02 

412, 500, 675, 862 

412-862 

Stone (2002) 


USA 

E, 3233 m a.m.s.l. 

GMD/NOAA) 

November 2010 






OPAR Institute, (Univ. de la 

75° 05' E, 123° 18' 

39 (Level 1.5 data, cloud- 

January, 2010, 

Hand-held 

379,441,674, 868 

441-868 

Smirnov et al. (2011) 


Reunion - CNRS, Saint Denis de la 

E, 3260 m a.m.s.l. 

screened by NASA/GSFC) 

2011, and 2012 

Microtops calibrated 





Reunion, France), and NASA/GSFC, 




at the NASA/GSFC 





(Greenbelt, Maryland, USA) 




Facility 




South Pole, Antarctic 

GMD/NOAA, Boulder, Colorado, 

90° 00' S, 00° 00' 

1279 (cloud-screened by 

November 2001- 

Carter Scott SP02 

412, 500, 675, 862 

412-862 

Stone (2002) 

Plateau 

USA 

E, 2835 m a.m.s.l. 

GMD/NOAA) 

March 2012 






AERONET, NASA/GSFC, Greenbelt, 

89° 59' S, 70° 18' 

147 (Level 2.0 data, cloud- 

November 2007- 

Cimel CE-318 

340, 380, 440, 500, 675, 870, 1020 

440-870 

Holben et al. (1998) 


Maryland, USA; GMD/NOAA, 

E, 2850 m a.m.s.l. 

screened by NASA/GSFC) 

December 2012 






Boulder, Colorado, USA 
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Act' /a', Gobbi et al. (2007) suggested using only observations of t(A) > 
0.1 5. We applied the same criterion as a threshold for accepting outputs 
from the Spectral Deconvolution Algorithm (SDA) of O'Neill et al. 
(2003) (notably the fine mode Angstrom exponent, which offers an 
alternative refinement to the calculation of a): our logic being that a 
and a' are input parameters to SDA and we did not want to introduce 
unacceptable processing errors to the extraction of a spectral exponent 
indicator. These limitations to the use of spectral values of a and a' were 
also applied by Yoon et al. (2012), who only considered observations 
with r(0.44 pm) > 0.15 in order to avoid relative errors >30% in 
a'. Tomasi et al. (2007) showed that for the period 1977-2006, 
r(0.50 pm) did not exceed 0.15 for background summer aerosol condi- 
tions at Barrow, Alert, Summit, Ny-Alesund, Hornsund, Sodankyla and 
Andenes/ALOMAR in the Arctic, and was greater than 0.15 only during 
very strong episodes of Arctic haze in late winter and spring, and 
BFF smoke transport in summer. Similarly, Tomasi et al. (2012) found 
that (i) the measurements of r(0.50 pm) recorded from 2000 to 
2012 at Ny-Alesund were estimated to exceed 0.15 in summer (June 
to September) in only a few cases of strong transport of BFF smoke 
from lower latitudes; even in winter (December to March), they were 
higher than 0.15 only for 10% of the cases, typically associated with 
Arctic haze transport episodes; (ii) measurements of r(0.50 pm) > 0.15 
recorded at Barrow over the same period were observed only in a 
few percent of cases, as a result of Arctic haze; and (iii) daily mean 
background summer values of r(0.50 pm) measured at Tiksi in Siberia 
were always lower than 0.08 during the summer 2010. With regard to 
Antarctic aerosol, Tomasi et al. (2007) estimated that the daily mean 
values of t(0.50 pm) were lower than 0.10 during the austral summer 
months at Marambio, Neumayer, Aboa, Mirny, Molodezhnaya, Syowa, 
Mario Zucchelli, Kohnen, Dome Concordia and South Pole. In particular, 
examining the sun-photometer measurements carried out from 2005 
to 2010, Tomasi et al. (2012) reported that the austral summer values 
of r(0.50 pm) measured at Neumayer and Mirny were <0.10 during 
the whole season, while those measured at South Pole never exceeded 
0.06. 

Therefore, since the values of r(0.44 pm) determined from the sun- 
photometer measurements conducted in polar regions are mostly lower 
than 0.15 and thus below the value recommended by Gobbi et al. 
(2007), we have decided not to determine the exponent a f using the 
O'Neill et al. (2001b) algorithm, but to calculate the best-fit value of 
the Angstrom exponent a over the spectral range 0.40 < \ < 0.87 pm 
using Eq. (2). In real cases, the exponent a provides by itself a first 
rough estimate of the optical influence of the fine particle component 
on r(\), since it gradually decreases on average from cases where fine 
particle extinction predominates to cases where coarse particles are 
optically predominant. 

In all cases with relatively high values of r(0.50 pm), the AERONET 
and SKYNET sun/sky radiometers can also be used to provide regular 
measurements of sky-brightness along the solar almucantar (an azi- 
muthal circle around the local normal whose zenith angle equals the 
solar zenith angle 0 O ) and also in the principal plane containing the 
direction to the Sun and the local normal. The analysis of these measure- 
ments enables better constraints for solving the inverse problem (while 
accounting for multiple scattering and surface reflectance effects) and 
allowing the following parameters to be derived: (i) the aerosol single 
scattering albedo <»(A), (ii) the phase function P(0), (iii) the particle 
number density size-distribution N(r) = dN(r)/d(ln r )) given as a func- 
tion of particle radius r, and (iv) the complex refractive index n(\)-i 
k(\) (King and Dubovik, 2013). 

The data analysis performed in this paper was subject to certain data 
processing constraints across networks of instruments. In the first in- 
stance, all network protocols differ in many (typically) minor details 
such as the means of estimating molecular optical thicknesses and 
solar air masses, the nominal time interval between measurements 
and calibration protocols. In general, all data-sets were cloud-screened 
using temporal-based criteria that were developed and rigorously 


tested by each network group. Only Level 2.0 AERONET data were used: 
this corresponds to temporal-based cloud-screened data (Smirnov et al., 
2000) that has undergone a final quality assurance step. The GMD/ 
NOAA data acquired at Barrow and Alert were further cloud-screened 
using a spectral criterion wherein r( A.) spectra were eliminated for a < 
0.38. This added cloud-screening feature was found to be necessary in 
order to eliminate the influence of homogeneous, thin cirrus clouds that 
has escaped the temporal cloud-screening step. Finally we note that the 
specific spectral ranges of the a computations, while being nominally 
limited to 0.40 < A < 0.87 tun are given in Tables 1 and 2 for each type 
of instrument (the a regression was carried out for all wavelength 
channels between and including the wavelength extremes given in the 
“Spectral interval” column). 

Before presenting the evaluations of the seasonal variations in 
parameters r(0.50 pm) and a measured at the various Arctic and 
Antarctic sites, it seems useful for the reader to give a measure of the ex- 
perimental errors and variability features associated with the aerosol 
optical characteristics varying as a function of aerosol origin and their 
chemical composition evolutionary patterns: 

(i) the mean experimental error of AOT measured with the most 
sophisticated sun-photometers leads to an uncertainty of -0.01 
in the visible and near-IR, mainly due to calibration errors (Eck 
et al., 1999); 

(ii) the relative errors of exponent a determined in terms of 
Angstrom's (1964) formula is about twice the relative error of 
AOT (Gobbi et al., 2007), therefore leading to obtain at polar 
sites relative errors close to 30% and, hence, absolute errors of 
-0.50 (Mazzola et al., 2012); and 

(iii) the spread of a arising from the natural variability of Arctic and 
Antarctic aerosol types has been estimated by Stone (2002), 
Tomasi et al. (2007) and Treffeisen et al. (2007) to yield average 
uncertainties of ± 0.4 for Arctic summer background aerosols, 
±0.6 for Arctic aerosols including particular cases (like Asian 
dust and boreal smoke particles), and ±0.5 for austral summer 
background aerosol observed at Antarctic sites. 

2.1.1. Measurements in the Arctic 

In order to define the seasonal variations in polar aerosol optical 
characteristics, sun-photometer measurements of r(A) conducted 
at various visible and near-infrared wavelengths can be conveniently 
examined to evaluate the exponent a in terms of Eq. (2). Such measure- 
ments have been carried out at various Arctic and Antarctic sites during 
the past decades, providing useful information on polar aerosol optical 
and microphysical features. In fact, t(\) gives a measure of the overall 
aerosol extinction along the vertical atmospheric path, while a depends 
on the combination of the different extinction effects produced by nu- 
clei, accumulation and coarse particles. Values of a higher than 1.3 are 
usually observed in air masses where Aitken nuclei and very fine parti- 
cles (having radii r< -0.12 pm) optically predominate, while relatively 
low values of a < 1.0 are observed when accumulation (over the 
0.12 < r < 1.25 pm range) and coarse (over the r > 1.25 pm range) 
mode particles produce stronger extinction effects. The vertical profile 
of the aerosol volume backscattering coefficient can be determined by 
means of lidar measurements, allowing the classification of aerosol 
layers in the troposphere and lower stratosphere. 

Direct solar radiation measurements have been regularly conducted 
for cloud-free sky conditions at numerous polar sites over the past 
decades, using multi-spectral sun-photometers of different design. 
Tables 1 and 2 list the 12 Arctic sites and 9 Antarctic sites, whose data 
were employed in this study together with their geographical coordi- 
nates, measurement periods, the peak wavelengths of the spectral 
channels used to measure AOT and determine the best-fit value of a, 
and the main references where the technical characteristics of the 
instruments are detailed. The geographical locations of these sites are 
separately indicated in Fig. 2 for the Arctic and Antarctic regions. 
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The individual measurements of the spectral values of r(A) and of 
the exponent a obtained from the analysis of the field data recorded 
for cloud-free sky conditions were then averaged to yield daily means. 
Since the present analysis is devoted to tropospheric aerosols, the sun- 
photometer measurements conducted in the presence of stratospheric 
layers of volcanic particles were removed for the following sites and 
intervals at: (a) the Arctic sites, in May 2006 (Soufriere Hills eruption), 
October 2006 (Tavurvur eruption) (Stone et al., 2014), from mid-August 
to late September 2008 (Kasatochi eruption) (Hoffmann et al., 2010), 
from early July to early October 2009 (Saiychev eruption) (O'Neill 
et al., 2012), and in April 2010 (Eyjafjallajokull eruption), as well as 
the sun-photometer data collected at Barrow during the periods that 
followed both the Okmok eruption in July 2008, and the Mt. Redoubt 
eruption in March 2009 (Tomasi et al., 2012); and (b) the Antarctic 
sites, for all data affected by volcanic features comparable to those 
of Mt. Pinatubo observed from late spring 1992 to late autumn 1994 
(Stone et al., 1993; Stone, 2002). Actually, the Stratospheric Aerosol 
and Gas Experiment (SAGE II) observations made since 2000 over 
Antarctica did not provide evidence of appreciable extinction fea- 
tures produced by volcanic particle layers at stratospheric altitudes 
(Thomason and Peter, 2006; Thomason et al., 2008), as also con- 
firmed by the analysis of sun-photometer measurements conducted 
at Mirny, Neumayer, Mario Zucchelli and South Pole from 2000 to 
2010 (Tomasi et al., 2012). 

The remaining daily mean “tropospheric” data collected at each site 
were then subdivided into monthly sub-sets consisting of data mea- 
sured in different years, for which the multi-year monthly mean values 
of r(0.50 pm) and a were determined. Relative frequency histograms 
(hereinafter referred to as RFHs) were defined separately using the 
daily mean values of r(0.50 pm) and a collected at Arctic sites during 
the following seasons: (i) winter-spring, from December to May, 
when Arctic haze events were most frequent, and (ii) the summer- 
autumn, from June to October, to characterise background aerosols in 
summer. For the Antarctic sites, RFHs were determined for austral 


summer (from late November to February), to define the mean optical 
characteristics of background aerosols. 

2.1. 1.1. Measurements in Northern Alaska. Two multi-year sets of sun- 
photometer measurements from Barrow, located on the Arctic Ocean 
coast, were analysed in the present study: (a) the first series, acquired 
with the Carter Scott sun-photometer, was conducted from March 
2000 to September 2012 by GMD/NOAA (see Table 1) and consisted 
of spectral r(A) measurements, taken every minute on apparently 
cloud-free days and then cloud-screened by applying the GMD/NOAA 
selection procedure (for a(0.412 pm/0.675 pm) < 0.38); and (b) the sec- 
ond series was AERONET data collected from March 2002 to September 
2013 (see Table 1 ). The results are shown in Fig. 3 for the GMD/NOAA 
and AERONET measurements, where the data coverage is 94% and 78% 
of the overall 14-year period (March 2000 to September 2013), respec- 
tively. The GMD/NOAA monthly mean values of r(0.50 pm) increase 
from about 0.07 in February to more than 0.12 in April and May, and 
then gradually decrease to less than 0.1 0 in June and July, and to around 
0.05 in October, with (grey toned) standard deviations a T > 0.05 from 
April to July, and <~0.03 in the other months. Similar results were ob- 
tained for the AERONET measurements, which exhibited monthly 
mean values of t( 0.50 pm) that increased from about 0.12 in March to 
more than 0.16 in April, and then decreased to 0.10 in June and August 
and 0.05 in September and October, with o T > 0.05 from April to August, 
and <~0.03 for the other months. The monthly mean values of a deter- 
mined from the GMD/NOAA and AERONET measurements were rather 
stable from February to April, varying from 1 .1 0 to 1 .20, with a standard 
deviation o a = 0.3 on average, followed by a convex cap from April to 
August, with values close to 1 .50 from May to July, and increasing values 
from 1.30 to -1.50 in August-October. Fig. 3 also shows the RFHs of the 
daily mean values of r(0.50 pm) and a measured during the winter- 
spring and summer-autumn seasonal periods. The analytical curves 
drawn to represent the RFHs are normal curves and are normalised to 
yield unit (100%) integration over the measured sampling intervals of 




Fig. 2. Part (a): Map of the Arctic with the geographical positions of the ground-based sun-photometer stations (solid stars), labelled with the following numbered circles: (1) Barrow, 
(2) Resolute Bay, (3) Eureka-OPAL, (4) Alert, (5) Thule, (6) Summit, (7) lttoqqortoormiit, (8) Ny-Alesund, (9) Barentsburg, (10) Hornsund, (11) Sodankyla, (12) Tiksi, (13) Andenes/ 
ALOMAR, and (14) Kiruna. Grey symbols indicate the geographical positions of ship-borne sun-photometer measurements made on clear-sky days during the AERONET/MAN cruises 
in three different sectors: (i) Northern Greenland-Norwegian Sea (GNS), between 20° W and 30° E (squares), (ii) Barents Sea and West Siberian Sea (BWS), between 30° E and 130° E 
(diamonds), and (iii) Eastern Chuckci Sea, Beaufort Sea and Amundsen Gulf (NAA), between 170° W and 110° W (triangles). Part (b): as in the upper part, for the Antarctic ground- 
based sun-photometer and/or lidar stations (crosses), labelled with the following numbered circles: (1) Marambio, (2) Neumayer, (3) Troll, (4) Novolazarevskaya, (5) Mirny, 
(6) Syowa, (7) Mario Zucchelli, (8) Dome Concordia, (9) South Pole, (10) McMurdo, and (11) Dumont d'Urville. Grey symbols indicate the geographical positions of the ship-borne 
sun-photometer measurement days during the AERONET/MAN cruises, subdivided into the four following oceanic sectors: (i) Southern Indian Ocean (1ND), between 20° E and 150° E 
(squares), (ii) Southern Pacific Ocean (PAC), between 150° E and 75° W (upward triangles), (iii) Southern Atlantic Ocean (ATL), between 50° W and 20° E (diamonds), and (iv) Antarctic 
Peninsula (APE), between 75° and 50° W (downward triangles). 
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Fig. 3. Left-hand side: Time-patterns of the monthly mean values and standard deviations (defined by the grey shaded areas) of aerosol optical thickness r(0.50 pm) (open squares) 
and Angstrom wavelength exponent a (solid circles) obtained from both series of multi-year sun-photometer measurements conducted at Barrow (Alaska) by (a) GMD/NOAA (Boulder, 
Colorado, USA) over the period from March 2000 to September 2012, using the Carter Scott sun-photometers listed in Table 1, and (b) NASA/GSFC (USA) (in cooperation with the 
Brookhaven National Laboratory (Upton, NY, USA)) over the period from March 2002 to September 2013, using the AERONET Cimel CE-318 sun-photometer having the spectral charac- 
teristics given in Table 1. Right-hand side: Relative frequency histograms of r( 0.50 pm) and exponent a obtained separately for the winter-spring (Arctic haze) period from December to 
May (black colour) and the summer-autumn (background aerosol) period from June to October (grey colour). The seasonal mean values and 25th, 50th, and 75th percentiles of t( 0.50 pm) 
and a are reported in the boxes inserted into the graphs, as obtained by examining the daily mean values of both optical parameters measured at this coastal site of the Arctic Ocean in the 
winter-spring (black italics) and summer-autumn (grey italics) periods. 


r(0.50 (im) and a shown in Fig. 3. The RFHs for both instruments were 
very similar, although showing appreciable discrepancies between the 
means and percentiles, which are in general lower than the correspond- 
ing standard deviations. The seasonal mean values of t( 0.50 pm) were 
equal to 0.12 and 0.13 in winter-spring, for the GMD/NOAA and 
AERONET data, respectively, and close to 0.08 in summer-autumn for 
both data-sets. The RFHs also have long-tails towards high values for 
winter-spring data, and larger kurtosis in summer-autumn. The long- 
tail features could in part be ascribed to larger r(0.50 pm) values in 
April and May (ranging from 0.12 to 0.16) attributable to the frequent 
Arctic haze cases observed in spring but would also be, in part, due to 
the (asymmetrical) log-normal distribution that is arguably a better fit 
to the AOT RFHs in general (c.f. O'Neill et al„ 2000). 

The small discrepancies found between the time-patterns of the 
monthly mean values of r(0.50 pm) and a defined for both data-sets 
as well as those between the RFHs of both parameters might well be at- 
tributable to slight differences in the total observation periods of both 
sun-photometers and/or differences in GMD/NOAA and AERONET 
cloud screening. The different seasonal features of r(0.50 pm) and a 
shown in Fig. 3 arise mainly from the origins of the aerosol load, associ- 
ated with the transport of continental polluted air masses mainly from 
North America and Asia, in winter-spring (Hirdman et al„ 2010). 
It can also be seen in Fig. 3 that the left-hand wings of the RFHs for 
a contain some values <0.75: these are probably due to an optical 


predominance of coarse mode sea-salt aerosols and/or local blown 
dust. Similarly, a fraction of values with a < 1 .20 are presumably related 
to Asian dust transport episodes (Di Pierro et al., 2011) that are most 
frequently observed in March and April, and which are generally 
characterised by persistent extinction features typical of coarse mineral 
dust particles (Stone et al., 2007). 

2.3. 3.2. Measurements in Northern Canada (Nunavut). The results derived 
from the AERONET/AEROCAN measurements conducted at Resolute 
Bay and Eureka-OPAL, and those carried out by GMD/NOAA at Alert 
(Canada) (see Table 1 ) over the past decade are shown in Fig. 4, as ob- 
tained for air masses containing aerosols mainly transported from the 
North American and Arctic Ocean areas (Hirdman et al., 2010) over 
these three Canadian sites. The monthly mean values of t( 0.50 pm) 
exhibit: (i) rather high values related to Arctic haze in March-April, 
vaiying between 0.10 and 0.15, and (ii) relatively low values in the sub- 
sequent months, decreasing to around 0.05 in September at all the three 
sites. The month-to-month differences, vaiying on average from 0.05 at 
Resolute Bay and Alert in March-June to 0.01 at Alert in September, are 
similar to or smaller than the monthly mean values of <j r . The winter- 
spring RFHs exhibit higher mean values of r(0.50 pm), ranging from 
0.09 to 0.12, and broader right-hand wings compared to the summer- 
autumn period, where the RFHs are narrower, giving seasonal mean 
values varying from 0.06 to no more than 0.08. 
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Fig. 4. As in Fig. 3, for the multi-year sun-photometer measurements of aerosol optical thickness r(0.50 pm) and exponent a conducted at: (a) Resolute Bay by Environment Canada 
(Ontario, Canada) over the period from July 2004 to October 2012, using the AERONET/AEROCAN Cimel CE-318 sun-photometer having the spectral characteristics given in Table 1; 
(b) Eureka-OPAL by CARTEL (Sherbrooke University, Canada) from April 2007 to September 2011, using the AERONET/AEROCAN Cimel CE-318 sun-photometer having the spectral 
characteristics given in Table 1 ; and (c) Alert by GMD/NOAA (Boulder, Colorado, USA) from August 2004 to September 2012, using the Carter Scott SP02 sun-photometer, with the char- 
acteristics given in Table 1. 


Relatively small differences were found between the monthly mean 
values of a in winter-spring and summer-autumn, with: (i) stable 
values close to 1.50 in March-June at Resolute Bay, gradually decreasing 
to 1.0 in July-September; and (ii) values increasing from 1.0 to 1.5 in 
March-June at Eureka and Alert, and remaining close to 1.50 from July 
to September (with a a < -0.3). The monthly mean values of a deter- 
mined at Resolute Bay and Eureka-OPAL are appreciably higher than 
those measured at Barrow, presumably as a result of the weaker extinc- 
tion produced by coarse sea-salt particles and/or local dust, and weaker 
contributions of Asian dust in the spring months. Fig. 4 shows that the 
seasonal RFHs of daily mean a values do not exhibit symmetrical 


shapes: (i) the winter-spring RFHs are rather wide, showing long- 
tailed left-hand wings and mean values varying from 1.29 (at Alert) 
to 1.52 (at Resolute Bay), and (ii) the summer-autumn RFHs are 
characterised by mean values varying from 1.46 (Alert) to 1.64 
(Eureka). Only moderate, relative increases in the coarse particle con- 
tent occurred at Resolute Bay from winter-spring to summer, while 
greater winter-spring to summer-autumn variations were measured 
at Eureka and Alert presumably due to fine particle smoke transported 
from North American boreal forest fires (Stohl et al., 2006). Finally, 
it is worth noting that a few cases were found with a < 0.40 at Alert, 
despite the a < 0.38 cloud-screening rejection criterion. They were 
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presumably associated with prevailing extinction by coarse mode sea- 
salt aerosols transported from the Arctic Ocean, especially in the spring 
months. 

2AA.3. Measurements in Greenland. Fig. 5 shows the results derived from 
the multi-year sets of sun-photometer measurements conducted at: 

(i) the Thule AERONET site, in the north-western corner of Greenland; 

(ii) the Summit PMOD/WRC site, located in the middle of the Central 
Greenland ice sheet; and (iii) Ittoqqortoormiit AERONET site, on the 
eastern coast of Greenland (see Table 1 for details on these sites). 


The monthly mean r(0.50 pm) values at Thule decreased slowly 
from about 0.10 in April-May to around 0.05 in June-September, with 
o T ~ 0.03, while a was rather stable with monthly mean values ranging 
from 1.40 to 1.50 from March to September, with u a ~ 0.3. The winter- 
spring r(0.50 pm) RFH is characterised by a mean value of 0.093 and an 
asymmetric shape whose long-tailed, right-hand wing is influenced by 
the frequent occurrences of Arctic haze episodes. The summer-autumn 
RFH for r(0.50 pm) was more symmetric, with mean and median values 
equal to 0.058 and 0.049, respectively, and values of the 25th and 75th 
percentiles relatively close to the median value, as can be seen in 
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Fig. 5. As in Fig. 3, for the for the multi-year sun-photometer measurements of aerosol optical thickness t(0.50 pm) and exponent a conducted at: (a) Thule (North-western Greenland) by 
NASA/GSFC (USA) over the period from March 2007 to September 2012, using an AERONET Cimel CE-318 sun-photometer having the spectral characteristics reported in Table 1 ; 
(b) Summit (Central Greenland) by PMOD/WRC (Switzerland) from January 2001 to October 2011, using a PFR (No. N34) sun-photometer of the GAW-PFR network, having the spectral 
characteristics given in Table 1 ; and (c) Ittoqqortoormiit (Eastern Greenland) by NASA/GSFC (USA) over the period from May 2010 to October 2013, using an AERONET Cimel CE-318 sun- 
photometer having the spectral characteristics reported in Table 1. 
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Fig. 5. Very similar shapes of both seasonal RFHs for a were obtained, 
with mean values of 1.38 in winter-spring and 1.41 in summer- 
autumn, and similar values of the main percentiles from winter-spring 
to summer-autumn, indicating no relevant seasonal changes in the 
aerosol size-distribution. 

Rather stable monthly mean values of r(0.50 pm) were obtained at 
the high-altitude Summit station, equal to 0.05 ± 0.03 from March to 
August, and about 0.03 ± 0.01 in September and October. Both seasonal 
RFHs of the daily mean r(0.50 pm) values assumed very similar shapes, 
with mean values close to 0.05, and values of the main percentiles 
differing by no more than 0.01 from season to the other. Greater differ- 
ences were determined between the two seasonal RFHs of a, with mean 
values equal to 1.27 in winter-spring and 1.52 in summer-autumn, and 
the main percentiles differing by no more than 0.3. These aerosol optical 
characteristics indicate that Summit is representative of the Arctic free 
troposphere, influenced mainly by particulate transport from North 
America and Europe, and only weakly by Siberian aerosols (Hirdman 
et al., 2010). 

The monthly mean values of r(0.50 tun) at Ittoqqortoormiit showed 
similar seasonal variations to those at Thule, gradually decreasing from 
-0.08 in March to less than 0.04 in September-October, with a T = 0.05 
in spring, gradually decreasing in summer until reaching values of -0.01 
in autumn. The monthly mean a values increase from - 1 .00 in March to 
-1.50 in summer, and slowly decrease in September and October, vary- 
ing around 1 .40, with cr a = 0.2 in March and close to 0.3 in the other 
months. The RFHs for t(0.50 pm) did not vary largely from winter- 
spring to summer-autumn. The seasonal mean values were equal to 
0.068 in winter-spring and 0.052 in summer-autumn, which were 
not considerably different from the median values in both seasons. 
They give a measure of the appreciable decrease in t(0.50 pm) observed 
from winter-spring to summer-autumn. The seasonal RFHs for a 
are more similar to those obtained at Summit than those of Thule. 
In fact, the winter-spring mean value of a was equal to 1.28, and the 
summer-autumn value was equal to 1.45. These results suggest that 
the atmospheric content of fine mode particles increases considerably 
from winter to summer at Ittoqqortoormiit. Such variations are proba- 
bly associated with the marked extinction effects produced by maritime 
accumulation/coarse mode particles in winter-spring and the predom- 
inant aerosol extinction effects produced by background continental 
particles mainly transported in summer-autumn from Europe and 
North America and containing in general significant loads of both 
anthropogenic and BFF particles. 

The seasonal changes shown in Fig. 5 at the three Greenland sites 
can be mainly attributed to the variations in aerosol transport processes 
from anthropogenic/polluted regions or remote oceanic mid-latitude 
areas, and only rarely to Asian dust. Actually, the transport processes 
of anthropogenic soot aerosols are known to appreciably enhance 
r(\), yielding rather high values of a in general (Tomasi et al., 2007; 
Stone et al., 2008). This may also occur in free-tropospheric layers, as 
observed during airborne measurements conducted at mid-altitudes 
over the Arctic Ocean (Stone et al., 2010). 

2.1. 1.4. Measurements in Spitsbergen (Svalbard). The results obtained at 
Ny-Alesund, Barentsburg and Hornsund (Svalbard, Norway) from four 
different series of measurements are shown in Fig. 6, as conducted 
by AW1 (Bremerhaven, Germany) and NILU (Kjeller, Norway) at Ny- 
Alesund, IAO-SB-RAS (Tomsk, Russia) at Barentsburg, and the Institute 
of Geophysics (Warsaw University, PAS, Poland) in cooperation with 
NASA/GSFC (USA) at Hornsund (see Table 1). 

The AWI monthly mean values of r(0.50 pm) varied from 0.07 to 
0.09 in winter-spring, and considerably decreased in summer-autumn 
from about 0.05 in June to 0.02 in October, with oy = 0.03 on average. 
The NILU monthly mean values of r(0.50 pm) were equal to -0.10 
from March to May, and varied from 0.06 in June-July to less than 
0.05 in August-September, with oy equal to 0.04 in spring and 0.03 
in summer and autumn. The comparison between the AWI and NILU 


results shows good agreement, although the NILU values were occa- 
sionally higher than those for AWI by no more than 10% on average in 
spring, such discrepancies probably arising from the slightly dissimilar 
measurement periods of 14 and 8 years, respectively. The AWI monthly 
mean values of a increased from less than 1 .30 in March to 1 .50 in May, 
and slowly decreased in summer-autumn until reaching a value <1.20 
in September and becoming nearly equal to 1.50 in October, with 
<j a = 0.30 on average, while the NILU values increased slowly from 
1.20 in March to 1.50 in June, and slowly decreased to -1.30 in Septem- 
ber, with <j a = 0.20 in all months. These discrepancies of no more than 
15% are in general smaller than the monthly values of o a . Therefore, it is 
not surprising that the RFHs found for the daily mean r(0.50 pm) values 
derived from the AWI and NILU data-sets differ considerably from one 
season to another: (i) the AWI and NILU winter-spring mean values 
were equal to 0.082 and 0.089, respectively, with the main percentiles 
differing by no more than 0.007; and (ii) the AWI and NILU summer- 
autumn mean values of r(0.50 pm) were equal to 0.052 and 0.059, 
respectively, with a T = 0.04 on average, and having differences be- 
tween the main percentiles no greater than 0.02. The AWI seasonal 
RFHs of a yielded mean values of 1.32 in winter-spring and 1.28 in 
summer-autumn, while NILU RFHs gave mean values of 1.35 in 
winter-spring and 1.38 in summer-autumn. Comparable values of 
the three main percentiles of the AWI and NILU RFHs were also ob- 
tained, differing by less than 0.1 in winter-spring and less than 0.2 
in summer-autumn. Therefore, the analysis of the AWI and NILU 
RFHs of a showed: (i) more dispersed features of a in winter-spring, 
presumably due to the larger variability of fine and coarse particle 
concentrations during the frequent Arctic haze episodes, and (ii) 
values of a mainly varying from 1.00 to 1.60 in summer-autumn, 
due to the large variability of the fine particle mode atmospheric 
content (mainly related to BFFs smoke particle transport) with re- 
spect to that of accumulation/coarse mode particles (mainly of oce- 
anic origin). 

The data-sets of t( 0.50 pm) and a derived from the Barentsburg IAO 
measurements consisted of a number of daily measurements smaller 
than 10% of that given by the AWI and NILU Ny-Alesund measurements. 
The monthly mean t( 0.50 pm) values varied from 0.07 to 0.10, with 
a T = 0.02 on average, and those of a increased from -0.90 to 1.40 
during summer, with a a = 0.2. Therefore, these measurements differ 
only slightly from those measured at Ny-Alesund over longer multi- 
year periods. The RFHs of the daily mean values of r(0.50 pm) and a 
were prepared only for the summer-autumn period, and were found 
to have a seasonal mean value of r(0.50 pm) = 0.078, with the main 
percentiles differing by less than 0.02 from the mean. These values are 
appreciably higher than those determined at Ny-Alesund from the 
AWI and NILU data-sets. A summer-autumn mean value of a = 1.29 
was obtained, with only slightly differing values of the main percentiles, 
and a considerably narrower RFH curve than those from the AWI and 
NILU data-sets measured at Ny-Alesund. 

The Ny-Alesund results can also be compared in Fig. 6 with the 
AERONET measurements recorded at Hornsund. The Hornsund results 
are in close agreement with those of Ny-Alesund, since the Hornsund 
monthly mean values of r(0.50 pm) varied from 0.10 to no more than 
0.12 in winter-spring, with oy = 0.04 on average, and from 0.06 to 
0.08 in summer-autumn, with cr T < 0.03. The monthly mean values of 
a were rather stable at Hornsund, mainly ranging from 1.20 to 1.50, 
and showing small differences with respect to the AWI and NILU results. 
The Hornsund RFHs of r(0.50 pm) show that the winter-spring daily 
mean values were on average higher than those obtained in summer- 
autumn, by more than 0.03, with differences between the main seasonal 
percentiles no higher than 0.04. Small variations in a < 0.07 were ob- 
served from one season to the other between the seasonal mean values 
and the main percentiles. The long-tailed left-hand wings of the RFHs of 
a determined during both seasons suggest that important extinction ef- 
fects were presumably produced by the sea-salt accumulation/coarse 
mode particles during both seasons. 
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The seasonal variations in t( 0.50 |im) shown in Fig. 6 are mainly due 
to the different aerosol extinction features produced over the Svalbard 
region by Arctic haze, especially in spring, and by background aerosol 
in summer. They are in part associated with the significant seasonal 
mean decrease in the mean concentration of sulphate particles mea- 
sured within the atmospheric ground-layer. For instance, on the basis 
of long-range routine measurements of particulate chemical composi- 
tion conducted at the Zeppelin station (78° 58' N, 11° 53' E, 474 m 
a.m.s.l.), near Ny-Alesund (Svalbard), Strom et al. (2003) estimated 
that the mean mass concentration of nss sulphate ions decreases on 
average with season, changing from about 3 x 10 _1 pg m~ 3 in March- 
April (for frequent Arctic haze episodes) to around 5 x 10 -2 pg m“ 3 in 
late summer (for background aerosol conditions). These features arise 
from the fact that the frequent Arctic haze episodes observed in winter 
and spring over the Svalbard region are mainly due to aerosol transport 
from the Eurasian area, rather than from North America or East Asia 
(Hirdman et al., 2010). The region north of 70° N is isolated in summer 
from the mid-latitude aerosol sources, as demonstrated by Stohl et al. 
(2006), who analysed aerosol transport patterns into the Arctic. 
BFF smoke particles are episodically transported over the Svalbard region 
in summer, from the Siberian region and sometimes from North America 
(Tomasi et al., 2007; Stone et al., 2008). For instance, huge emissions 
from BFFs in North America reached Svalbard (Stohl et al., 2006) in July 
2004, while agricultural fires in Eastern Europe caused very strong pollu- 
tion levels in the Arctic during spring 2006 (Stohl et al., 2006; Lund 
Myhre et al„ 2007). 

2.1.1. 5. Measurements in Scandinavian and Siberian regions. Fig. 7 shows 
the time-patterns of the monthly mean values of r(0.50 pm) and a 
and the winter-spring and summer-autumn RFHs of both parameters, 
derived from the sets of FMI/PFR and AERONET sun-photometer mea- 
surements carried out at Sodankyla (Northern Finland), and the set of 
AERONET measurements conducted atTiksi in Northern-Central Siberia 
(Russia) (see Table 1). 

The FMI/PFR monthly mean values of r(0.50 pm) slowly increased 
from -0.05 in February to 0.08 in May, remained quite stable from 
June to August, then slowly decreased to 0.05 in September-October, 
with comparable values of a T ranging mainly from 0.04 to ± 0.06, with- 
out showing clear variations from winter-spring to summer-autumn. 
The monthly mean values of a increased from about 1.10 in February 
to over 1.50 in July, and then gradually decreased to 0.75 in November. 
Different time-patterns of the monthly mean values of r(0.50 pm) and 
a were obtained from the AERONET measurements conducted at 
Sodankyla over a shorter 7-year period, including only about a third 
of the daily PFR observations, and giving monthly mean values of 
r(0.50 pm) varying from 0.05 to 0.09 in winter-spring, and from 0.06 to 
0.11 in summer, which then decreased to less than 0.04 in September- 
October. The monthly mean values of a varied from 1 .20 to 1 .40 in 
winter-spring, increasing to more than 1.70 in July and decreasing to 
nearly 1.00 in October. 

To provide a more complete picture of the atmospheric turbidity 
features over Northern Scandinavia, the time-patterns of the PFR and 
AERONET monthly mean values of r(0.50 pm) and a obtained at 
Sodankyla are compared in Fig. 7 with those determined by Toledano 
et al. (2012) analysing the r(0.50 pm) and a data-sets collected at; 
(i) Kiruna (67° 51' N, 20° 13' E, 580 m a.m.s.l.) in Northern Sweden 
(270 km WNW from Sodankyla) using the GAW-PFR sun-photometer 
of the Swedish Meteorological and Hydrological Institute (SMHI) from 
2007 to 2010 by, and (ii) Andenes/ALOMAR (69° 18' N, 16° 01' E, 
380 m a.m.s.l.) in Northern Norway, using the AERONET/R1MA Cimel 


CE-318 sun-photometer from 2002 to 2007. The Kiruna PFR results 
are compared in Fig. 7 with those recorded at Sodankyla. The Kiruna 
monthly mean values of r(0.50 pm) were rather stable over the whole 
measurement period, with cr T varying from 0.01 to 0.04, while the 
monthly mean values of a increased from 1.10 in February to nearly 
1.60 in July, and then decreased gradually to -1.00 in August and 1.20 
in October, with values of u a varying from 0.1 0 to 0.25. Thus, the Kiruna 
monthly mean values of r(0.50 pm) closely agree with those measured 
at Sodankyla and only exhibit small differences between the August- 
October monthly mean values of a. 

A similar comparison is also made in Fig. 7 between the AERONET/ 
FM1 results obtained at Sodankyla and the ALOMAR results derived 
from the AERONET/RIMA Cimel CE-318 sun-photometer measurements 
at Andenes, from 2002 to 2007. The ALOMAR monthly mean values of 
r(0.50 pm) increased from about 0.04 in February to 0.13 in May, and 
then slowly decreased to around 0.11 in September-October, having 
values of a T mainly varying from 0.04 to 0.06. Therefore, the ALOMAR 
evaluations of r(0.50 pm) were in general considerably higher than 
those measured at Sodankyla, with differences comparable to the stan- 
dard deviations. The ALOMAR monthly mean values of a varied at 
Andenes from about 0.85 to 1.05 in February-May, increased in the fol- 
lowing months to -1.30, and subsequently decreased in late summer 
and autumn to reach a value close to 1 .00 in October, with o a varying 
mainly from 0.20 in winter-spring to 0.40 in summer-autumn. These 
findings indicate that the ALOMAR monthly mean values of a were con- 
siderably higher than the AERONET evaluations obtained at Sodankyla, 
by about 1 5% on average, presumably because of the more pronounced 
extinction effects by maritime particles. 

The time-patterns of the monthly mean values of r(0.50 pm) and 
a shown in Fig. 7, as obtained by us at Sodankyla, and at Kiruna 
and Andenes by Toledano et al. (2012) differ appreciably from those 
typically observed at the other Arctic sites located at higher latitudes 
and shown in Figs. 3-6. The reason for such differences is that the air 
masses reaching Northern Scandinavia during the year originate from 
the Eurasian continent and mid-latitude Atlantic Ocean in 56% of 
cases, and from the Arctic Basin and Northern Atlantic Ocean in the re- 
maining 44% (Aaltonen et al., 2006). Due to the alternation of polluted 
air masses from Eurasia with sea-salt particles from ocean areas, the 
monthly mean values of r(\) were rather stable over the entire year, 
while the monthly mean values of a were higher in early summer, 
when the Arctic Basin was the principal aerosol source. Because of the 
efficient transport processes taking place during the year from conti- 
nental polluted or ocean areas, the FMI/PFR and AERONET Sodankyla 
RFHs of r(0.50 pm) did not exhibit significant differences between the 
seasonal mean values and the main percentiles defined in winter- 
spring and summer-autumn. These had monthly mean values of 
around 0.08 and 0.07 in winter-spring, respectively, with the main per- 
centiles differing by no more than 0.01, and summer-autumn monthly 
mean values close to 0.07 on average, with differences of about 0.01 be- 
tween the main percentiles. The FMI/PFR and AERONET Sodankyla 
monthly mean values of a decreased by 0.10-0.20 on average, from 
winter-spring to summer-autumn. Clearer discrepancies over both sea- 
sonal periods were found, with FMI mean values of about 1.32 and 1.52 
in winter-spring and summer-autumn, respectively, and AERONET 
mean values equal to 1.23 and 1.44 in the same two seasons, providing 
similar values of the main seasonal percentiles. 

A comparison between the winter-spring and summer-autumn 
estimates of t( 0.50 pm) and a was not made at Tiksi, since AERONET 
sun-photometer measurements have been routinely conducted at this 
remote Siberian site only over the May-October period. The monthly 


Fig. 6. As in Fig. 3, for the sun-photometer measurements of aerosol optical thickness r(0.50 pm) and exponent a conducted on Spitsbergen Island (Svalbard, Norway) at: (a) Ny-Alesund 
by AWI (Bremerhaven, Germany) from April 2000 to September 2013, using the sun- and star-photometers listed in Table 1 ; (b) Ny-Alesund by NILU (Kjeller, Norway) from April 2002 to 
September 2004 and from March 2006 to September 2013, using a PFR (No. N18) sun-photometer of the GAW network; (c) Barentsburg by IAO-SB-RAS (Tomsk, Russia) in the April- 
August months of 2011 and 2012, using the SPM portable sun-photometer having the characteristics reported in Table 1, and (d) Hornsund by NASA/GSFC (USA) in cooperation with 
the Warsaw University (PAS, Poland from April 2005 to August 2013), using an AERONET Cimel CE-318 sun-photometer having the characteristics reported in Table 1. 
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Fig. 7. As in Fig. 3, for the sun-photometer measurements of aerosol optical thickness 7(0.50 pm) and exponent a conducted at: (a) Sodankyla (Northern Finland) by FMI (Helsinki, 
Finland) in the winter-spring and summer-autumn periods from late May 2004 to March 2013, using a PFR (No. N32) sun-photometer having the characteristics reported in Table 1, 
(b) Sodankyla by NASA/GSFC (USA) in cooperation with FMI (Helsinki, Finland) from February 2007 to November 2013, using an AERONET Cirnel CE-318 sun-photometer having the 
characteristics reported in Table 1 ; and (c) Tiksi (North-central Siberia, Russia) by NASA/GSFC (USA) in the summer-autumn months (June-October) of 201 0, 201 1 and 2012, using an 
AERONET Cimel CE-318 sun-photometer having the characteristics given in Table 1. In the first graphs, the time-patterns of the monthly mean values of 7(0.50 pm) and a (solid and 
open diamonds, respectively) estimated by Toledano et al. (2012) at Kiruna from PFR measurements conducted over the 2007-2010 period are shown with their standard deviations 
(vertical bars) for comparison with the FMI results found at Sodankyla. The time-patterns of the monthly mean values of 7(0.50 pm) and a (solid and open triangles, respectively) 
estimated by Toledano et al. (2012) at Kiruna from the Cimel CE-318 measurements conducted from 2002 to 2007 are shown in the seconda graph for comparison with those derived 
from the AERONET measurements carried out at Sodankyla. 


mean values of r(0.50 pm) exhibit a clear increase from -0.13 in May 
to 0.16 in July, with large average values of cr T = 0.14, followed by 
a marked decrease in the subsequent months to about 0.05 ± 0.03 in 
October. Such large variations in r(0.50 |im) were associated with a 
slow decrease in the monthly mean values of a from 1.9 in May to 
-1.5 in September, for which o a < 0.2. The rather high values of 
r(0.50 pm) and a determined in summer were probably due to the 


frequent BFF smoke transport episodes from the inner regions of Siberia. 
In fact, the summer-autumn RFH of t( 0.50 pm) exhibits a mean value 
close to 0.09 and a value of the 75th percentile equal to 0.12, which 
are both appreciably higher than those measured at Sodankyla in sum- 
mer. The RFH of a yields a mean value of 1.57, and shows a long-tailed 
left-hand wing with 25th percentile of 1.39, and a long-tailed right- 
hand wing with 75th percentile of 1.74. This rather high value is 


C. Tomasi et al. / Earth-Science Reviews 140 (2015) 108-157 


123 


probably associated with small fine particles generated by combustion 
processes, which dominate extinction effects. 

in summary, the Arctic results provide evidence of the seasonality 
of r(0.50 pm) and a. A scatter plot of the median values of a versus 
those of r(0.50 pm) is shown in Fig. 8a, separately for the winter- 
spring and summer-autumn seasonal periods. Fig. 8a shows that the 
median values of a vary from 1.10 to 1.70 over the whole year, with: 
(i) winter-spring median values of r(0.50 pm) ranging from 1.20 to 
1.50 over the 0.04-0.12 range, and (ii) summer-autumn median 
r(0.50 pm) values all smaller than 0.08 and mainly ranging from 1.30 
to 1.70. These features suggest that appreciable differences characterise 
aerosol extinction in: (a) winter-spring, when the median values of 
t( 0.50 pm) vary greatly from one Arctic site to another. This results 
from their dependency on the importance of particulate transport 
from the most densely populated mid-latitude regions towards the 
Arctic, which is particularly strong in late winter and early spring; and 
(b) summer, when the background aerosol composition varies from 



t( 0.50 |im) 



Fig. 8. Upper part (a) : Scatter plot of the seasonal median values of the Angstrom exponent 
a versus the corresponding seasonal median values of aerosol optical thickness r(0.50 pm) 
determined from the sun-photometer measurements listed in Table 1. Acronym key: 
Barrow (BRW-NOAA for the GMD/NOAA measurements, and BRW-AER for the NASA/ 
GSFC AERONET measurements), Resolute Bay (RES), Eureka-OPAL (E0P), Alert (ALE) 
Thule (THU), Summit (SUM), Ittoqqortoormiit (ITT), Ny-Alesund (NYA-AW1 and NYA- 
N1LU, for the AW1 and N1LU measurements, respectively), Barentsburg (BAR), Homsund 
(HOR), Sodankyla (SOD-PFR and SOD-AER, for the PFR and AERONET measurements, 
respectively) and Tiksi (T1K). The median values are represented using grey-solid and 
grey-and-white open symbols to represent the summer-autumn results (as shown in 
the legend) and solid or black-and-white open symbols to represent the winter-spring 
results, showing the background aerosol and Arctic haze optical characteristics, respec- 
tively. The 25th and 75th percentiles are used to define the limits of the vertical and hor- 
izontal dashed bars. Part (b), as in the upper part (a), for the winter-spring (grey symbols) 
and summer-autumn (open symbols) median values of a plotted versus the correspond- 
ing median values of 7(0.50 pm), as obtained from the seasonal sets of ship-borne sun- 
photometer measurements collected over the Arctic Ocean sectors: (i) GNS (Greenland 
Sea and Norwegian Sea) (diamonds), (ii) BWS (Barents Sea and West Siberian Sea) 
(circles), and (iii) NAA (North American Arctic ocean, including the East Chukchi Sea, 
the Bering Strait, the Beaufort Sea and the Amundsen Gulf) (squares). 


one site to another, as a result of different extinction characteristics of 
fine and coarse mode particles transported from remote regions. 

2.1.2. Measurements in Antarctica 

Ground-based sun-photometer measurements of aerosol optical 
parameters have been conducted in Antarctica during the short austral 
summer period. In the present study, nine multi-year sets of measure- 
ments made since 2000 have been analysed (see Table 2), collected at 
six coastal sites (Marambio, Neumayer, Novolazarevskaya, Mirny, 
Syowa, and Mario Zucchelli), a mid-altitude station (Troll) and two 
high-altitude sites (Dome Concordia and South Pole). 

2. 1.2.1. Measurements at coastal and mid-altitude sites. The results 
obtained from the sun-photometer measurements carried out at 
Marambio, Neumayer and Troll are presented in Fig. 9. The Marambio 
measurements were conducted from September to April, and provided 
monthly mean values of r(0.50 pm) varying from -0.03 in January to 
0.06 in March (with o T varying from 0.02 to 0.04) and a ranging from 
0.50 in March to 1.50 in November-January (with o a =0.50 on aver- 
age). The austral summer RFHs exhibit regular features with mean 
values of r(0.50 tun) = 0.038 and a = 1.20, probably due to sea-salt 
particles, which dominate extinction. The Neumayer measurements 
were conducted over the September-April period, showing rather sta- 
ble time-patterns of the monthly mean values of r(0.50 pm), ranging 
from 0.04 to 0.06 (with a T = 0.03 on average) and associated with 
very stable values of a varying from 0.50 to 1.00 (with u a = 0.3 on 
average), which indicate that aerosols are mostly of oceanic origin. 
The RFHs of both optical parameters are similar to those determined 
at Marambio, showing mean values of r(0.50 pm) = 0.041 and a = 
0.82, confirming that these stable extinction features are mainly pro- 
duced by sea-salt particles. The time-patterns of the monthly mean 
values of r(0.50 pm) and a measured at Troll, about 235 km from the 
Atlantic Ocean coast in the Queen Maud Land, were also quite stable 
from September to April, yielding values of r(0.50 pm) varying from 
-0.02 to 0.03 (with cr T < 0.005), and values of a slowly increasing 
from 1.25 in September to -1.50 in January, and then decreasing to 
1.00 in April. The RFH of r(0.50 pm) exhibits nearly symmetrical fea- 
tures with little dispersion, with a mean value of 0.023, and 25th and 
75th percentiles equal to 0.019 and 0.026, respectively, while the RFH 
of a was also quite symmetrical over the 0.60-2.10 range, with a 
mean value of 1.42 and 25th and 75th percentiles differing by less 
than 0.4 one from the other. These estimates of r(0.50 pm) and a differ 
appreciably from those obtained at Marambio and Neumayer, showing 
that the aerosol extinction features are only in part produced by mari- 
time aerosols at this mid-altitude site, and in part by fine mode particles, 
such as nss sulphate aerosols. 

Fig. 1 0 shows the results obtained from the measurements conduct- 
ed at Novolazarevskaya (Queen Maud Land) and Mirny (on the Davis 
Sea coast), using various sun-photometer models during different 
years, as reported in Table 2. The monthly mean values of r(0.50 tun) 
obtained at Novolazarevskaya were very close to 0.02 ± 0.01 over the 
whole period, while a varied from about 0.80 in November to 1.10 
in January, with o a equal to 0.10 in the first two months and 0.50 in 
January and February. The RFH of r(0.50 pm) exhibited a leptokurtic 
curve, with mean value close to 0.02, while a more dispersed distribu- 
tion curve was shown by the RFH of a, with the mean value close to 
unity, and 25th and 75th percentiles differing by less than 0.20 from 
it. The monthly mean values of r(0.50 pm) determined at Mirny varied 
from 0.02 to 0.03 (with a T = 0.01), and those of a from 1.50 to 2.00 in 
September-January, which then slowly decreased to -1.20 in April. 
The RFH of t( 0.50 pm) exhibited features which had a nearly symmetri- 
cal peak, with a mean value of 0.025, only slightly differing from that 
obtained at Novolazarevskaya, while the RFH of a was found to be 
dispersed and platykurtic, having a mean value of 1.60, and 25th and 
75th percentiles differing by more than 0.40 from the mean value. 
Therefore, it can be concluded that the aerosol extinction features 
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Fig. 9. As in Fig. 3, for sun-photometer measurements of aerosol optical thickness t(0.50 pm) and exponent a conducted during the September-April period at: (a) Marambio in the 
Seymour-Marambio Island (Antarctic Peninsula) by FMI (Helsinki, Finland), using a PFR (No. N29) sun-photometer from August 2011 to March 2013; (b) Neumayer on the Akta Bay 
(Weddel Sea coast) by AWI (Bremerhaven, Germany), using the sun-photometers SP1 A and SP2H, and the star-photometer STAR01 listed in Table 2 from September 2000 to April 
2007; and (c) Troll located at Jutulsessen (in the Queen Maud Land, 235 km from the coast, at 1309 m a.m.s.l.) by N1LU (Kjeller, Norway), alternately using two PFR (Nos. N40 and 
N42) sun-photometers from January 2007 to April 2013. 


shown in Fig. 10 are predominantly produced by sea-salt particles gen- 
erated by winds over the ocean, yielding values of a mainly ranging 
from 0.50 to 1.30. 

The results derived from the measurements conducted at Syowa 
and Mario Zucchelli are shown in Fig. 11, as obtained using various 
sun-photometer models over the different periods reported in Table 2. 
The monthly mean values of r(0.50 pm) varied from less than 0.02 
to -0.04 over the September-April period (with oy < 0.03), while 
the monthly mean values of a were also very stable, and described a 
large maximum of around 1.30 in December, with minima of -1.00 
in September and -0.90 in April. The RFH of r(0.50 pm) assumed a 


mesokurtic shape, skewed to the right, with the mean value close to 
0.02, while the RFH of a exhibited mesokurtic and symmetrical features, 
arising from sea-salt particles, which caused the predominant extinc- 
tion. The ISAC-CNR measurements conducted at Mario Zucchelli provid- 
ed the time-patterns of the monthly mean values of t( 0.50 pm) shown 
in Fig. 11, which were close to 0.02 in November and December (with 
oy = 0.01) and then increased to -0.04 in January and February (with 
cr T = 0.02), while the monthly mean values of a were stable and very 
close to 1.00 from November to January, and equal to -0.80 in February 
(with o a not exceeding 0.30). The RFHs of r(0.50 pm) and a were found 
to exhibit more dispersed features than those determined at Syowa and 
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Fig. 10. As in Fig. 3, for sun-photometer measurements of aerosol optical thickness 7(0.50 pm) and exponent a conducted at the Antarctic sites of: (a) Novolazarevskaya in the Schirmaker 
Oasis (Quenn Maud Land), 75 km from the coast by AARI (St. Petersburg, Russia), using a hand-held Microtops calibrated at GSFC (USA) over the periods from December 2008 to February 
2009 and from November 2009 to February 2010, and obtaining Level 1 .5 cloud-screened data; and (b) Mirny on the Davis Sea coast by AARI (St Petersburg, Russia), using the AARI, SPM 
and Microtops sun-photometers, having the spectral characteristics given in Table 2, over the period from March 2000 to October 2013. 


both Russian stations, giving mean values equal to 0.04 and 0.96, 
respectively. However, they clearly indicate that aerosol extinction is 
mainly due to sea-salt particles at this coastal site, yielding values of a 
ranging in general from 0.50 to 1.30. 

2. 3 .2.2. Measurements at the high-altitude sites on the Antarctic Plateau. 
The results obtained analysing the sun-photometer measurements con- 
ducted since 2000 at the Dome Concordia and South Pole high-altitude 
sites are shown in Fig. 12. The measurements were conducted by five 
groups using different instruments over distinct periods, as reported 
in Table 2. Due to the background transport of aerosols from very 
remote sources and the predominant role of subsidence processes 
on the aerosol load, the time-patterns of the monthly mean values 
of t( 0.50 pm) determined at Dome Concordia with different sun- 
photometers were found to be very stable, mainly ranging from 0.02 
to 0.04 (with a T evaluated to be <0.01 ) in September-April. The corre- 
sponding monthly mean values of a mainly varied from 1.00 to 2.00, 
with o a = 0.20 on average. The RFH of r(0.50 pm) exhibited a well- 
marked leptokurtic shape, with a mean value close to 0.02, while the 
RFH of a showed dispersed features over the 0.5-2.2 range, with a 
mean value close to 1.40, and 25th and 75th percentiles equal to 
about 1.0 and 1.80, respectively. 

The South Pole multi-year measurements conducted by GMD/NOAA 
at the Amundsen-Scott base were found to provide very stable time- 
patterns of the monthly mean values of r(0.50 pm), mainly ranging 
from September to March from -0.02 and 0.04 (with o T = 0.01 on 
average), and monthly mean values of a varying from 1.00 to 2.00, 
with the highest values in October and November (with o a = 0.10 on 


average). The AERONET monthly mean values of t( 0.50 pm) were 
found to increase from about 0.01 to 0.02 over the November-February 
period, with a T = 0.01 on average, the uncertainty of these measure- 
ments being primarily due to both calibration (estimated by Eck et al. 
(1999) to be of -0.01 in the visible), and forward scattered light enter- 
ing the instrument (Sinyuk et al., 2012). Very stable monthly mean 
values of a - 1 .00 were correspondingly found, with o a = 0.50 on aver- 
age. Both RFHs of t( 0.50 pm) derived from the GMD/NOAA and 
AERONET measurements assumed very narrow and “peaked" curves, 
with mean values of around 0.02, which appeared to be slightly skewed 
to the right. The corresponding RFHs of a presented dispersed features 
over the 0.50-2.00 range, with mean values of 1.54 and 1.06. Although 
such clearly dispersed results suggest the presence of an important 
fraction of large-size particles, it is important to take into account that 
a predominant particulate mass fraction of around 66% was estimated 
to consist of nss sulphates at South Pole, with lower concentrations of 
nitrates and sea-salt particles (Arimoto et al., 2004; Tomasi et al. 
2012), mainly associated with the background transport of aerosols 
from very remote sources and the strong effects exerted by subsidence 
processes. A few cases showing values of a < 0.50 were probably caused 
by consistent loads of diamond dust ice-crystals within the lower plan- 
etary boundary layer (hereinafter referred to as PBL), due to wind 
mobilisation. 

The covariance of the median values of t( 0.50 pm) and a obtained at 
the nine Antarctic sun-photometer stations is shown in Fig. 13a. There 
are two clusters: (i) the first at coastal sites, where the median values 
of r(0.50 pm) were in general considerably higher (at least by a factor 
two) than those found at the high-altitude Antarctic sites, and the 
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Fig. 11. As in Fig. 3, for sun-photometer measurements of aerosol optical thickness r(0.50 pm) and exponent a conducted at the Antarctic sites of: (a) Syowa (East Ongul Island, Lutzow-Holm 
Bay) over the period from January 2000 to December 201 1 by the Office of Antarctic Observation (Japan Meteorological Agency, Tokyo, Japan), using the EKO MS-1 1 0 sun-photometer, having 
the spectral characteristics given in Table 2; and (b) Mario Zucchelli on theTerra Nova Bay (Ross Sea, Victoria Land) during the austral summer periods of 2001/2002 and 2005/2006 by ISAC- 
CNR (Bologna, Italy), using the PREDE POM-01 L and the ASP-15WL sun-photometers, having the spectral characteristics given in Table 2. 


median values of a are in general smaller than 1 .20, and (ii) the second 
cluster, with lower median values of r(0.50 pm) and median values of 
a> 1.40, for which aerosol extinction is predominantly due to sea-salt 
coarse particles. It is interesting to note that lower median values of 
r(0.50 pm) associated with higher values of a were observed at the 
coastal sites of Mirny and Marambio, probably because of aeolian fine 
particles, which dominated extinction effects when compared to sea- 
salt particles. The median values of r(0.50 pm) and a determined at 
the Troll mid-altitude station indicate that the extinction effects exhibit 
intermediate features at this non-coastal site. This can probably be 
attributed in part to sea-salt aerosol and in part to background nss 
sulphate aerosol transported over large distances, which themselves 
contribute towards an appreciable enhancement in a. As can be seen 
in Fig. 13a, the weak median values of r(0.50 pm) and relatively strong 
median values of a at Dome Concordia and South Pole (the latter being 
partially obscured by the DMC-OPAR symbol) show that atmospheric 
aerosol extinction is rather weak above the Antarctic Plateau and mainly 
produced by fine particles consisting of nss sulphates. The lower ex- 
tremes of the a error bars and the low a value of the AERONET South 
Pole point are likely due to near-surface diamond dust layers that are 
observed on windy, cloud-free days at these sites. 

3. Ship-borne measurements 

A large number of ship-borne measurements of the main columnar 
aerosol radiative parameters r(A) and a were made over the Arctic 
and Antarctic oceanic regions during the past 10 years. The cruises 
were conducted by research groups and institutions from different 


countries, using hand-held Microtops II sun-photometers calibrated 
at the NASA/GSFC calibration facility. As can be seen in Fig. 2, large 
areas of the polar oceans were studied during these cruises, yielding 
an exhaustive picture of aerosol optical characteristics in these high- 
latitude remote regions. These ship-borne sun-photometer measure- 
ments were conducted as a part of the activities promoted and devel- 
oped in the framework of the Maritime Aerosol Network (MAN), 
which is a component of the AERONET network (see http://aeronet. 
gsfc.nasa.gov/new_web/maritime_aerosol_network.html). The mea- 
surements of r{\) were carried out at visible and near-infrared wave- 
lengths. Arctic Ocean cruises are listed in Table 3, while those taken in 
the oceanic regions around Antarctica are listed in Table 4. The hand- 
held Microtops II sun-photometers are equipped with five narrow- 
band interference filters with peak wavelengths of 0.340 (or 0.380), 
0.440, 0.675, 0.870 and 0.936 pm for most cruises, or 0.440, 0.500, 
0.675, 0.870 and 0.936 pm in other cruises, as reported in Tables 3 
and 4. The first four filters of each spectral set were used to measure 
the spectral values of r(\) and determine the best-fit value of exponent 
a over the wavelength range from 0.440 to 0.870 pm, while the mea- 
surements taken at 0.936 pm were combined with those made within 
the nearby channel centred at 0.870 pm to evaluate precipitable water 
(Smirnov et al., 2009). Among the data-sets available at the MAN 
website, we only selected measurements conducted since 2006 at lati- 
tudes higher than 67° N in the Arctic and within ocean areas far by no 
more than 1 000 km from the Antarctic coasts. In addition, the measure- 
ment sets obtained by Tomasi et al. (2007) during the RV Oceania 
cruises made in the summer months of 2003 and 2006 were included 
in the present analysis. 
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Fig. 12. Upper part: as in Fig. 3, for sun-photometer measurements of aerosol optical thickness r(0.50 pm) and exponent a conducted at: (a) the high-altitude site of Dome Concordia 
(DomeC), on the Eastern Antarctic Plateau, over the period from September to April, by (i) GMD/NOAA (Boulder, Colorado, USA) from January to November 2010, using a Carter Scott 
SP02 sun-photometer (squares); (ii) NASA/GSFC (USA) in cooperation with LGGE/CNRS (Grenoble, France) in January and December of 2003, and in January 2004, using an AERONET 
Cimel CE-318 sun-photometer having the spectral characteristics reported in Table 2 (circles); and (iii) OPAR Institute (University of Reunion, St. Denis, France) in January of 2010, 
201 1 and 2012, using an hand-held Microtops II sun-photometer calibrated at GSFC (USA), obtaining Level 1.5 cloud-screened data (triangles). Lower part: as in Fig. 3, for sun-photometer 
measurements of aerosol optical thickness r(0.50 pm) and exponent a conducted at the high-altitude site of South Pole (SPO) by (i) GMD/NOAA (Boulder, Colorado, USA), using a Carter 
Scott SP02 sun-photometer in the September-March period from November 2001 to March 2012; and (ii) NASA/GSFC (USA) in cooperation with GMD/NOAA (Boulder, Colorado, USA) 
using an AERONET Cimel CE-318 sun-photometer having the spectral characteristics reported in Table 2 in the November- February period from November 2007 to December 2012. 


3.1. Aerosol measurements in the Arctic Ocean 

Fourteen Arctic Ocean cruises were selected for this study, as 
summarised in Table 3. The data-set collected during the cruise of RV 
Polarstern (Alfred Wegener Institute (AWI), Bremerhaven, Germany) 
in August-September 2009 was not examined to evaluate the back- 
ground tropospheric aerosol extinction features, because these mea- 
surements were found to be strongly affected by the extinction of 


stratospheric volcanic aerosols generated by the Sarychev eruption in 
July-October 2009. No evidence of Sarychev volcanic particle extinction 
at stratospheric levels was found in CALIPSO (Cloud-Aerosol Lidar and 
Infrared Pathfinder Satellite Observations) data recorded over the 
Beaufort Sea (O'Neill et al., 2012) during the CCGS Amundsen cruise of 
August 2009. Therefore, this set of measurements was used to calculate 
the extinction effects associated in August 2009 with the BFF smoke par- 
ticles transported from the North American forests over the Arctic 
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Fig. 13. Part (a): as in Fig. 8, but for the sun-photometer measurements listed in Table 2. 
Acronym key: Marambio (MAR), Neumayer (NEU), Troll (TRO), Novolazarevskaya 
(NOV), Mirny (MIR), Syowa (SYO), Mario Zucchelli (MZS), Dome Concordia (DMC- 
NOAA, DMC-AER, and DMC-OPAR, for the measurement sets collected by GMD/NOAA, 
AERONET and OPAR groups, respectively (see Table 4)) and South Pole (SPO-NOAA and 
SPO-AER, for the measurement sets collected by GMD/NOAA and AERONET groups, 
respectively). Part (b): as in part (a), for the Microtops sun-photometer measurements 
performed in Antarctic coastal (open symbols) and off-shore areas (solid symbols) during 
the cruises conducted in the four following oceanic sectors: (i) IND (Southern Indian 
Ocean, squares), (ii) PAC (Southern Pacific Ocean, upward triangles), (iii) ATL (Southern 
Atlantic Ocean, diamonds), and (iv) APE (Antarctic Peninsula, downward triangles). 


Ocean. In order to study the aerosol extinction features of the Arctic 
Ocean, the 14 Microtops sun-photometer data-sets collected during 
the above-selected MAN cruises were subdivided into three sub-sets 
pertaining to the following large-scale oceanic sectors: (a) Northern 
Greenland Sea and Norwegian Sea (hereinafter referred to as GNS 
sector), between 20° W and 30° E longitudes; (b) Barents Sea and 
West Siberian Sea (BWS sector), between 30° E and 130° E longitudes; 
and (c) the North American Arctic Ocean (NAA sector), including the 
East Chukchi Sea, the Bering Strait, the Beaufort Sea and the Amundsen 
Gulf, between 170° W and 110° W longitudes. No MAN cruises were 
carried out in the East Siberian Sea and West Chukchi Sea sectors, at 
latitudes higher than 67° N and longitudes ranging from 130° E to 
170° W, as well as in the Canadian sector from 110° W to 20° W 
longitudes, including the Baffin Bay, the Davis Straits and the Northern 
Atlantic Ocean. 

The same criteria as used for the ground-based stationary sun- 
photometer measurements taken at the various Arctic sites were 
also adopted to examine the ship-borne measurements, including the 
selection of data made to reject all the measurements affected by strato- 
spheric volcanic particle extinction. The daily mean values of r(0.50 pm) 
and a obtained from the data-sets collected within each of the three 
oceanic sectors were subdivided into monthly sub-sets to determine 
the monthly mean values of both parameters, and define their RFHs 
separately for the winter-spring (Arctic haze) and summer-autumn 
(background aerosol) periods. The main purpose of the analysis was 


to provide evidence of the seasonal characteristics of aerosol optical 
parameters in the Arctic Ocean areas far from anthropogenic aerosol 
sources. 

3.1.3. Northern Greenland Sea and Norwegian Sea 

Monthly mean values of r(0.50 pm) and a derived from the 
Microtops measurements performed during nine cruises made in the 
Norwegian Sea, west of Spitsbergen, and one cruise in the northern 
Greenland Sea and Norwegian Sea are shown in Fig. 14a. Values of 
r(0.50 pm) gradually decreased from less than 0.10 in April to nearly 
0.05 in August, with values of a T < 0.03 over the entire period. The 
monthly mean values of a were rather stable from April to August, rang- 
ing from 1.20 to 1.30, with values of o a equal to 0.10 in spring and 0.40 
in summer. These features suggest that the columnar aerosol content 
was mixed in spring, consisting of anthropogenic (Arctic haze) and nat- 
ural sea-salt aerosols, and contained sea-salt aerosol mixed with BFF 
smoke particles in summer. In fact, the aerosol extinction parameters 
measured in spring and summer with the Microtops sun-photometers 
turn out to be very similar to those determined at Ittoqqortoormiit 
and Hornsund from the AERONET sun-photometer measurements, dur- 
ing spring (when polluted air masses are mainly transported from 
Europe) and summer (when polar air masses are transported from the 
Arctic land and oceanic regions towards the Svalbard Archipelago) 
(Rozwadowska and Sobolewski, 2010; Rozwadowska et al„ 2010). 
These variable features can be more clearly seen by considering the 
seasonal RFHs of r(0.50 pm) and a shown in Fig. 14b, derived from 
the measurements taken on 77 measurement days in the Greenland 
Sea and Norwegian Sea (GNS) sector. The RFH of t( 0.50 pm) exhibits 
a large peak over the 0.02-0.20 range, with a mean value equal to 
0.075, and 25th and 75th percentiles equal to 0.054 and 0.094, respec- 
tively, indicating significant variations even in the remote Arctic 
Ocean. The corresponding RFH of a is also rather wide, showing a 
mean value of 1.23, and 25th and 75th percentiles equal to 1.02 and 
1.47, respectively, as a result of the combination of natural maritime 
and anthropogenic/continental aerosol loads. It is worth noticing that 
the GNS average values of r(0.50 pm) and a agree fairly well with 
those determined at Ittoqqortoormiit (Fig. 5) and Hornsund (Fig. 6) 
and exhibit similar relationships between a and r(0.50 pm) to those 
of the above-mentioned AERONET stations (compare Fig. 8a with 
those of Fig. 8b). 

3.3.2. Barents Sea and West Siberian Sea 

Only the RV Polarstern 2012 cruise crossed the Barents Sea and 
West Siberian Sea (BWS) sector, in August 2012, at latitudes varying 
from 80° N to 84° N and longitudes ranging from 28° E to 120° E (see 
Fig. 2). Thus, the monthly mean values of r(0.50 pm) and a were only 
calculated for three measurement days, giving statistically poor values 
of t(0.50 tun) = 0.054 ± 0.010 and a = 1.42 ± 0.33. As noted in 
Fig. 14a, the above value of r(0.50 pm) closely agrees with that mea- 
sured within the GNS sector on the same month, and a is slightly higher 
than that measured in August within the same sector. Such results are 
typical of an oceanic area during relatively calm wind periods. 

3.3.3. North American Arctic Ocean 

Within the North American Arctic (NAA) oceanic sector, including 
the East Chukchi Sea, Bering Strait, Beaufort Sea and Amundsen Gulf, 
only four cruises were considered among those reported in Table 3. 
The time-patterns of the monthly mean values of t( 0.50 pm) and a 
calculated over the March-September period are shown in Fig. 14a. 
Monthly mean value of r(0.50 pm) were about 0.10 in March, peaked 
in April and May at 0.17 and 0.20, respectively, and decreased in the 
summer months to 0.10 in June and 0.05 in September. No pronounced 
variations in a were obtained, with values from 0.90 to 1.20 from March 
to September, and a a varying from 0.20 (in April) to 0.50 (in July). An 
overall number of 58 measurement days was collected in this oceanic 
sector over the 7-month period from March to September, for which 
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Table 3 

List of the 14 cruises undertaken in the Arctic oceanic regions by various managing institutions in different geographical areas (at latitudes >67° N), where the ship-borne Level 2.0 sun- 
photometer measurements of the Maritime Aerosol Network (MAN) programme were conducted from 2003 to 2012, using hand-held Microtops sun-photometers calibrated at the NASA/ 
GSFC calibration facility (Smirnov et al., 2009, 201 1 ). The peak wavelengths of the narrow-band interference filters mounted on the instruments are given, together with the names of the 
sun-photometer measurement P.I.s. 


Cruises and sun-photometer 
managing institutions 

Geographical area 

Overall 
number of 
measurement 
days 

Measurement 

period 

Peak wavelengths (nm) 
of the spectral channels 
used to measure aerosol 
optical thickness t(A) 

Principal investigator 
and references 

RV Oceania 2003 (IOPAS, Sopot, Poland) 

Norwegian Sea, west 
Spitsbergen 69°-79° N, 2°-14° E 

5 

June-July 2003 

440, 500, 675, 870 

T. Zielinsky 
(Tomasi et al., 2007) 

RV Oceania 2006 (IOPAS, Sopot, Poland) 

Greenland Sea and Norwegian 
Sea, 69°-79° N, 15° W-14° E 

6 

June-July 2006 

440, 500, 675, 870 

T. Zielinsky 
(Tomasi et al., 2007) 

RV Oceania 2007 (IOPAS, Sopot, Poland) 

Norwegian Sea, 69°-78° N, 
1°-16°E 

6 

July 15-August 11,2007 

440, 500, 675, 870 

T. Zielinsky 

CCGS Louis St. Laurent 2007 
(Wood Hole Oceanographic Institution, 
Woods Hole, Massachusetts, USA) 

Beaufort Sea, 70°-79° N, 
125°-150° W 

9 

July 28-August 25, 2007 

440, 500, 675, 870 

A. Proshutinsky 

RV Knorr 2008 (PMEL, NOAA, Seattle, 
Washington, USA) 

Norwegian Sea, 70°-71° N, 
19°-31°E 

4 

April 6-10, 2008 

440, 500, 675, 870 

P. K. Quinn 

CCGS Amundsen 2008 (Institut Maurice- 
Lamontagne, Mont-Joli, Quebec, Canada) 

Beaufort Sea, 70°-73° N, 
121°-131° W 

43 

March 16-August 2, 2008 

440, 500, 675, 870 

P. Larouche 

RVJan Mayen 2009 (IMEDEA, Esporles, 
Mallorca, Spain) 

Norwegian Sea, 76°-80° N, 
11°-28°E 

4 

June 18-22, 2009 

440, 500, 675, 870 

C. Duarte 

RV Oceania 2009 (IOPAS, Sopot, Poland) 

Norwegian Sea, 73°-79° N, 
3°-16°E 

11 

June 23-July 13, 2009 

440, 500, 675, 870 

T. Zielinsky 

CCGS Amundsen 2009 (Universite du 
Quebec, Rimouski, Quebec, Canada) 

Beaufort Sea, 70°-72° N, 
127°-135° W 

9 

July 31 -August 25, 2009 

440, 500, 675, 870 

S. Belanger 

RV Oceania 2010 (IOPAS, Sopot, Poland) 

Norwegian Sea, 68°-80° N, 
0°-13° E 

9 

July 6-August 17, 2010 

380, 440, 500, 675, 870 

T. Zielinsky 

RV Oceania 2011 (IOPAS, Sopot, Poland) 

Norwegian Sea, 69°-79° N, 
3°-19°E 

10 

June 17-August 13, 2011 

440, 500, 675, 870 

T. Zielinsky 

USCGC Healy 201 1 (Naval Research 
Laboratory, Monterey, California, USA) 

Eastern Chukci Sea, Beufort 
Sea and Arctic Canadian Ocean, 
71°-79° N, 142°-165°W 

6 

August 18-September 24, 
2011 

440, 500, 675, 870 

E. A. Reid 

RV Oceania 2012 (IOPAS, Sopot, Poland) 

Norwegian Sea, 73°-78° N, 
10°-20° E 

7 

June 27-July 31, 2012 

440, 500, 675, 870 

T. Zielinsky 

RV Polarstern 2012 
(Institute of Environmental Physics, 
University of Bremen, Bremen, Germany) 

Norwegian Sea, Barents Sea and 
West Siberian Sea, 70°-84° N, 
18°-109°E 

5 

August 2-22, 2012 

440, 500, 675, 870 

L. Istomina 


calculations of the RFHs were made, illustrating a statistically robust 
data-set. The RFH of 7(0.50 |im) showed a more dispersed distribution 
curve over the 0.02-0.25 range than that determined for the GNS sector, 
with a mean value of 0.16, which was considerably higher than that 
measured in the GNS sector, and also had wider right-hand wings. The 
RFH of a exhibits a similar shape to that determined in the GNS sector, 
with a mean value close to 1.10. 

The median values and the 25th and 75th percentiles of parameters 
r(0.50 pm) and a were calculated for the sets of daily mean values de- 
rived from the Microtops sun-photometer measurements conducted 
over the GNS, BWS and NAA sectors. The results are shown in Fig. 8b 
for an easier comparison with those derived from the ground-based 
sun-photometer measurements conducted at Arctic sites. The ship- 
borne measurements provided median values of r(0.50 pm) ranging 
from 0.04 to about 0.13, which are therefore comparable with the 
mean background extinction parameter measured in summer at the 
ground-based stations. The ship-borne median values of a varied main- 
ly from 1.10 to 1.50, and were appreciably lower than those measured in 
summer-autumn at the land stations: the scatter plot shown in Fig. 8b 
provides clear evidence of the prevailing oceanic origins of these aerosol 
loads. 

3.2. Aerosol measurements in the Antarctic Ocean 

Eighteen AERONET/MAN cruises were conducted from 2006 to early 
2013 over the Antarctic Oceans. They are listed in Table 4, together with 
the geographical areas covered by the various cruises, the number of 
measurements days (see also Fig. 2), the measurement periods and 
the spectral characteristics of the portable Microtops sun-photometers 
onboard the vessels. To analyse the aerosol optical parameters more 


homogeneously, the sun-photometer data-sets recorded during the 
18 cruises were subdivided into the following four sub-sets: (1) the 
Southern Indian Ocean (IND) sector, between 20° E and 150° E longi- 
tudes; (2) the Southern Pacific Ocean (PAC) sector, between 150° E 
and 75° W longitudes; (3) the Southern Atlantic Ocean (ATL) sector, be- 
tween 50° W and 20° E longitudes; and (4) the oceanic region around 
the Antarctic Peninsula (APE), between 75° and 50° W longitudes. The 
daily mean values of r(0.50 pm) and a measured with the portable 
Microtops sun-photometers calibrated at NASA/GSFC facility (Smirnov 
et al., 2009, 2011) were analysed separately for coastal data (for dis- 
tances smaller than 300 km from the Antarctic coast) and off-shore 
data (for distances of 300 to 1000 km from the coast), according to 
the criteria used by Wilson et al. (2010) to distinguish the data related 
to a significant landmass from those of pure oceanic origin (i.e„ given 
for the large part by maritime aerosols generated by wind-related 
sources). The overall set of Microtops data was then examined following 
the same criteria adopted to analyse the Arctic ship-borne measure- 
ments, separately for the above-mentioned four oceanic sectors. 

3.2.1. Southern Indian Ocean 

As can be seen in Table 4, eleven cruises were conducted in the 
Southern Indian Ocean, during periods from December to March in 
the years from 2005/2006 to 2011/2012, collecting an overall set of 
226 coastal measurement days, and only nine off-shore measurement 
days (see Fig. 2). The time-patterns of the monthly mean values of 
r(0.50 pm) and a are shown in Fig. 15a for coastal and off-shore data. 
AOT values of r(0.50 pm) were rather stable from December to April 
over the coastal areas, varying from 0.010 to 0.025, with cr T = ±0.01 
on average, and were appreciably greater over the off-shore area, 
being close to -0.06 in January, 0.04 in November and December, and 
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Table 4 

List of the 18 cruises undertaken in the Antarctic oceanic regions by various managing institutions in different geographical areas (at latitudes >62° S), where the ship-borne sun- 
photometer measurements of the Maritime Aerosol Network (MAN) programme were conducted from 2005/2006 to 2012/2013 in the austral summer months, using hand-held 
Microtops sun-photometers calibrated at the NASA/GSFC calibration facility (Smirnov et al., 2009, 2011). The peak wavelengths of the narrow-band interference filters mounted on 
the instruments are given, together with the names of the sun-photometer measurement P.I.S. 


Cruises and sun-photometer 

Geographical area 

Overall 

Measurement period 

Peak wavelengths (nm) 

Principal 

managing institution 


number of 


of the spectral channels 

investigator 



measurement 


used for measuring aerosol 




days 


optical thickness r(A) 


RV Akademik Fedorov 2005/2006 

Southern Indian Ocean 

20 

December 20, 2005-January 26, 

340, 440, 500, 675, 870 

S. M. Sakerin 

(IAO-SB-RAS, Tomsk, Russia) 

65°-70° S, 44°-93° E 


2006 



RV Akademik Fedorov 2006/2007 

Southern Indian Ocean 

40 

December 13, 2006-March 4, 2007 

340, 440, 500, 675, 870 

S. M. Sakerin 

(IAO-SB-RAS, Tomsk, Russia) 

65°-69° S, 46°-93° E 





MV SA Agulhas 2007/2008 (Climatology 

Southern Atlantic Ocean, 

15 

December 18, 2007-January 10, 

440, 500, 675, 870 

S. Piketh 

Research Group, University of Witwatersrand, 

68°-70° S, 2° W-4° E 


2008 


(Wilson 

Johannesburg, South Africa) 





et al., 2010) 

RV Akademik Fedorov 2007/2008 

Southern Indian Ocean, 

18 

December 14, 2007-January 4, 

440, 500, 675, 870 

S. M. Sakerin 

(IAO-SB-RAS, Tomsk, Russia) 

66°-69° S, 45° E-95° E 


2008 




Southern Pacific Ocean, 
68°-72° S, 90°-163° W 

14 

January 24-February 17, 2008 




Antarctic Peninsula, 63°- 

3 

February 20-24, 2008 




65° S, 60°-45° W 
Southern Atlantic Ocean, 
65°-70° S, 50°W-15°E 

13 

March 25-April 22, 2008 



RV Akademik Fedorov 2008/2009 

Southern Indian Ocean, 

11 

December 22, 2008-January 6, 

440, 500, 675, 870 

S. M. Sakerin 

(IAO-SB-RAS, Tomsk, Russia) 

66°-69° S, 76°-93° E 


2009; January 25-March 23, 2009 




Southern Pacific Ocean, 
67° S, 161° E 

1 

January 17, 2009 



Vessel RV Astrolabe 2009/2010 

Southern Indian Ocean 

5 

January 5-23, 2010 

440, 500, 675, 870 

Y. Courcoux 

(OPAR, Univ. de Reunion, Reunion, France) 

66°-67° S, 140°-141°E 





RV Akademik Fedorov 2009/2010 

Southern Indian Ocean, 

12 

December 1 5, 2009-January 4, 

440, 500, 675, 870 

S. M. Sakerin 

(IAO-SB-RAS, Tomsk, Russia) 

65°-70° S, 70°-100° E 


2010; 

March 30-April 6, 2010 




Southern Pacific Ocean 

7 

January 19-28, 2010 




65°-75° S, 87°-172° W 
Antarctic Peninsula, 62° S, 
59° W 

1 

February 6, 2010 



Prince Albert II 2010 (Space Physics Laboratory, 

Antarctic Peninsula 62°- 

11 

January 9-February 17, 2010 

440, 500, 675, 870 

G. Milinevsky 

University of Kyiv, Kyiv, Ukraine) 

67° S, 58°-67° W 





NP Almirante Maximiano 2010/2011 (Rio de 

Antarctic Peninsula 62°- 

7 

January 7-March 21, 2011 

440, 500, 675, 870 

H. 

Janeiro State University, Rio de Janeiro, 
Brazil) 

64° S, 56°-61° W 




Evangelista 

RV Akademik Fedorov 2010/201 1 

Southern Indian Ocean 

18 

December 18, 2010-February 2, 

440, 500, 675, 870 

V. F. 

(AARI, St. Petersburg, Russia) 

65°-69° S, 45°-93° E 


2011 


Radionov 

RV Astrolabe 2011 (OPAR Institute, Univ. de la 

Southern Indian Ocean 

24 

January 5-March 1, 2011 

340, 440, 500, 675, 870 

Y. Courcoux 

Reunion, Saint Denis de la Reunion, France) 

67° S, 140° E 





MV SA Agulhas 2011/2012 ( Climatology 

Southern Atlantic Ocean 

17 

December 19, 2011 -February 13, 

340, 440, 500, 675, 870 

S. Broccardo 

Research Group, University of Witwatersrand, 
Johannesburg, South Africa) 

67°-71°S, 0°-9° W 


2012 



NP Almirante Maximiano 201 1/2012 (Rio de 

Antarctic Peninsula, 61 °- 

6 

October 24, 2011 -November 14, 

440, 500, 675, 870 

H. 

Janeiro State University, Rio de Janeiro, Brazil) 

62° S, 55°-59° W 


2011 


Evangelista 

RV Akademik Fedorov 2011/2012 

Southern Indian Ocean 

19 

December 12, 2011 -March 21, 

440, 500, 675, 870 

V. F. 

(AARI, St. Petersburg, Russia) 

66°-69° S, 37°-93° E 


2012 


Radionov 


Southern Atlantic Ocean 

8 

March 22-April 7, 2012 




62°-70° S, 18° W-20° E 
Antarctic Peninsula 62.2° 
S, 58.9° W 

2 

April 21 and 22, 2012 



RV Astrolabe 2012 (OPAR Institute, Univ. de la 

Southern Indian Ocean 

14 

February 20-March 2, 2012 

380, 440, 500, 675, 870 

Y. Courcoux 

Reunion, Saint Denis de la Reunion, France) 

66° 40' S, 140° E 





RV SA Agulhas 11 2012/2013 (Climatology 

Southern Atlantic Ocean 

11 

December 19, 2012-February 5, 

380, 440, 500, 675, 870 

S. Broccardo 

Research Group, University of Witwatersrand, 
Johannesburg, South Africa) 

70°-71°S, l°-8° W 


2013 



RV Akademik Fedorov 2012/2013 

Southern Indian Ocean 

24 

December 13, 2012-March 15, 

440, 500, 675, 870 

V. F. 

(AARI, St. Petersburg, Russia) 

66°-70° S, 46°-93° E 


2013 


Radionov 


Southern Atlantic Ocean 
67°-70° S, 10°-12°E 

3 

March 22-24, 2013 



RV Akademik Treshnikov 2012/2013 

Antarctic Peninsula, 64°- 

13 

February 11 -March 5, 2013 

440, 500, 675, 870 

V. F. 

(AARI, St. Petersburg, Russia) 

70° S, 57°-72° W 




Radionov 


0.02 in February. The coastal values of a varied from 1 .20 to 1 .40 in the 
December-April period, while the off-shore values decreased from 1.20 
to less than 0.60 from November to December, and then slowly in- 
creased to around 0.90 in February. These findings give a measure of 
the increase in r(A) and the variations in a, which are observed passing 
from coastal to off-shore areas as a result of the stronger production of 
sea-salt coarse particles. The daily mean values of t(0.50 pm) and a 


were derived over the Southern Indian Ocean on 226 measurement 
days near the Antarctic coasts, and on only nine days in off-shore 
areas. Therefore, only the RFHs of the coastal daily mean values of 
r(0.50 pm) and a have been analysed: the RFH of r(0.50 pm) shown 
in Fig. 1 5b presents a very narrow curve, with a mean value of 0.024, 
and 25th and 75th percentiles very close to the mean, while the RFH 
of a exhibits a broad curve, with a mean value of 1.20 and 25th and 
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Fig. 14. Part (a): as in Fig. 3, for Microtops sun-photometer measurements of aerosol optical thickness 7(0.50 pm) and exponent a performed during the cruises conducted from 2003 to 
2012 in the GNS (Greenland Sea and Norwegian Sea) sector (solid circles), from 2006 to 201 2 in the BWS (Barents Sea and West Siberian Sea) sector (solid squares), and from 2008 to 201 1 
in the NAA (Eastern Chuckci Sea, Beaufort Sea and Amundsen Gulf) sector (open triangles) (see also Table 3). Part (b): relative frequency histograms of the daily mean values of aerosol 
optical thickness r( 0.50 pm) and exponent a determined from the Microtops measurements conducted from March to September over the GNS (Greenland Sea and Norwegian Sea) sector, 
and the NAA (North American Arctic Ocean) sector, including the Eastern Chuckci Sea, Beaufort Sea and Amundsen Gulf. 


75th percentiles equal to 0.84 and 1.51. These findings clearly indicate 
that a large variability characterises the sea-salt accumulation and 
coarse mode particle concentrations in coastal areas. 

3.2.2. Southern Pacific Ocean 

As reported in Table 4, only three AERONET/MAN cruises were con- 
ducted in the Southern Pacific Ocean over the December-January 
periods of 2007/2008 to 2009/2010, over the area defined in Fig. 2, 
giving an overall number of 20 coastal and 5 off-shore measurement 
days. The monthly mean values of r(0.50 pm) and a varied from 0.03 
to 0.04 near the coasts, and from 0.04 to 0.05 in the off-shore areas, 
while the coastal monthly mean values of a decreased from -1.45 in 
December to -1.00 in January, and the off-shore values from 1.40 to 
0.70. The values of r(0.50 pm) and a result to closely agree with those 
measured at Neumayer, Novolazarevskaya, Mirny, Syowa, and Mario 
Zucchelli, shown in Figs. 9-12. A coastal median value of r(0.50 pm) 
equal to -0.04 was obtained, with the main percentiles differing by 


less than 0.02, while the median value of a was equal to 1.20, with the 
main percentiles differing by about 0.50. The off-shore data-set provid- 
ed a similar value of r( 0.50 pm) and more distant main percentiles, with 
a median value of a = 1.15, and main percentiles differing by more than 
0.40. 

3.2.3. Southern Atlantic Ocean 

A large set of AERONET/MAN sun-photometer measurements was 
collected in the Southern Atlantic Ocean during five cruises conducted 
in 2007/2008, 201 1/2012 and 2012/2013, as reported in Table 4, giving 
63 coastal and only 8 off-shore measurement days. The monthly mean 
values of r(0.50 pm) varied from 0.03 in December to less than 0.02 in 
April, with a T = 0.01 on average, while those of a increased from 
-1.10 in December to 1,50 in February, and then decreased to 1.30 in 
April. The RFH of the coastal r(0.50 pm) values was found to exhibit sim- 
ilar features to those determined in the coastal area of the Antarctic 
Indian Ocean, with a leptokurtic shape having a mean value of 0.025 
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Fig. 15. Part (a): as in Fig. 3, for Microtops sun-photometer measurements of r(0.50 pm) and a performed during the cruises conducted from late 2005 to spring 2013 in the frame of the 
Maritime Aerosol Network (MAN) activities listed in Table 4. The data refer to coastal (open symbols) and off-shore (solid symbols) measurements, carried out in theIND sector (Southern 
Indian Ocean, circles), PAC sector (Southern Pacific Ocean, upward triangles), ATL sector (Southern Atlantic Ocean, squares), and APE sector (Antarctic Peninsula, downward triangles). 
Part (b): Relative frequency histograms of r(0.50 pm) and a obtained for the coastal data-sets collected (i) on 226 measurement days from November to February over the Southern 
Indian Ocean (IND) sector; (ii) 63 measurement days from December to April over the Southern Atlantic Ocean (ATL) sector; and (iii) on 49 measurement days from October to April 
over the coastal area surrounding the Antarctic Peninsula (APE sector). 


and the 25th and 75th percentiles differing by less than 0.01, while the 
RFH of a showed a rather broad shape over the 0.5-2.4 range, with a 
mean value of 1.39 and the 25th and 75th percentiles differing by 
about 0.30, thus being similar to that determined for the Pacific Ocean 
coastal data-set. 

3.2.4. Around the Antarctic Peninsula 

Numerous aerosol optical data were collected in coastal areas sur- 
rounding the Antarctic Peninsula during the seven AERONET/MAN 
cruises conducted from 2007/2008 to 2012/2013, as reported in 
Table 4, collecting an overall number of 49 measurement days, not 
far from the coasts. The time-patterns of monthly mean values of 
r(0.50 pm) and a obtained from these measurements are shown in 
Fig. 15b. A very high value of r(0.50 pm) close to 0.12 was obtained in 
October, followed by decreasing monthly mean values equal to 0.08 in 
November and lower than 0.06 in December, and then varying from 
0.04 to 0.06 over the January-March period, with a T = 0.02 on average. 
The monthly mean values of a were all rather low, vaiying from 0.50 to 
1.00 over the October-February period, with u a = 0.20 on average. The 


exceptionally high monthly mean values of r(0.50 pm) were deter- 
mined in October and November for measurement days giving values 
of a varying from 0.24 to 0.69, and presenting prevailing transport of 
air masses from the off-shore areas of the Drake Passage. Similar aerosol 
optical characteristics were observed by Posyniak and Markowicz 
(2009) for flows from the North-Eastern quadrant, accompanied by 
low visibility conditions because of the transport of large amounts 
of sea-salt particles from the Drake Passage, with wind speed >10 m/s. 
Relatively low daily mean values of a were measured on these days, 
ranging from 0.43 to 0.73, due to maritime aerosols generated by the 
strong winds, which dominated extinction. For southerly wind circula- 
tion, considerably lower values of r(0.50 pm) were measured in this 
coastal area. Therefore, because of the large variability of wind direc- 
tions, the RFH of r(0.50 pm) determined over the October-March period 
exhibited rather dispersed features, with a mean value of -0.05 and the 
25th and 75th percentiles differing by less than 0.02. These evaluations 
are higher by about a factor 2 than those obtained over the coastal areas 
of the Indian and Atlantic Oceans, because of the more significant contri- 
bution of sea-salt accumulation and coarse mode particles generated by 
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stronger winds. The RFH of a shows a rather broad shape over the range 
a < 1 .5, with a mean value close to 0.70, which is about half of that mea- 
sured over the Indian and Atlantic Ocean coastal areas, and the 25 th and 
75th percentiles differing by no more than 0.30, due to the predominant 
extinction by maritime particles. 

The above results obtained from the coastal and off-shore data 
recorded over the four sectors of the Southern Antarctic Ocean are 
characterised by features of t(0.50 pm) and a typical of sea-salt mari- 
time particles, as can be clearly seen in Fig. 13b, in which the median 
values of a are plotted versus the median values of r(0.50 pm), sepa- 
rately for the coastal and off-shore data. The comparison between 
these results and the scatter plot shown in Fig. 13a shows that they 
are very similar to those derived from the sun-photometer measure- 
ments conducted at the coastal Antarctic stations, since the cluster 
derived from the AERONET/MAN measurements of a and t(0.50 tun) 
essentially covers the same domain of the ground-based measurements 
recorded at Marambio, Novolazarevskaya, Syowa, Neumayer, and Mario 
Zucchelli, although exhibiting slightly lower values of a. 

4. Aerosol backscattering coefficient profiles from 
lidar measurements 

Lidar sends a light pulse through the atmosphere and the telescope 
collects the backscattered lidar return signal. As the speed of light is 
known, one can easily calculate the exact atmospheric position, from 
which the lidar signal has originated. Therefore, the determination of 
aerosol vertical structure is possible, at least over the altitude range in 
which the image of the laser is completely within the field-of-view of 
the recording telescope. Lidar pulses can be generated at several wave- 
lengths and polarisation elements can also be used at the point of light 
entrance to the telescope. This enables a determination of the shape of 
scatterers because spherical particles do not depolarise light, while 
non-spherical particles normally do. 

Overcoming the considerable operational and logistic difficulties, 
various lidar systems have been deployed in the polar regions to 
measure: (i) the vertical profiles of backscatter and extinction coeffi- 
cients of the various aerosols, cloud droplets and ice-crystals, as well 
as their microphysical and radiative characteristics; (ii) the vertical dis- 
tribution of temperature and water vapour mixing ratio; and (iii) the 
ozone concentration at stratospheric levels, to more thoroughly investi- 
gate the role of Polar Stratospheric Clouds in destroying ozone molecules 
at polar latitudes. Actually, most of these measurements were conducted 
to study the dynamic features of the ozone hole in Antarctica and ozone 
depletion in the Arctic region associated with the polar vortex 
occurrences, and the physical characteristics of tropospheric clouds, 
while relatively few measurement campaigns were specifically conduct- 
ed with lidar techniques to investigate the optical characteristics of 
tropospheric aerosols at polar latitudes. An exhaustive description of 
the lidar measurements conducted in polar regions over the past few 
decades was recently made (Nott and Duck, 2011). Lidar activities 
were conducted in the Arctic at the following sites: 

(i) Barrow (Northern Alaska), where long-term lidar measurements 
were recorded and a micropulse lidar is currently used to carry 
out cloud climatology studies that have evaluated average sea- 
sonal cloud occurrences of 65% in winter, 68% in spring, 83% in 
summer, and 89% in autumn, with an annual mean decrease of 
— 4.8% per year over the past 10 years (Dong et al„ 2010). 

(ii) Eureka (Nunavut, Canada), where an elastic lidar was used from 
1993 to 1997 to carry out winter-time aerosol measurements for 
studying the occurrences of clouds and Arctic haze. Backscatter 
peaks due to haze particle layers were frequently observed at 
altitudes lower than 3 km, and more occasionally at altitudes of 
3-5 km, for relative humidity <60% over ice, while clouds domi- 
nate at >80% over ice. Spectral measurements of the aerosol 
backscatter coefficient fibsW and depolarisation ratio tv(A) 


were regularly conducted from 2005 to 2010, over the height 
range from 70 m up to approximately 15 km, using a High Spec- 
tral Resolution Lidar (HSRL) of the University of Wisconsin 
(USA), at the wavelengths 0.532 and 1.064 pm, for the aerosol 
backscatter coefficient /3 bs (A.), and 0.532 pm only for 6 V (A.) and 
PextW. This lidar system was used together with a millimetre 
cloud radar to classify the various cloud particle types. Since 
2008, a troposphere ozone DIAL (Differential Absorption Lidar) 
and a Raman lidar have been operated simultaneously, equipped 
with aerosol channels in the visible and ultraviolet light and 
water vapour and rotational Raman temperature channels. 

(iii) Alert (Nunavut, Canada), where an elastic lidar was operated in 
1 984- 1 986 to measure the vertical profiles of aerosol backscatter 
coefficient. 

(iv) Summit (Central Greenland), where a depolarisation lidar and a 
micropulse lidar were used since 2010, to characterise the aero- 
sol types suspended over the Greenland ice sheet. 

(v) Ny-Alesund (Spitsbergen, Svalbard), where various lidar models 
were used, such as: (a) a stratospheric lidar since 1988; (b) the 
KARL (Koldewey-Aerosol-Raman Lidar) since 1999, this system 
being equipped with aerosol channels, depolarisation and ultravi- 
olet and visible water vapour channels, and subsequently rebuilt 
in 2008 with enhanced multi-wavelength aerosol channels to 
cover the 0.45-30 km altitude range; and (c) an automated 
micropulse lidar (Hoffmann et al., 2009) since 2003, as a part of 
the MPLNET and NDACC networks, which provided cloud and 
aerosol measurements in March and April 2007, useful to charac- 
terise air masses as a function of depolarisation and backscatter 
ratios. 

(vi) Hornsund (Spitsbergen, Svalbard), where an automated lidar has 
been operated since 2009, to measure the vertical profiles of vol- 
ume extinction coefficient fi ex[ ( 0.532 pm) produced by aerosols 
during transport episodes of unpolluted air masses from the 
Greenland Sea and Norwegian Sea areas; and 

(vi) Andoya Rocket Range (Arctic Lidar Observatory for Middle 
Atmosphere Research, ALOMAR), near Andenes (Northern 
Norway), where a tropospheric lidar was installed in 2005 and 
a polarisation-sensitive bistatic lidar system has been used 
since 2006 to study the vertical distribution features of polar 
aerosols and their optical characteristics. 

The above-mentioned field measurements provided vertical profiles 
of the most significant aerosol scattering parameters, illustrating that 
the Arctic haze particles are in general present below the 3 km level dur- 
ing severe haze events, and occasionally at 3-5 km altitudes, as shown 
by the measurements conducted at Eureka (Nott and Duck, 201 1 ). Dur- 
ing the dense haze episode observed with the KARL lidar at Ny-Alesund 
on May 2, 2006, strong extinction features were recorded from the 
ground up to 2.5 km altitude, which showed multi-layered profiles of 
the aerosol backscattering coefficient fi bs (0.532 tun) decreasing on 
average from 3.5 km -1 near the ground to below 0.5 km -1 at 2.5 km 
altitude. Conversely, for unpolluted air conditions, the scale height of 
/3 bs (0.532 pm) was often found to vary from 1.0 to 1.3 km at Ny-Alesund, 
as estimated for instance by Hoffmann et al. (2012) who examined the 
backscatter ratio measurements conducted with the KARL lidar over the 
whole of 2007. 

However, so far our knowledge of Arctic haze from a lidar point of 
view is still incomplete, because the majority of published results are 
based on case studies. Therefore, the yearly cycle of aerosols has been 
studied in this paper on the basis of the lidar measurements conducted 
at Ny-Alesund over the period from 1 November, 2012, to 31 October, 
2013, examining the overall set of KARL lidar measurements taken 
without cloud interferences and with resolutions of 60 m in height 
and 10 min in time, according to Ansmann et al. (1992). To guarantee 
a homogeneous data-set, only lidar profiles with identical technical set- 
tings have been considered. Moreover, cloud screening was conducted 
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and the lowest 800 m were removed from the analysis due to overlap 
(Weitkamp, 2005). Finally, the data derived from the Vaisala RS-92 
radiosonde data collected at the site have been used to subtract the 
Rayleigh backscatter coefficient from the lidar data. The fact that 
the time of the lidar observation does not coincide with the launch of 
the radiosonde should cause errors no higher than 5% in evaluating 
the aerosol backscatter coefficient. This is probably even an upper 
bound, since Rogers et al. (2011) only estimated a 3% bias for the 
CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) lidar on 
board the CALIPSO satellite. Fig. 1 6 shows the monthly averaged profiles 
of the aerosol volumetric backscatter coefficient /3 bs (0.532 pm) from 
January to May and the three bi-monthly averaged profiles from June 
to December, providing evidence of significant variations in this optical 
parameter throughout the year. A clear annual cycle of the aerosol back- 
scatter coefficient can be seen, which increased in February mainly in 
the low troposphere and reached its maximum values during March 
and April in the whole troposphere, subsequently decreasing apprecia- 
bly in May. The months June to January were relatively clear, always 
showing appreciably lower values at all altitudes. Another view of the 
same data is shown in Fig. 17, where the monthly averaged values of 
fi bs ( 0.532 pm) were integrated over five partial height-intervals and 
the total range from 0.8 to 7.0 km. The results show that the Arctic 
haze layer starts to form at low altitudes early in the season and lasts 
until May not only in the low troposphere but also at high altitudes. 
In fact, the maximum values below 2.5 km altitude have been found 
in March, while the largest values of fi bs { 0.532 pm) occur in April at 
1. 5-2.5 km altitudes. It can be also noted in Fig. 17 that the haze season 
starts more rapidly than it disappears, since the increase in backscatter 
observed between January and March is steeper, whereas the decline 
lasts from April till August. 

The annual cycle of the lidar ratio S a (0.532 pm) between the aerosol 
extinction coefficient fi ext ( 0.532 pm) and the aerosol backscatter coeffi- 
cient j3i, s (0.532 pm) is presented in Fig. 18, separately calculated over 
the altitude sub-ranges z < 3.5 km and z > 3.5 km and over the whole 
altitude range, for the KARL lidar measurements conducted from early 
November 2012 to late October 2013. It shows that S a ( 0.532 pm) takes 
the highest values in June-July and generally increases with altitude. 
However, its annual cycle is not very pronounced below 3.5 km altitude, 
where most of the aerosol is located, presenting values varying from 30 
to 35 sr. This is typical of clean continental aerosol (Winker et al., 2009). 
Conversely, the monthly mean values of S a (0.532 pm) calculated above 



Fig. 16. Monthly and bi-monthly averaged vertical profiles of aerosol volume backscatter 
coefficient j3j, s (0.532 pm) obtained from the KARL lidar measurements conducted at 
Ny-Alesund (AWIPEV station) from 1 November, 2012, to 31 October, 2013. 



Fig. 17. Time-patterns of the monthly average values of the aerosol volume backscatter 
coefficient /3 bs (0.532 pm) integrated over the different altitude ranges reported in the 
legend, as obtained from the KARL lidar measurements performed at Ny-Alesund 
(AWIPEV station) from 1 November, 2012, to 31 October, 2013. 


3.5 km were estimated to vary from 30 to 40 sr in winter, i.e. assuming 
values typical of clean continental aerosol and dust, and to increase until 
reaching surprisingly high values of 60-70 sr in summer, which are nor- 
mally typical of BFF smoke and polluted continental aerosol. However, 
one has to keep in mind that the lidar ratio depends not only on the 
chemical composition but also on the shape and size-distribution of 
aerosols. It also agrees well with the evaluations of S a (0.532 pm) made 
by Ackermann (1998) for different aerosol types suspended in air for 
relative humidity conditions ranging from 40% to 50%, estimated to be 
of 27 sr for maritime aerosols and around 55 sr for continental aerosol 
suspended in relatively dry air masses, like those frequently observed 
in the polar atmosphere over land (which probably indicate a decrease 
in the influence of marine air masses with altitude). Due to the contribu- 
tions of both low and high tropospheric regions, the overall values of 
S„(0.532 pm) were estimated to range mainly from 30 to 50 sr. 



Fig. 18. Time-patterns of the monthly average values of the lidar ratio S„(0.532 pm) = 
j3 e *t( 0.532 pm)/Pi, s (0.532 pm) calculated over the whole altitude range (triangles) and 
the altitude sub-ranges z < 3.5 km (squares) and z > 3.5 km (+), as obtained from the 
KARL lidar measurements conducted at Ny-Alesund (AWIPEV station) from 1 November, 
2012, to 31 October, 2013. 
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To quantify the direct radiative forcing of Arctic haze, it is very im- 
portant to know the morphological and optical features of atmospheric 
aerosols, as well as the multimodal size-distribution, the refractive 
index and the particle shape. The ratio of the perpendicular to parallel 
polarised backscatter returns from aerosols at a certain wavelength A 
gives the volume depolarisation ratio 6y(\), which provides an effective 
range-resolved method of determining whether the laser pulse has 
been backscattered by spherical or non-spherical aerosols (Winker 
et al„ 2009). Bearing in mind that spherical particles do not depolarise 
the incoming solar radiation, the scatter plot of the monthly (from 
January to May) and bi-monthly (from June to December) averaged 
values of depolarisation ratio percentage are shown in Fig. 19 versus 
the aerosol backscatter coefficient j3t, s (0.532 pm), as obtained from the 
KARL lidar measurements conducted from early November 2012 to 
late October 2013. Generally quite low values of aerosol depolarisation 
ratio 5y(0.532 pm) have been found (appreciably lower for instance 
than those of desert dust), with the lowest values injanuary, and values 
ranging from 2.0% to 2.6% in October-December, over the range 
of j8 bs (0.532 pm) < 4.5 x 10" 4 km -1 sr" 1 . Remarkably higher 
depolarisation values, greater than 2.6%, have been found during the 
maximum of the haze season, in March and April, for values of 
/3 ta (0.532 pm) ranging from 4.0 x 10" 4 to 5.5 x 10" 4 km -1 sr" 1 . This 
means that the Mie theory, while being a fair assumption to model 
the forcing of Arctic aerosol during most of the year, is less accurate 
during spring time. 

Lidar measurements were taken at the coastal Antarctic stations of 
Syowa, Mario Zucchelli, McMurdo and Dumont d'Urville (66° 40' S, 
140° 01' E), and at the Antarctic Plateau sites of Dome Concordia and 
South Pole over the past decades (see the Antarctic map in Fig. 2), to 
monitor the microphysical parameters of clouds at various altitudes 
(for cloud climatology studies), and only occasionally to measure the 
vertical profile of the tropospheric aerosol scattering coefficient 
/3 ext ( 0.532 pm). It was found in general that the cloud-free vertical pro- 
file of aerosol scattering does not exhibit particularly dense layers near 
the surface, but shows that fi ext ( 0.532 pm) decreases rapidly with 
height, to reach the free troposphere “background aerosol” conditions 
immediately above the marine PBL. The first lidar measurements were 
conducted during the 1974/1975 austral summer at South Pole, where 
regular measurements have also been carried out with a micropulse 
system since 1999, mainly to study the microphysical characteristics 



BACKSCATTER @532nm (10' 3 km' 1 sr‘ 1 ) 


Fig. 19. Scatter plots of the monthly and bi-monthly averaged values of the depolarisation 
ratio (%) versus the aerosol backscatter coefficient /VfO.532 pm), as obtained from the 
KARL lidar measurements conducted at Ny-Alesund (AWIPEV station) from 1 November, 
2012, to 31 October, 2013. 


of diamond dust and blowing snow events. Ground-based lidar mea- 
surements have also been regularly conducted at Dome Concordia 
since 2007, by using the elastic-backscatter and depolarisation lidar 
system of the IFAC-CNR Institute (Florence, Italy). Analysis of these 
data revealed that “diamond dust” ice-crystals are often present during 
windy conditions at this remote site, within the boundary layer of 1 00- 
200 m depth, while aerosols contribute to yield slowly decreasing 
values of /3 sca (0.532 pm) with height, until rapidly reaching the back- 
ground conditions of the free troposphere. 

Diamond dust episodes were also frequently observed over the 
Arctic Ocean. Ground-based lidar and radar were used for this purpose 
during the Surface Heat Budget of the Arctic Ocean (SHEBA) pro- 
gramme, determining the physical characteristics of diamond dust ice- 
crystals and assessing the surface radiative effects induced by such par- 
ticles under cloud-free sky conditions. Examining a set of 188 diamond 
dust or ice-crystals episodes over the western Arctic Ocean between 
November 1997 and May 1998, Intrieri and Shupe (2004) found that di- 
amond dust episodes covered about 13% of the time between November 
and mid-May over the Arctic Ocean and were never observed from mid- 
May to October. Lidar measurements highlighted that the diamond dust 
vertical depth varied from 100 m to 1000 m on the various cloud-free 
sky days, although it was most frequently observed to involve the 
lower 250 m of the troposphere on average, thus contributing to induce 
only small radiative effects on the sea ice surface. 

Airborne lidar measurements of vertical profiles of aerosol extinc- 
tion and backscatter coefficients were conducted in the Arctic, over 
the Svalbard area in spring 2000, 2004 and 2007, during the ASTAR 
(Arctic Study of Aerosol Clouds and Radiation) campaigns, together 
with complementary ground-based lidar measurements at Ny-Alesund 
(Spitsbergen). The ASTAR 2000 campaign ran from 12 March until 25 
April 2000 with extensive flight operations over the Svalbard region 
using the AWI aircraft POLAR 4. Simultaneous ground-based measure- 
ments were conducted at Ny-Alesund (at both German Koldewey and 
Japanese Rabben stations as well as at the Zeppelin station). Vertical 
profiles of various aerosol parameters were measured during the cam- 
paign and provided evidence of the strong temporal variability of 
the Arctic spring aerosol with height. A strong haze event occurred be- 
tween 21 and 25 March, in which AOT measured from ground-based 
observation was 0.18, i.e. largely greater than the aerosol background 
value measured in summer. The airborne measurements made on 23 
March showed a high aerosol layer with an extinction coefficient 
of 0.03 km" 1 or more up to 3 km. The chemical analyses of airborne 
measurements showed that such an aerosol transported from anthro- 
pogenic sources mainly consisted of sulphate, soot and sea-salt particles 
(Yamanouchi et al., 2005). Lidar measurements of 532/355 nm colour 
ratio were obtained by Lampert et al. (2009) during the ASTAR 2007 ex- 
periment, indicating the presence of particles with an effective diameter 
<5 pm. However, Lampert et al. (2010) found rather low Arctic haze 
levels in spring 2007, examining lidar measurements of cloud and 
aerosol layers. As pointed out by Hoffmann et al. (2009), the Arctic 
haze optical thickness t( 0.50 pm) measured during such a seasonal 
period was evaluated to vary mainly from 0.05 to 0.08, compared to a 
14-year average value of 0.10. The PBL extended up to 2.5 km altitude 
and predominantly contained well-hydrated particles, such as local 
sea-spray derived sulphates. 

Various airborne lidar measurements were conducted as part of the 
POLARCAT programme during the 2007-2008 International Polar Year, 
to define the aerosol radiative characteristics. In particular, analysing 
the backscatter measurements provided by the LEANDRE multi- 
channel lidar mounted on the ATR-42 aircraft, cases of aerosol transport 
from Europe and Eastern Asia were monitored in April 2008 by de 
Villiers et al. (2010). These studies showed that anthropogenic aerosols 
originating in Europe exhibited in general smaller sizes than Asian 
particles within the PBL below a height of 0.8 km and within layers at 
heights ranging from 2.8 to 4.5 km, consisting mainly of biomass burn- 
ing substances and urban/industrial particulate matter, with aerosol 
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backscatter coefficient /3 te (0.532 pm) ranging from 1.2 x 10" 3 to 2.2 x 
10“ 3 km" 1 sr" 1 , and depolarisation ratio varying from 1.6% to 2.1% 
(with standard deviations of 0.3% on average). Conversely, the Asian 
dust plume observed in April 2008 was estimated to contain larger 
particles, presenting some well-defined aerosol layers at altitudes 
from 1.7 to 5.6 km, with /3(, s (0.532 pm) varying from 1.6 x 10" 3 to 
3.9 x 10“ 3 km" 1 sr" 1 , and depolarisation ratio values equal to 2.0% 
(smoke particles), 3. 0-5.0% (clean continental aerosol), 8.0% (dust), 
and 4.0-10.0% (polluted dust). 

The above results were obtained using also the Level 2.0 products 
given by the CALIOP lidar mounted on the CAL1PSO satellite, which 
was launched in April 2006 (Winker et al., 2009, 2010) as a part of the 
A-train constellation of satellites, and has an orbit inclination of 98.2°, 
thus providing coverage to 82° latitude. This wide coverage has proven 
immensely important in polar regions, where ground-based installa- 
tions are few and far between. The CALIOP lidar operates at wavelengths 
of 0.532 and 1.064 pm, with a depolarisation ratio measured at 
0.532 pm, and was used to characterise the various aerosol types by as- 
suming in such an approach a set of values of lidar ratios S a ( 0.532 pm) 
and S„( 1.064 tun) equal to 20 and 45 sr for clean marine aerosol, 
35 and 30 sr for clean continental aerosol, 40 and 30 sr for dust, 65 
and 30 sr for polluted dust, 70 and 30 sr for polluted continental aerosol, 
and 70 and 40 sr for combustion smoke particles, respectively (Winker 
et al., 2009). 

A first airborne campaign referred to as the Polar Airborne Mea- 
surements and Arctic Regional Climate Model Simulation Project 
(PAMARCMiP) took place during April 2009 in the Arctic, using a pair 
of downward-looking lidars, a depolarisation aerosol lidar and a tropo- 
spheric ozone DIAL mounted onboard the Polar 5 research aircraft of 
the Alfred Wegener Institute (AWI) for Polar and Marine Research 
(Germany). On the basis of the evaluations made by Rogers et al. 
(201 1 ), it can be assumed that a relative 3% bias of aerosol backscatter 
coefficient j3 bs (0.532 pm) was made on average, due to a wrong assump- 
tion of the Rayleigh backscatter coefficient used to correct the overall 
lidar backscatter coefficient. A series of flights was conducted from 
Longyearbyen (Svalbard) across the western Arctic to Barrow (Alaska), 
reaching Station Nord (Greenland), Alert, NP-36 (at 88° 40' N latitude), 
Eureka, Resolute Bay, Sachs Harbor, Inuvik, and Fairbanks, while com- 
plementary ground-based lidar measurements were regularly conduct- 
ed at Ny-Alesund, Eureka and Barrow (Stone et al., 2010). Significant 
aerosol backscatter and extinction effects were measured over the en- 
tire region, mainly caused by haze particles with diameters ranging 
from 0.13 to 0.20 tun. The PAMARCMiP lidar measurements showed 
that the aerosol extinction coefficient decreased in general with height 
from the surface to the top-level of the thermal inversion layer, found 
to be on average equal to 4 km, and often presented a secondary peak 
above. A second PAMARCMiP airborne campaign was conducted with 
the Polar 5 aircraft during April 2011, following the route from 
Longyearbyen to Barrow, flying over Station Nord, Alert, Eureka, Reso- 
lute Bay and Inuvik (Herber et al., 2012). During both campaigns, 
black carbon was observed at all altitudes sampled but at relatively 
low concentrations compared with historical values. 

A climatological study of Arctic aerosols was made by Devasthale 
et al. (2011 ) examining the CALIOP measurements recorded from June 
2006 to May 2010. The study showed that the bulk of the aerosol parti- 
cles are confined within the first 1000 m of the troposphere, from about 
65% of cases in winter to 45% in summer, while the rest is mainly 
suspended within the middle troposphere, at altitudes ranging from 3 
to 5 km, especially in spring. The relative occurrences of aerosol types 
in the vertical profiles show that clean continental aerosol is the largest 
contributor in all seasons except in summer, while polluted continental 
aerosols are the second largest contributor to the total number of 
observed aerosol layers in winter and spring, and clean marine aerosol 
is the second largest contributor in summer and autumn. The average 
vertical profile of aerosol extinction coefficient was found to exhibit a 
pronounced peak at heights of about 0.4-0.7 km. Associated with the 


intrusions from mid-latitudes, polluted continental and smoke aerosols 
presented much broader distributions of optical and geometrical thick- 
nesses, appearing to be more often optically thicker and higher up in the 
troposphere. 

5. Airborne and satellite measurements 

5.3. Basic remarks 

Airborne sun-photometry employs ground-based techniques 
adapted to the special requirements (notably solar tracking) of an 
aircraft environment. Relatively frequent measurements can be made 
during ascending or descending profiles to achieve vertically stratified 
AOT measurements, which are sensitive to the aerosol vertical structure 
and can be converted to volume extinction coefficient profiles, which 
are analogous to the aerosol backscatter coefficient profiles derived 
from lidar measurements and shown in Fig. 16 (see Section 4). It 
is also possible to carry out measurements onboard airplanes flying 
above and below pollution plumes. These aircraft measurements offer 
the possibility of surveying large areas in relatively short time-periods, 
which is not possible, e.g., for ship-borne measurements. In particular, 
airborne measurements of this kind were made during the PAMARCMiP 
airborne campaign of April 2009 (Stone et al., 2010), in which various 
airborne lidar systems and a multi-wavelength sun-photometer were 
mounted onboard the AWI Polar 5 research aircraft to perform numer- 
ous measurements of volume aerosol scattering and extinction coeffi- 
cients were conducted over the wide area from Svalbard across the 
western Arctic to Barrow, from the European to the Alaskan Arctic, 
and from sub-Arctic latitudes to near the Pole. The sun-photometer 
measurements were performed using a 8-channel sun-photometer sys- 
tem manufactured by NOAA/GMD (Boulder, Colorado, USA) in coopera- 
tion with ISAC-CNR (Bologna, Italy), and equipped with narrow-band 
visible and near-IR channels centred at wavelengths ranging from 
0.355 to 1.064 pm, while the lidar measurements were conducted 
using a pair of downward-looking lidars, a depolarisation aerosol lidar 
and a tropospheric ozone DIAL. Additional independent in-situ mea- 
surements of particle size-distribution were also made onboard the 
AWI Polar 5 aircraft, and light extinction data were derived from the air- 
borne lidar measurements to investigate the spectral effects produced 
by haze particles. 

The measurements revealed that the spatial variations in r(A) ob- 
served during the most turbid period of the 2009 haze season were 
closely related to the atmospheric circulation patterns regulated by a 
dominant airflow from Eurasian anthropogenic sources. In addition, 
lidar observations showed the presence of elevated aerosol layers, as 
shown at Ny-Alesund by Treffeisen et al. (2007) in early May 2006 dur- 
ing a transport event of biomass burning aerosols from agricultural fires 
in eastern Europe, and by Stock et al. (2012) in March 2008 during a BFF 
smoke pollution event over the European Arctic region. 

Background values of r(0.50 pm) not exceeding 0.06 were deter- 
mined in the most remote areas, while haze values of r(0.50 pm) vary- 
ing from 0.12 to around 0.35 were measured, with the highest values 
found in the Beaufort Sea region towards the end of April. Such values 
of t(A) were anomalously high compared with those measured in the 
previous three years, because of the transport of haze particles from 
the industrial regions in Europe and Northern Asia. Arctic haze particles 
were frequently found to be concentrated within and just above the 
surface-based temperature inversion layer, showing in general bimodal 
size-distribution features, consisting of an accumulation mode of mod- 
erately small (water-soluble) particles and an additional mode mainly 
composed of insoluble coarse particles. In addition, the in-situ sampling 
and optical measurements revealed a marked decrease in the mean par- 
ticle size with increasing altitude from the surface to 4 km, yielding 
values of the Angstrom exponent a(0.41 2-0.675 pm) varying from 
1.40 to 1.70. The airborne measurements showed that black carbon 
(BC) was highest near the North Pole, suggesting long-range transport 
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from combustion sources. However, the BC concentration measure- 
ments performed near the surface were nearly an order of magnitude 
lower than those reported from similar campaigns in the 1980s. En- 
hanced opacity at higher altitudes during the campaign was attributed 
to an accumulation of industrial pollutants in the upper troposphere, 
consisting of residual aged aerosol and soot particles originating from 
coal burning in China, in combination with volcanic aerosol resulting 
from the March-April 2009 eruptions of Mount Redoubt in Alaska and 
perhaps minor contributions from aircraft emissions. The Arctic haze 
particles observed during April 2009 were estimated to have reduced 
the net short-wave irradiance by — 2 to — 5 W m~ 2 , resulting in a slight 
cooling of the surface. Examining the data recorded during the second 
PAMARCMiP airborne campaign in April 2011, and some aerosol mea- 
surements coordinated with satellite flyovers of NASA's CALIPSO mis- 
sion, Herber et al. (2012) found that such satellite-borne data can be 
useful to validate aerosol retrievals from the lidar and sun-photometer 
measurements made onboard the Polar 5 aircraft. 

Even larger areas than those covered by airborne sun-photometer 
campaigns can be monitored by means of satellite observations. The 
retrieval of aerosol properties from satellite-based measurements is an 
ill-posed problem. The satellite retrieved aerosol optical thickness over 
bright surfaces such as snow and ice is often biassed and leads to errors 
when estimating aerosol-radiation interactions (IPCC, 2013). Polar- 
orbit satellite-mounted sensors like MODIS (Moderate Resolution 
Imaging Spectroradiometer) and MERIS (Medium Resolution Imaging 
Spectrometer) utilise multi-spectral information to characterise surface 
properties. The NASA/MODIS “dark-target” (DT) approach developed 
for the retrieval of aerosol properties uses the 2.1 pm band to estimate 
the reflectance in visible bands (Kaufman et al., 1997). The empirical 
relationships between the surface reflectance of visible channels and 
2.1 pm channel were improved by considering the effect of geometry 
and surface types (Levy et al„ 2007), while information on global aero- 
sol properties was determined as a function of location and season by 
performing cluster analysis of in-situ measurements (Levy et al., 
2010). The major limitation of the MODIS DT algorithm is that no 
retrievals are performed when the surface reflectance of the 2.1 pm 
channel is higher than 0.25 (Levy et al., 2007). For the retrieval over 
bright surfaces, such as desert and urban regions, Hsu et al (2004) de- 
veloped the Deep Blue (DB) algorithm utilizing the fact that the surface 
is much darker in the short blue spectral channels compared with lon- 
ger wavelengths. The second generation of DB algorithm, the so-called 
Enhanced Deep Blue algorithm, improved estimates of the surface 
reflectance, aerosol model selection and cloud-screening schemes 
(Sayer et al., 2012; Hsu et al., 2013). The current MODIS Collection 6.0 
AOT product is created from three separate retrieval algorithms for 
different surface types (Levy et al., 2013): they are two DT algorithms 
for dark ocean surface as well as the vegetated/dark-soiled land and 
a DB algorithm over desert/arid land. MODIS Collection 6 also provides 
a 3 km AOT product compared to Collection 5 (Remer et al., 2013). 
The Bremen AErosol Retrieval (BAER) algorithm for MERIS (von 
Hoyningen-Huene et al., 2003, 2011) utilises the Normalized Difference 
Vegetation Index (NDVI). Polar-orbit satellites with multi-view obser- 
vations such as AATSR (Advanced Along Track Scanning Radiometer) 
(Curieret al., 2009) and MISR (Multi-angle Imaging Spectroradiometer) 
(Diner et al., 2005a) make use of the Bidirectional Reflectance Distribu- 
tion Function (hereinafter referred to as BRDF) surface properties 
to constrain the ill-posed inverse problem. A time series method is 
the most popular aerosol retrieval method for geostationary satellite 
based on the assumption that the surface reflectance does not change 
significantly during a short period of time (Knapp et al., 2002; 
Govaerts et al„ 2010; Mei et al., 2012). The retrievals based on 
POLDER (Polarization and Directionality of the Earth's Reflectances) 
(Deuze et al., 2001; Dubovik et al., 2011) intensity and degree of 
polarisation measurements provide more accurate AOT products with 
the advantage that the contribution of land surfaces to the degree of 
polarisation at the TOA is generally smaller in the visible as compared 


to the contribution of the underlying surface to the light intensity as 
measured on a satellite. 

Among the various satellite-based observations, those provided by 
MODIS (Justice et al., 1998) mounted on the Terra and Aqua platforms 
are particularly useful for evaluating the AOT features in the polar re- 
gions, over the ocean and ice-free land surfaces, and studying aerosol 
climatology associated with the sources, transport and sinks of specific 
aerosol types (e.g., sulphates and biomass burning smoke particles). 
The twin MODIS sensors have been flying on the Terra platform since 
2000 and on the Aqua platform since 2002, giving a large data-set of 
aerosol products using a number of algorithms over the past decades 
to retrieve columnar aerosol parameters over sea and land (Kaufman 
et al., 1997; Tanre et al., 1997; Hsu et al„ 2004; Remer et al., 2005; 
Levy et al., 2013; Xue et al., 2014). To define the average seasonal 
maps of r(0.55 pm), we downloaded data-sets of MODIS/Terra and 
MODIS/Aqua monthly Level-3 data (Atmosphere Monthly Global lxl 
Degrees Products MOD08_M3.051 and MYD08_M3.051, for the Terra 
and Aqua platforms, respectively) from the Giovanni website (http:// 
disc.sci.gsfc.nasa.gov/giovanni). These had been collected for cloud- 
free atmospheric conditions over the eight years from January 2005 to 
December 2012, at Arctic latitudes >67° N and Antarctic latitudes 
>62° S. The cloud-screened Arctic AOT data-set was divided into 
seasonal sub-sets of r(0.55 pm) for spring (March to May) and summer 
(June to August), separately for Terra and Aqua observations. The sea- 
sonal average maps of Level-3 AOT r(0.55 pm) derived from MODIS/ 
Aqua and MODIS/Terra are separately shown in Fig. 20, for 1 ° x 1 ° pixels 
over oceans and land areas not covered by snow and ice. The values of 
t( 0.55 pm) in the summer months are mainly lower than 0.15, with 
some areas characterised by average values ranging from 0.15 to 0.30 
over the Central and Western Siberian Sea, and the Chukchi Sea, as 
well as over the Beaufort Sea, presumably due to the transport of BFF 
aerosol from Siberia and North America, respectively. The spring aver- 
age values of t( 0.55 pm) range mainly from 0.10 to 0.25, with peaks of 
more than 0.30 in the Siberian, North American and North European 
sectors, associated with dense Arctic haze transported from the anthro- 
pogenic mid-latitude sources. These findings closely agree with the re- 
sults obtained from the ground-based sun-photometer measurements 
conducted at coastal sites (Fig. 8a) and from the ship-borne AERONET/ 
MAN measurements (Fig. 8b). Fig. 20 also shows that AOT cannot be 
retrieved from the MODIS data using the operational MODIS algorithm 
over Greenland and the North Pole in both spring and summer, and 
over large regions of Northern America and Siberia in spring, because 
of the high reflectance of the surfaces covered by sea ice and snow 
and also because of generally small values of AOT in these areas. 

Reliable maps of exponent a cannot be retrieved from MODIS data 
over both land and ocean surfaces because of the low values of AOT in 
the polar regions and the relevant uncertainties affecting the MODIS 
aerosol products (Mishchenko et al„ 2010; Kahn et al., 201 1 ). In partic- 
ular, Levy et al. (2013) determined the Collection 6 MODIS aerosol prod- 
ucts to retrieve AOT and aerosol size-parameters from MODIS-observed 
spectral reflectance data and found that the Angstrom exponent prod- 
uct over land cannot be reliably used. These findings were confirmed 
by Mielonen et al. (2011) who stated that MODIS data do not provide 
quantitative information about aerosol size and parameter a over 
land. Similarly, as a result of a validation study of Collection 5 MODIS 
Level-2 Aqua and Terra AOT and a products over ocean, Schutgens 
et al. (2013) found that these products exhibit significant biases 
due to wind speed and cloudiness of the observed scene, being sig- 
nificantly affected by AOT and a random errors due to cloud fraction 
contributions. 

On the basis of these results, we have decided to examine here only 
the MODIS AOT data and exclude from the present analysis the MODIS a 
products obtained over snow- and ice-free land and ocean regions, since 
they were presumably affected by considerable uncertainties arising 
from: (i) the presence of cloud-fractions within the 1° x 1° pixels that 
cannot be correctly evaluated, and (ii) the variability of AOT at visible 
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Fig. 20. Arctic maps of the seasonal average Level-3 aerosol optical thickness 7(0.55 pm) derived from MODIS/Aqua (upper part) and MODISTerra (lower part) satellite data recorded from 
2005 to 2012 during the spring (left-hand side) and summer (right-hand side) 3-month periods. 


and near-infrared wavelengths, due to cloud extinction effects, which 
may efficiently contribute to lower exponent a. On the basis of the 
AERONET/MAN measurements, this parameter was evaluated to range 
from less than 1.00 to around 1.75 over the Arctic Ocean sectors 
in spring and summer, and from less than 0.50 to about 2.00 in 
Antarctica, with in general lower values over the off-shore areas far 
from the Antarctic coasts. 

The Antarctic aerosol data-set downloaded from the Giovanni 
website over the latitude range >62° S covers not only the land region 
of the Antarctic continent but also the off-shore ocean areas no further 
than 1000 km from the coasts. Only the MODIS/Aqua and MODIS/ 
Terra data collected from December to February were separately exam- 
ined to determine the maps of average r(0.55 pm), limited in practice to 
the oceanic areas only (Fig. 21). The results indicate that such seasonal 
average values of r(0.55 pm) are lower than 0.10 over all ocean areas 
close to the Antarctic coasts and sometimes can exceed 0.10 over the 
off-shore areas, in close agreement with the values derived from 


ground-based sun-photometer measurements at coastal Antarctic 
sites (Fig. 13a) and from AERONET/MAN measurements (Fig. 13b). 

Fig. 21 also shows that no useful information is available from 
MODIS observations over the interior of Antarctica when using the tra- 
ditional retrieval procedures. The retrieval of r(A) over a bright surface 
is indeed a very difficult task, because it is hard to separate the radiance 
contributions by the reflecting surface and atmospheric aerosol back- 
scattering to the overall radiance observed by a satellite-borne sensor 
at the TOA-level, especially for the very large solar zenith angles typical 
of polar latitudes. Mei et al. (2013a) overcame the above difficulties by 
following a synergetic approach based on the use of both MODIS/Terra 
and MODIS/Aqua data, together with: (i) the a-priori assumption of 
aerosol optical parameters retrieved over snow made with the Aerosol 
Properties Retrieval over Snow (APRS) algorithm, and (ii) an appropri- 
ate model of the BRDF reflectance representing snow-covered surfaces. 
The APRS algorithm was based on the operational bi-angle approach 
proposed by Xue and Cracknell (1995) to retrieve non-absorbing 
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Fig. 21. Antarctic maps of the austral summer Level-3 aerosol optical thickness r( 0.55 pm) derived from the MODIS/Terra (left-hand side) and MODIS/Aqua (right-hand side) satellite data 
recorded over the 2005-2012 period. 


aerosol extinction parameters over land surfaces (Tang et al., 2005; 
Wang et al., 2012), in which particulate absorption was considered, 
and a two-stream approximation was adopted. This new algorithm 
was found reliable by Mei et al. (2013a) by means of an extended com- 
parison between the values of r(0.55 pm) retrieved from MODIS data 
and simultaneous AERONET measurements made at six high-latitude 
Arctic stations (Andenes, Barrow, Ittoqqortoormiit, OPAL, Thule, and 
PEARL) in April 201 0 and April 20 1 1 ( i.e., during periods in which Arctic 
haze frequently occurs). A regression line with slope coefficient equal to 
+ 0.764 was found by applying the least squares method to a set of 70 
MODIS retrieved values of r(0.55 pm) obtained with a 10 km x 10 km 
resolution plotted versus the corresponding AERONET measurements, 


with a regression coefficient equal to +0.81 and a root mean square 
eiTor RMSE = 0.079. In particular, comparing the MODIS retrieved values 
of r(0.55 pm) with the AERONET measured values at the above stations, 
Mei et al. (2013a) found relative discrepancies between them ranging 
from a few percent to more than 30% in cases with r(0.55 pm) > 0.20, pre- 
sumably associated with strong Arctic haze extinction. 

These results demonstrate the considerable potential of the APRS 
algorithm to retrieve r(\) at visible wavelengths over the Arctic, for 
highly reflective snow/ice surfaces and large solar zenith angles, as 
can be seen in the maps of r(0.55 pm) derived from the MODIS/Aqua 
observations made on 29 March and 3 May in 2006 (Fig. 22), using 
the method of Mei et al. (2013a). 


MODIS/AQUA AOT [550nm] 29 March, 2006 MODIS/AQUA AOT [550nm] 3 May,2006 
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Fig. 22. Maps of aerosol optical thickness 7(0.55 pm) derived over the Arctic region from the MODIS/Aqua observations made on 29 March, 2006 (left-hand side) and 3 May, 2006 (right- 
hand side) using the method of Mei et al. (2013a). 
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5.2. Aerosol optical thickness retrievals over snow and ice using 
backscattered solar light 

Of particular relevance for remote sensing applications in polar re- 
gions, are the multiple-view (MISR) and double-view (AATSR) sensors. 
The measurements made with these instruments allow the contribution 
from underlying bright surfaces (e.g., snow and ice) to be removed, as- 
suming that the atmosphere is the same for observations made at differ- 
ent angles, in all cases where the atmosphere is not filled with broken 
clouds. MISR is a sensor built by the Jet Propulsion Laboratory (JPL) 
and hosted on the NASA Terra platform, which was launched on 18 
December 1999, and became operational in February 2000. It was 
designed to measure the intensity of solar radiation reflected by the 
surface-atmosphere system in various directions and spectral bands, 
with the main mission of measuring the intensity of solar radiation 
reflected and absorbed by the Earth. The device is composed of nine 
separate digital cameras, each with four spectral bands (blue, green, 
red, and near-infrared), that view the Earth at nine different angles 
(Diner et al., 1998). The algorithms for AOT retrieval from the MISR ob- 
servations are discussed in Martonchik et al. (1998a,b) and Diner et al. 
(2005a), and further assessed by Martonchik et al. (2002) and Kahn 
et al. (2005, 2010). Examining the MISR monthly Level-3 data (Monthly 
Global 0.5° x 0.5° Aerosol Product, MIL3MAE.004) recorded during the 
2005-2012 period, the regional maps of monthly average AOT were 
determined in the present study, over the Arctic region for latitudes 
>67° N, and the Antarctic region for latitudes >62° S. The monthly 
average values of r(0.55 pm) collected over the Arctic region under 
cloudless sky conditions were separated into a pair of seasonal sub- 
sets for spring (March-April-May) and summer (June-July-August), 
from which the seasonal average maps of t(0.55 pm) shown in Fig. 23 
were obtained. It can be seen that these results are very similar to 
those retrieved from the MODIS data in Fig. 20, with spring average 
values of t(0.55 pm) varying mainly from 0.10 to 0.25, and being higher 
than 0.20 (and sometimes exceeding 0.30) over large areas of Central 
and Eastern Siberia as in the North American sector, presumably associ- 
ated with dense Arctic haze transport episodes. The average summer 
values of r(0.55 pm) were found to mainly range from 0.05 to 0.15, 
showing similar features to those detected by MODIS observations. 
With regard to this, it is worth mentioning that Campbell et al (2012) 
examined (i) MODIS and MISR observations made over the Arctic 
regions in 2007, obtaining estimates of r(0.55 pm) all lower than 
0.10, and (ii) a CALIOP data-set collected in 2007, determining mean 
day-time and night-time values of r(0.532 pm) all lower than 0.15. 
In addition, examining the CALIOP data recorded over the Arctic in 
2008, Winker et al. (2013) obtained mean cloud-free day-time values 
of 7(0.532 pm) mainly lower than 0.05 and mean cloud-free night- 
time values of r(0.532 pm) not exceeding 0.10 in March-May and 
September-February. In particular, they determined average values of 
7(0.532 pm) from the CALIOP observations made from 2007 to 2011 
that were lower than: (i) 0.02 over inner Greenland, (ii) 0.08 over the 
GNS sector surrounding the Svalbard region, (iii) 0.12 over the Scandi- 
navian area, (iv) 0.12 over the Western Siberian Sea sector, (v) 0.05 
over the Eastern Siberian Sea sector, and (vi) 0.06 over the North 
American sector of the Arctic Ocean. These data are in good agreement 
with ground-based sun-photometer measurements (Figs. 3-8) and 
AERONET/MAN ship-borne measurements (Figs. 8b and 14), which 
showed that 7(0.50 pm) mainly varies from 0.08 to 0.16 in spring and 
from 0.04 to 0.10 in summer. 

The same procedure was also followed in the analysis of MISR data 
recorded over Antarctica during the austral summer from 2005 to 
2012 (Fig. 23), which indicates that the seasonal average values of 
7(0.55 pm) are lower than 0.10 over the greater part of the coastal 
ocean areas around Antarctica and vary mainly from 0.10 to 0.25 over 
the off-shore areas very far from the coasts, due to effective sea-salt 
production by strong winds. These features are very similar to those ob- 
tained from MODIS observations (Fig. 21) and those derived from the 


AERONET/MAN measurements carried out in the coastal and off-shore 
ocean areas surrounding the Antarctic continent (Fig. 13b). 

The AATSR sensor is a dual-view (nadir and 55 forward) radiometer 
on board the European Space Agency (ESA) Environmental Satellite 
(ENVISAT), which was launched in March 2002 and is one of the satel- 
lite instruments designed for providing a well-calibrated long-term 
global data-set of satellite data for climate research (Grey and North, 

2009) . ENVISAT was on a sun-synchronous orbit with an equator pass- 
ing time of about 10 am. AATSR has a ground pixel resolution which is 
equal to 1 x 1 km 2 at nadir and a swath width of 500 km, thus taking 
about six days to achieve full global coverage (three days at mid- 
latitudes). Observations in both views are made at seven wavelengths: 
0.55 pm, 0.66 pm, 0.87 pm, 1 .6 pm, 3.7 pm, 1 1 pm and 1 2 pm, the spectral 
resolution of visible channels is approximately 20 nm, which avoids at- 
mospheric water vapour absorption regions in the electromagnetic 
spectrum. Contact was lost with ENVISAT on April 2012 after 10 years 
of service. The follow-up SLSTR mission as currently foreseen on 
Sentinel-3 will ensure the continuity of the multi-view-angle method 
into the future. 

The reflectance at the TOA-level can be described as the sum of 
the atmospheric path reflectance (black surface reflectance) and the 
atmosphere-surface interaction (Kaufman et al., 1997). With the use 
of AATSR dual-view observations, the ratio between observations in 
the forward and nadir directions can be determined by calculating the 
difference between TOA reflectance and atmospheric path reflectance 
in both directions and taking features of the total atmospheric transmis- 
sion into account (Istomina et al., 2010; Mei et al., 2013b). It is worth 
mentioning that the ratio of both AATSR directional observations is 
not sensitive to the snow albedo, which reduces the high snow albedo 
effect on AOT retrieval over the Arctic region. For pure snow, Istomina 
et al. (2009) used an analytical snow BRDF model to estimate the ratio 
of both directional observations with a correction for TOA reflectance, 
where this approach was affected by the shape of the BRDF and not by 
its magnitude (Vermote et al., 1997a). Mei and Xue (2013) used an 
equivalent snow and ice mixture pixel for the spectral surface reflec- 
tance: it is constructed using a linear mixing model of snow and ice 
spectra tuned by the Normalized Difference Snow Index (NDSI) as an 
indicator of snow cover. Before calculating NDSI, cloud and free water 
were masked and a Look-Up-Table method was used to retrieve AOT. 

Current operational AOT retrieval algorithms, such as MODIS DT 
(Levy et al„ 2013), MISR algorithm (Diner et al., 2005b; Kahn et al., 

2010) and AATSR (Holzer-Popp et al., 2013), can determine AOT over 
the Arctic Ocean (Glantz et al., 2014). Evaluation of the MODIS/Aqua 
AOT product over Svalbard using a 9-year set of AERONET observations 
gives an accuracy of ±0.03 ± 0.05 AOT (Glantz et al., 2014). The AOT 
maps defined by Holzer-Popp et al. (2013) in September 2008 qualita- 
tively show an overestimation over the Arctic Ocean compared with 
the 9-year mean values from the MODIS Aqua product, AERONET 
and Alfred Wegener Institute (AWI) observations, which have been 
separately estimated, giving values of AOT equal to ± 0.031 ± 0.021, 
±0.031 ± 0.022, and ±0.024 ± 0.052, respectively (Glantz et al., 
2014). However, due to the high reflectance of surface encountered for 
large solar zenith angles in the Arctic (Mei et al., 2013a; Macdonald 
et al., 2014), it is very difficult to quantitatively retrieve AOT over snow 
and ice. Stamnes et al (2004) pointed out that the TOA reflectance is 
sensitive to AOT over both very bright and dark surfaces, demonstrating 
the possibility of retrieving AOT over Arctic regions. Recent publications 
(Istomina et al., 2011; Mei et al., 2013a, b) show promising AOT results 
over snow/ice-covered regions, as described below. 

Preliminary attempts based on the above concepts were tried by 
Istomina et al. (2010) to individualise cloud-free snow-covered areas 
using AATSR measurements, and to discriminate clear snow fields for 
the retrieval of AOT. In addition to the aerosol retrieval method at the 
above-mentioned visible and near-infrared wavelengths (Istomina 
et al., 2009), the AATSR measurements in the 3.7 pm channel were 
also utilised to retrieve AOT (Istomina et al., 2011). Radiative transfer 
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Fig. 23. Upper part: Arctic maps of the seasonal average Level-3 aerosol optical thickness r(0.55 pm) derived from the MISR satellite data recorded from 2005 to 2012 during the spring 
(left-hand side) and summer (right-hand side) 3-month periods. Lower part: As in the upper part, for the austral summer average Level-3 aerosol optical thickness 7(0.55 pm) derived 
from MISR satellite data recorded from 2005 to 2012 over oceans and land areas not covered by snow and ice. 


simulations for the accumulation and coarse particle modes of four 
main aerosol components were conducted in order to represent the re- 
trieved r(A.) in the visible region of the spectrum. The advantage of this 
algorithm is that it can be used over any blackbody-like surface (open 
ocean, sea ice, snow-covered land). An example of its usage is shown 
in Fig. 24. 

Improved results were obtained by Mei et al. (2013b) using the 
Dual- View Multi-Spectral (DVMS) approach, in which the dual view is 
used to separately evaluate the contributions made by atmospheric 
aerosols and underlying surfaces to the reflectance measured by this 
satellite sensor at the TOA-level. The algorithm uses an analytical 
snow BRDF model to estimate the ratio of snow reflectances in the 
nadir and forward views, and the atmospheric contribution to the 
TOA-level reflectance, obtained using the dark pixel method over 
the adjacent ocean surface, by assuming that this value applies over 
nearby land surfaces in the absence of significant sources across the 


coastline. An iteration involving all four AATSR spectral channels in 
the visible and near-infrared is used to retrieve the most relevant 
t( 0.55 pm) information. The method was illustrated for AATSR over- 
passes over Greenland in April 2009, during cloud-free sky periods. 

Some examples of the results achieved by Mei et al. (2013b) in 
Fig. 25 show the daily maps of r( 0.55 pm) retrieved from AATSR data 
recorded on four days in April 2009 with 1 km x 1 km resolution over 
the western part of Greenland, under cloud-free sky conditions. It can 
be seen that r(0.55 pm) had very low values smaller than 0.10, with 
little spatial variation for moderate contributions of sea spray aerosol 
associated with the very low wind speeds at the surface, as observed 
on the chosen days. 

A comparison test between the values of r(0.55 pm) retrieved 
using the DVMS approach for high-quality AATSR observations and 
the corresponding AERONET values measured at Thule (North-western 
Greenland) during April 2009 exhibited a good correlation, with a 
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Fig. 24. Time sequence of aerosol optical thickness r(0.55 pm) retrieved using AATSR data over sea ice and snow-covered land surfaces with the algorithm described by Istomina et al. 
(2011). Left panel refers to orbit No. 31673 on 21 March, 2008; middle panel to orbit No. 31687 on 22 March, 2008; and right panel to orbit No. 31773 on 28 March, 2008. The increase 
in r(0.55 pm) might be connected to ozone depletion/bromine explosions observed during March 2008 in the region (Nghiem et al., 2012; Moore et al., 2014). 



regression coefficient equal to +0.76 (Mei et al., 2013b). However, 
it should be mentioned that such good results were obtained for a 
selected set of measurements, containing a high number of cases with 
r(0.55 (tm) > 0.12, due to the frequent occurrence of Arctic haze epi- 
sodes in April. A synergetic approach to retrieve aerosol optical charac- 
teristics over the Arctic, using data from the MODIS sensor mounted 
onboard the Terra and Aqua platforms and prior knowledge of aerosol 
optical parameters retrieved over snow (as assumed in the APRS proce- 
dure), was presented by Mei et al. (2013a). Bearing in mind that cloud 
contamination can cause abnormally high retrieved values of AOT, the 
analysis of MODIS data was limited to only clear-sky pixels selected 
through MODIS cloud products and visual inspection. In applying the 
APRS procedure, it was realistically assumed that the aerosol optical 
parameters do not change during the time between the overpasses of 
Terra and Aqua platforms. A comparison was also made between the 
values of r(0.55 pm) retrieved over Greenland from MODIS observa- 
tions following the APRS procedure and those derived from AATSR 
data with the DVMS method. This comparison showed that the APRS 
values of r(0.55 pm) ranged from 0.07 to 0.09 over the monitored 
areas, while the DVMS ones were mainly -0.09, and exhibited a maxi- 
mum value just 0.6% lower than that derived from MODIS data. 

It has to be mentioned that the global application of the above-cited 
AOT retrieval algorithms for both MODIS and AATSR is still very much 
limited due to several challenges which are listed below. For instance, 
one of the major issues is the cloud masking over snow and ice. Cloud 
screening for AOT retrieval in the Arctic is a challenge due to the similar- 
ity of snow/ice and cloud reflective properties in the visible region 
of the spectrum. The sensitivity of AOT retrievals over snow to atmo- 
spheric particulate matter implies a potentially high contamination 
of the aerosol product with eventually unscreened clouds. It is for 
these two reasons that many cloud masking algorithms developed for 
a general case of multiple surface types, such as open water, vegetation 
and various soils, are of questionable quality over snow and ice. An 
example of this is the operational AATSR cloud mask provided with 
the level lb AATSR product (Istomina et al., 2010). The main approach 
to distinguishing clouds from the underlying surface is that clouds are 
generally brighter and colder than the underlying surface. However, 
this is not always the case over snow/ice-covered regions such as the 
Arctic. A prominent spectral signature of snow and ice makes it possible 
to discriminate clouds using NIR and TIR channels. At the same time, not 
all sensors have such channels. Therefore, cloud-screening algorithms 


over snow differ from sensor to sensor in order to most effectively uti- 
lise the available features (which are not only the spectral and spatial 
resolution, but also different observation directions and polarisation 
measurements). The resulting cloud mask over snow and ice obtained 
from different sensors also differs in quality. Employing absolute or 
relative thresholds (e.g. NDSI or NDVI) in the available VIS, NIR and IR 
channels to separate clouds from snow and ice is a commonly used ap- 
proach which gives a robust cloud mask over snow for common applica- 
tions (MERIS, AATSR, SCIAMACHY) (Allen et al., 1990; Istomina et al., 
2010; Schlundt et al., 2011). However, as the thresholds are often em- 
pirically derived and based on a limited data-set, this approach might 
have drawbacks when it comes to particular surface or cloud types, 
e.g. ice clouds over snow. Possible mis-classification here could be: 
( 1 ) mistaking heavy aerosol loads for clouds, as e.g. might be the case 
of MODIS test BT11-BT3.9 (Ackerman et al., 2010), or (2) mistaking 
snow surfaces for clouds as in MODIS test BT11-BT12, which is used 
to detect cirrus over many surface types also covered by snow and ice 
(Ackerman et al., 2010). Mis-classification may occur due to snow emis- 
sivity being sensitive to grain size at these wavelengths, and creating a 
difference between the brightness temperatures BTs recorded in the 
11 and 12 pm channels (Hori et al., 2006). 

A snow BRDF model can also be used to derive the threshold value 
for the NIR spectral channel, as shown by Spangenberg et al. (2001 ) 
and Trepte et al. (2001). This approach might have difficulties in 
distinguishing cloudy conditions from cloud-free but aerosol-polluted 
scenarios. The best results can be obtained by combining the various 
approaches such as thresholding, image clustering, supervised and un- 
supervised classification, usage of oxygen and water vapour absorption 
bands if available (e.g. MERIS) (Gomez-Chovaetal., 2007) together with 
usage of external topography data, land/water masks, ecosystem maps 
and daily operational snow/ice products (MODIS, MISR) (Ackerman 
et al., 1998; Shi et al., 2002, 2007; Liu et al., 2004). The MODIS cloud 
mask product with confidence cloudy flag generally works well for 
thick warm clouds over the Arctic by using 1 9 channels and auxiliary in- 
puts such as the topography, but will have difficulties with optically thin 
cold clouds such as cirrus due to the similarity of both cloud types. 
A comprehensive combined cloud mask is also available for AVHRR 
(Saunders and Kriebel, 1988; Gesell, 1989; Key and Barry, 1989; Han 
et al., 1999; Mei et al., 2014). However, the challenge of applying a 
cloud mask over snow is still current even for radiometers equipped 
with IR and NIR bands. 
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Fig. 25. Maps of aerosol optical thickness t(0.55 pm) retrieved from AATSR data taken with 
Greenland using the method of Mei et al. (2013b). 



km x 1 km resolution on 9, 15, 18, and 21 April, 2009, over snow-covered surfaces in west 


Another challenge of any AOT retrieval over snow is the variability of 
snow and sea ice types in the field, the effect of surface melt or ageing, 
pollution, and regarding sea ice the effect of varying sea ice concentra- 
tions. All these effects make it difficult to reliably account for the surface 
signal contribution to AOT for retrievals over snow surfaces even when 
double-viewing observational techniques (e.g. AATSR or simultaneous 
MODIS Terra and Aqua usage) are used. At the same time, as the AOT 
in the Arctic is rather low (background values -0.05, up to 0.1 -0.3 
during haze events with occasionally higher peaks), the surface signal 
is responsible for the greatest part of the TOA signal measured by satel- 
lites. Multi-spectral AOT retrievals are in addition challenged by a high 
variability in ice and snow surface types, depending on grain size, pollu- 
tion, density, etc. (Warren and Wiscombe, 1980; Negi et al., 2010). The 
variability of the aerosol types, which includes a varying scattering 


phase function shape, spectral behaviour and particle size-distribution, 
makes it even more challenging to reliably retrieve the AOT, especially 
when spectral information on both the surface and atmosphere is 
used simultaneously. The sensitivity to complex polar atmospheric con- 
ditions and temperature inversions (especially for IR channels used for 
the retrieval or the cloud mask) are currently not taken into account 
and also need to be investigated. The satellite-retrieval algorithms 
proposed by Istomina et al. (2011); Mei et al. (2013a) and Mei et al. 
(2013b) can only be used for an initial qualitative evaluation of the 
atmospheric aerosol load and only for selected apparently cloud-free 
scenes. They are able to give a correct impression of spatial and tempo- 
ral AOT distributions on a regional scale under the assumption that the 
surface and aerosol properties do not change, but need to be further im- 
proved with more flexibility regarding surface types, their extensive 
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validation and specifically developed cloud masking being the first pri- 
orities. Due to the above listed issues and challenges, application on a 
global basis is currently not possible. 

6. Optical characteristics of Arctic and Antarctic aerosols 

To obtain a realistic representation of the optical properties of 
atmospheric aerosol derived from ground-based and ship-borne sun- 
photometer measurements, a set of aerosol extinction models has 
been defined here by taking into account the evaluations of aerosol 
radiative parameters determined from sun-photometric and in-situ op- 
tical measurements made at various Arctic and Antarctic sites. The main 
results obtained for the nuclei, accumulation and coarse particles are 
presented in the following two sub-sections for four Arctic and four Ant- 
arctic aerosol types. 

6.1. Arctic aerosol particle size-distributions and optical characteristics 

Multi-year measurements of aerosol chemical composition and light 
scattering and absorption coefficients were conducted at Barrow by 
Quinn et al. (2002), separately for the sub-micron and super-micron 
particle modes, showing that: (i) extinction effects are dominated by 
sulphate fine particles during the spring-time Arctic haze episodes, 
and by sea-salt accumulation mode particles formed from wind- 
driven sea spray in winter, and (ii) sub-micron sulphate and sea-salt 
particles efficiently contribute during summer to attenuate the incom- 
ing solar radiation. Both sulphate and sea-salt particle concentrations 
followed well-defined annual cycles at Barrow, Alert and Arctic EMEP 
sites (Quinn et al., 2007). Bimodal features of the aerosol size- 
distribution were reported at Barrow by Delene and Ogren (2002), 
who found that the overall visible light scattering and absorption at 
the ground-level is given over the entire year by two particle modes: 

(i) an accumulation mode (with mode diameter D c < 1 pm) yielding 
annual mean values of volume scattering coefficient [i sca = 6.17 ± 
3.61 Mm~\ volume absorption coefficient fi abs = 0.36 ± 0.38 Mm -1 , 
ground-level Angstrom exponent a„ = 1.67 ± 0.36, and ground-level 
single scattering albedo co 0 = 0.954 ± 0.028; and (ii) a coarse mode 
(with D c < 10 pm), giving annual mean values of /3 sca = 9.76 ± 
5.20 Mm~\ p abs = 0.39 ± 0.41 Mm -1 , a 0 = 1.11 ± 0.39, and = 
0.965 ± 0.023. 

Multimodal features of the fine particle size-distribution were also 
reported by Strom et al. (2003) on examining particulate matter 
sampled at the Zeppelin station (near Ny-Alesund) with a Differential 
Mobility Particle Sizer (DMPS) over the D < 1 tun diameter range. 
For such multimodal characteristics of the sub-micron aerosols and con- 
sidering that a varies throughout the year from more than 1.8 (for 
prevailing extinction by fine particles) to less than 1 (for predominant 
attenuation by coarse particles), a set of particulate extinction models 
is proposed here, each given by a linear combination of a nuclei mode 
with an accumulation or coarse particle mode, to represent the optical 
properties of airborne aerosols. In this approach, the overall size-distri- 
bution can be defined for a certain value of a by vaiying the percentage 
number density concentrations of both modes until the measured value 
of a is fitted. 

Each aerosol unimodal size-distribution curve was represented by a 
log-normal curve having the analytical form, 

N()(r) = dN(r)/d( Log r) 

= ^ N ° exp [- 1 f L ° g , r — Log T ' )*\ , (3) 

\/2n( In 10)( Log a) L 2 V Log u )\ 

where N 0 is the total particle number concentration (measured in 
cm' 3 ), In 10 is a constant approximately equal to 2.3026, Log is the 
decadal logarithm (with base = 10), cris the geometric standard devia- 
tion, and r c is the mode radius (measured in pm). Thus, the following 


size-distribution curves were defined to represent four different Arctic 
aerosol types: 

(1 ) The size-distribution curve of winter-spring Arctic haze parti- 
cles, consisting of: (i) a nuclei mode of dry-air Arctic water- 
soluble aerosol, and (ii) an accumulation particle mode of sea- 
salt particles, according to the Hess et al. (1998) OPAC models. 

(2) The average size-distribution of summer background aerosols, 
consisting of: (i) a fine particle mode predominantly composed 
of mineral dust nuclei and water-soluble substances, and (ii) a 
coarse particle mode of sea-salt particles, as defined by Hess 
et al. (1998). 

(3) The average size-distribution curve of Asian dust, consisting 
of: (i) a mineral dust nuclei mode, and (ii) a mineral coarse 
particle mode, having features similar to those represented 
by Hess et al. (1998) in the mineral dust OPAC models, and 
chosen to simulate an episode of Asian dust transport over a 
ground-layer of unpolluted aerosol observed by Stone et al. 
(2007) at Barrow. 

(4) The overall size-distribution curve of a summer background 
aerosol model containing BFF smoke, assumed to consist of: 
(i) a combustion nuclei mode, and (ii) an accumulation parti- 
cle mode. Both modes were assumed to consist of a mixture of 
combustion dust and soot particles giving a ratio of 12.86 
between the scattering and absorption coefficients, leading 
to a single scattering albedo value co = 0.928, as indicated by 
the in-situ measurements conducted by Mielonen et al. 
(2013) for Russian wildfire particles. 

Each of the above eight modes was represented in terms of Eq. (3) 
for the values of shape-parameters r c and a given in Table 5, while 
the unimodal optical/chemical characteristics were represented assum- 
ing the mass percentages of the 6S dry-air components defined by 
Vermote et al. (1997a) and provided in Table 5 for each pair of modes, 
according to the chemical composition estimates made by Quinn et al. 
(2007) and Tomasi et al. (2012). 

The spectral values of the real n(A) and imaginary fc(\) parts of 
the particle refractive index obtained for the Arctic aerosol models 
given in Table 5 were found to decrease very slowly with wavelength 
over the 0.40-1.0 pm spectral range. They are reported in Table 5 for 
A = 0.55 tun, together with values of the single scattering albedo 
w(0.55 pm), asymmetry factor g( 0.55 pm), volume extinction coefficient 
/3 ex[ (0.55 pm), and exponent a. It can be seen in Table 5 that the single 
scattering albedo d)(0.55 pm) of the Arctic haze nuclei mode is equal 
to 0.86 and that of the accumulation mode is equal to 0.937. Therefore, 
using the present bimodal model, <u(0.55 pm) is made to vary over the 
0.86-0.94 range as a function of the mass fractions of nuclei and accu- 
mulation particles. The assumptions made for this optical parameter 
agree very well with ground-level measurements of <n o (0.55 pm) ob- 
tained from in-situ nephelometer and aethalometer 2 measurements 
carried out by: (i) Bodhaine (1995), who estimated a monthly mean 
value of 0.928 in March at Barrow, (ii) Sharma et al. (2006), who deter- 
mined an average value of 0.94 from February to May at Barrow, and 
(iii) Strom et al. (2003), who estimated a monthly mean value 
&> o (0.55 pm) = 0.94 in spring at Zeppelin, and monthly mean values 
vaiying from 0.84 to 0.91 in late autumn and winter. 

The calculations presented in Table 5 provide a single scattering 
albedo value co(0.55 pm) = 0.93 for the Arctic summer background 
fine particle mode, and 0.81 for the corresponding coarse particle 
mode, suggesting that w(0.55 pm) varies mainly over the 0.81-0.93 
range during the Arctic summer. The results do not differ considerably 
from those of: (i) Strom et al. (2003), who found that w o (0.55 pm) as- 
sumes values vaiying mainly from 0.94 to 0.98 in spring and summer 


2 for detailed technical characteristics of the two instruments, see Bodhaine (1995) and 

Anderson and Ogren (1998). 
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Table 5 

Values of shape-parameters a and r c of the log-normal curves adopted to represent the nuclei, accumulation and coarse particle modes determined at Arctic and Antarctic sites, all nor- 
malised to give the value of overall particle number concentration N 0 = 1000 cm“ 3 . They are given together with the spectral values of real part n and imaginary part k of the complex 
refractive index at the 0.55 pm wavelength, single scattering albedo io(0.55 pm), asymmetry factor g(0.55 pm), volume extinction coefficient fi ext ( 0.55 pm) and Angstrom's exponent a 
calculated over the 0.40-0.87 tun wavelength range for the spectral evaluations of /W(A) made using the 6S radiative transfer code of Vermote et al. (1997a). The 6S aerosol components 
are labelled according to Vermote et al. (1997 b): oceanic (OC), water-soluble (WS), dust-like (DL) and soot (SO). 


Aerosol type 

Log-normal 

mode 

Mass percentages of the basic 
6S dry-air components 

Geometric 
standard 
deviation a 

Mode 

radius 

r c 

(pm) 

Particulate matter 
refractive index parts 

(O 

(0.55 pm) 

g 

(0.55 pm) 

fiext 

(0.55 pm) 
(km -1 ) 

Exponent 

a 

OC 

WS 

DL 

SO 

n 

(0.55 pm) 

k 

(0.55 pm) 

Winter-spring 

Nuclei 

35 

39 

22 

4 

2.24 

2.1 10“ 2 

1.487 

2.17 10“ 2 

0.864 

0.637 

1.221 10“ 8 

1.581 

(Arctic haze) aerosol 

Accum. 

58 

4 

38 

- 

2.03 

3.0 10“’ 

1.444 

3.28 10“ 3 

0.937 

0.745 

4.088 10“ 3 

-0.117 

Arctic summer 

Fine 

18 

35 

45.8 

1.2 

1.95 

3.5 10“ 2 

1.506 

1.10 10“ 2 

0.930 

0.605 

4.353 10“ 8 

1.871 

background aerosol 

Coarse 

71 

1 

28 

- 

2.03 

1.75 

1.424 

2.30 10“ 3 

0.813 

0.846 

3.200 

-0.054 

Asian dust 

Fine 

- 

24 

76 

- 

1.95 

7.0 10“ 2 

1.530 

7.52 10“ 3 

0.956 

0.666 

5.094 10“ 6 

1.039 


Coarse 

4 

6 

90 

- 

2.15 

1.30 

1.552 

7.56 10“ 3 

0.673 

0.860 

1.385 

-0.058 

Boreal forest fire smoke 

Nuclei 

- 

- 

97.9 

2.1 

2.00 

3.9 10“ 2 

1.535 

1.71 10“ 2 

0.906 

0.627 

1.698 10“ 7 

1.556 


Accum. 

- 

- 

98.4 

1.6 

2.00 

1.210“’ 

1.534 

1.49 10“ 2 

0.892 

0.700 

9.878 10“ 5 

0.345 

Antarctic austral summer 

Fine 

52 

40 

7.5 

0.5 

2.03 

7.0 10“ 2 

1.454 

6.30 10“ 3 

0.958 

0.712 

5.811 10“ 6 

0.997 

coastal aerosol 

Coarse 

61 

16.5 

22.5 

- 

2.10 

1.50 

1.439 

2.79 10“ 3 

0.795 

0.844 

2.117 

-0.058 

Austral summer 

Fine 

11.2 

86.4 

2.2 

0.2 

2.24 

2.1 10“ 2 

1.514 

6.24 10“ 3 

0.959 

0.620 

1.320 10“ 8 

1.603 

Antarctic Plateau aerosol 

Coarse 

86 

6 

8 

- 

2.03 

1.75 

1.402 

1.00 10“ 3 

0.900 

0.831 

3.204 

-0.056 

Antarctic austral 

Nuclei 

- 

90 

10 

- 

2.03 

7.0 10“ 2 

1.530 

6.20 10“ 3 

0.960 

0.671 

7.477 10“ 6 

0.857 

winter aerosol (*) 

Accum. 

95 

- 

5 

- 

2.03 

5.0 10“’ 

1.388 

4.00 10“ 4 

0.986 

0.774 

2.724 10“ 2 

-0.157 


(*) for its use at both Antarctic coastal and high-altitude sites. 


at Zeppelin, and (ii) Tomasi et al. (2012), who obtained monthly mean 
values of<u o (0.55 pm) equal to 0.93 ± 0.04 in June, 0.86 ± 0.09 in July, 
0.91 ± 0.08 in August, and 0.89 ± 0.09 in September, yielding a seasonal 
average value of 0.90 ± 0.07 at Ny-Alesund during summer 2010. The 
values of tt>(0.55 pm) assumed in Table 5 for the Asian dust are close 
to 0.96 for the fine particle mode and 0.67 for the coarse particle 
mode, indicating that this parameter is presumably subject to decrease 
gradually as a assumes lower values as a result of the increase in the rel- 
ative super-micron dust content. The nuclei and accumulation particle 
components of the BFF smoke were both found to yield values of 
&>(0.55 pm) close to 0.90, according to Mielonen et al. (2013). 

The angular distribution curves of phase function P(0) determined 
for the eight unimodal size-distribution curves defined in Table 5 to 
represent the four Arctic aerosol types are shown on the left-hand 
side of Fig. 26. The calculations clearly show that the most pronounced 


forward scattering lobe is produced by the Asian dust coarse mode and, 
to a lesser extent, by the coarse mode of Arctic summer background 
aerosol, which yields the lowest lateral scattering (at 0 = 90°) and 
the strongest backward scattering (at 0 = 180°). 

The linear combinations of nuclei and accumulation or coarse parti- 
cle modes considered in Table 5 for the four Arctic aerosol types were 
made to vary until best-fit values of a were obtained from field sun- 
photometer measurements, which provide the unimodal column parti- 
cle number concentrations for each field-measured value of r(0.50 pm). 
The following four cases were examined: (a) the winter-spring aerosol 
case recorded at Eureka during the Arctic haze periods of 2007 and 
2008, with average values of r(0.50 pm) = 0.12 and a = 1.48, (b) the 
summer background aerosol case determined at Tiksi, giving mean 
values of r(0.50 pm) = 0.08 and a = 1.60, (c) the Asian dust transport 
episode observed at Barrow on 16 April, 2002, yielding average daily 
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Fig. 26. Angular distribution curves of phase function P(0) as a function of the scattering angle 0 for the 8 unimodal Arctic aerosol extinction models (left) and the 6 unimodal Antarctic 
aerosol extinction models defined in Table 5. 
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Fig. 27. Left-hand side: examples of bimodal particle size-distribution curves obtained as best-fit linear combinations of aerosol unimodal models for nuclei and accumulation or coarse 
particles defined in Table 3 in the following four case studies: (a) the average winter-spring aerosol case determined at Eureka (Nunavut, Northern Canada) for the mean values 
7(0.50 pm) = 0.12 and a = 1.48; (b) the summer background aerosol case determined atTiksi (Russia) in North-central Siberia for the mean values t( 0.50 pm) = 0.08 and a = 1.60; 
(c) the Asian dust episode observed at Barrow on 16 April, 2002, giving the mean daily values 7(0.50 pm) = 0.22 and a = 0.26; and (d) the BFF smoke episode observed by Stock 
etal. (2012) at Ny-Alesund on 23 March, 2008, for the daily mean values 7(0.50 pm) = 0.22 and a = 1.50. Right-hand side: as on the left, for the following four case studies: (e) the austral 
summer coastal aerosol case determined at Mario Zucchelli (MZS) for the mean values 7(0.50 pm) = 0.03 and a = 0.90; (f) the austral summer coastal aerosol case determined at 
Neumayer (NEU) for the mean values of 7(0.50 pm) = 0.045 and a = 0.78; (g) the austral summer Antarctic Plateau aerosol case determined at South Pole for the mean values 
7(0.50 pm) = 0.018 and a = 1.49; and (h) the austral winter aerosol case assumed at Neumayer for the seasonal average values 7(0.50 pm) = 0.035 and a = 0.65. 


values of r(0.50 pm) = 0.22 and a = 0.26, and (d) the BFF smoke case 
observed by Stock et al. (2012) at Ny-Alesund on 23 March, 2008, giving 
average daily values of r(0.50 pm) = 0.22 and a = 1.50. The best-fit 
bimodal size-distribution curves determined for the four cases are 
presented on the left-hand side of Fig. 27, together with the correspond- 
ing column particle number concentration of the nuclei/fine mode (N/) 
and accumulation or coarse mode ( N a/C ). The bimodal size-distribution 
curves obtained for the winter-spring case (a) and the BFF smoke case 
(d) exhibit in practice a nearly unimodal shape, since N a/C is lower 
than JV/by about 7 orders of magnitude in case (a) and by more than 
4 orders of magnitude in case (d). Conversely, the bimodal size- 
distribution curves obtained in cases (b) and (c), each using one of the 
coarse particle modes defined in Table 5, exhibit rapidly decreasing 
right-hand wings of the columnar number content over the whole 
super-micron radius range, associated with sea-salt and dust particle 
loads. 

6.2. Antarctic aerosol particle size-distributions and optical properties 

Multimodal characteristics of the aerosol size-distribution curves 
were also reported at Antarctic coastal and high-altitude sites. Analysing 
a set of size-segregated particle samples collected at Mario Zucchelli 
(Terra Nova Bay) in austral summer 1995, over the 0.035-16 pm aero- 
dynamic diameter range, Hillamo et al. (1998) found that the particle 
size-distribution usually consists of four modes: (i) a fine particle 
mode, with peak diameter D c = 0.07 pm, which contributes -1% to 
the total particulate mass content, (ii) an accumulation particle mode 
with D c « 0.30 pm, contributing 14%, (iii) a coarse particle mode with 
D c ~ 2 pm (22%), and (iv) a coarse particle mode with D c w 6.5 pm 
(63%). Similar results were obtained by Weller and Lampert (2008) 
at Neumayer, revealing a seasonal cycle characterised by nss sul- 
phate particles in austral summer, which dominated optical effects, 
and prevailing extinction by sea-salt aerosol in the other seasons. 
They found an average aerosol mass concentration at ground-level 


close to 1.1 pg m -3 in austral winter and 1.3 pg m -3 in austral summer. 
The average chemical composition of particulate mass consisted of: 
(i) 48% sea-salt, 33% nss sulphate, 12% MSA, and ~7% nitrates mixed 
with mineral dust and ammonium during the austral summer, and (ii) 
93% sea-salt particles, with a few percent of nss sulphate, nitrate, MSA, 
ammonium and mineral dust in the austral winter. Minikin et al. 
(1998) found that the nss sulphate and MSA concentrations measured 
in the austral summer at the coastal stations of Neumayer, Halley and 
Dumont d'Urville are closely correlated, showing a regular sequence 
of pronounced peaks of both such biogenic sulphur aerosol compo- 
nents, mainly formed from dimethyl sulphide (DMS). 

The aerosol chemical composition was investigated at the Finnish 
station of Aboa, about 1 50 km from the Atlantic Ocean coast in Queen 
Maud Land, by Kerminen et al. (2000); Teinila et al. (2000) and 
Koponen et al. (2003), who analysed sets of regular aerosol sampling 
measurements conducted from December 1997 to February 1998. The 
sea-salt particle concentration was estimated to be considerably lower 
than those measured at Mario Zucchelli and Neumayer. Multimodal 
features of the particle size-distribution curve were detected in most 
cases over the 0.045 < D < 15 pm range, with five principal modes: 

(i) a first mode of nuclei, over the 0.03-0.10 pm range, consisting of 
63% nss sulphate, 29% ammonium, and 8% MSA mass percentages, 

(ii) two accumulation particle modes, with average values of D c close 
to 0.30 and 0.65 pm, both consisting on average of 61% nss sulphate, 
22% MSA, and -14% ammonium, besides a few percents of sea-salt par- 
ticles, (iii) a coarse particle mode, with D c varying from 1.4 to 1.9 pm, 
and consisting of 52% sea-salt, 27% nss sulphate, 12% MSA, and minor 
percentages of nitrates and ammonium ions, and (iv) a coarse particle 
mode, with daily mean values of D c varying from 2 to 5 pm, and contain- 
ing on average 63% sea-salt, 16% nitrate, 14% nss sulphate, 5% MSA and 
2% ammonium. These compositional characteristics indicate that bio- 
genic sulphur compounds constitute on average more than 90% of the 
aerosol mass concentration of sub-micron particles, -50% of the accu- 
mulation mode particle mass, and -30% of the coarse mode particle 
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mass. Correspondingly, sea-salt particles were found to mainly contrib- 
ute to the super-micron aerosol mode by more than 60% on average. The 
above results were substantively confirmed by Virkkula et al. (2006), 
who analysed a large set of particulate samples collected in January 
2000 using a 12-stage low-pressure impactor, and observed a pair of 
sea-salt particle modes centred at D c « 0.8 pm and D c « 3 pm, respec- 
tively, and a pair of nitrate particle modes centred at diameters D c of 
-1.2 pm and -3 pm. 

The Antarctic Plateau aerosol size-distribution was found to consist 
of a fine particle mode mainly composed of nss sulphate and MSA sub- 
stances at South Pole during the austral summer (Hara et al., 2004), 
when aerosols are predominantly due to strong subsidence effects 
from the free troposphere. However, sea-salt accumulation and coarse 
particles can be often transported over the Antarctic Plateau in air 
masses with relatively high mass concentrations on days characterised 
by intense advection of oceanic air masses over the interior of Antarctic 
continent, associated with large storm systems (Shaw, 1988). Therefore, 
in the particular cases yielding values of a appreciably lower than 1.8, 
the overall particle size-distribution was represented using a bimodal 
size-distribution model consisting of a fine particle mode mainly con- 
taining nss sulphates and a sea-salt coarse particle mode. 

On the basis of the above measurements, the following size- 
distribution curves were considered: 

(a) The average size-distribution of Antarctic austral summer coastal 
aerosol was assumed to be bimodal, since the extinction effects 
produced by the two Aitken nuclei modes reported by Hillamo 
et al. (1998) were neglected. It was represented by the linear 
combination of: (i) a first particle mode consisting of fine nss sul- 
phate and sea-salt accumulation particles, with low contents of 
mineral nuclei and soot particles, and (ii) a coarse mode mainly 
composed of sea-salt particles. The composition of both modes 
was defined using the 6S (Vermote et al., 1997b) mass percent- 
ages reported in Table 5, according to the composition data calcu- 
lated by Tomasi et al. (2012) at Mario Zucchelli and Neumayer. In 
both particle modes, the soot mass concentration was assumed 
taking into account the low values of BC mass concentration 
measured by Wolff and Cachier (1998) at Halley during the aus- 
tral summer, ranging in general from 1.0 to 2.0 ng m -3 . 

(b) The average size-distribution of austral summer background 
aerosol over the Antarctic Plateau was represented by: (i) a 
fine particle mode mainly containing nss sulphates and minor 
percentages of sea-salt and dust particles, with a soot content 
equal to 40% of that estimated by Wolff and Cachier (1998) at a 
coastal site, and (ii) a coarse particle mode composed mainly 
of sea-salt particles, and containing lower mass percentages of 
water-soluble and dust-like 6S (Vermote et al., 1997b) compo- 
nents, as given in Table 5. 

(c) The average size-distribution of Antarctic austral winter aerosol 
at coastal sites was not based on sun-photometer measurements, 
as such measurements have neither been conducted during the 
austral winter at coastal sites (for 66°-75° S latitudes) nor at 
Antarctic Plateau stations. In order to achieve reliable evaluations 
of complex refractive index n(\)-i k( A) and the other radiative 
parameters during the austral winter, the chemical composition 
data determined by Minikin et al. (1998) at Neumayer and 
other Antarctic coastal sites over the 14-year period from 1983 
to 1996 were taken into account. These in-situ measurements 
showed a regular sequence of pronounced mass concentration 
minima of biogenic sulphur aerosol components originating 
from DMS in austral winter and some marked peaks in summer. 
In particular, they found that: (i) the nss sulphate mass concen- 
tration decreased on average from 50 to 4 ng m“ 3 (i.e. by 92%) 
from austral summer to winter, and (ii) MSA mass concentration 
decreased from 17 to 2 ng m -3 (i.e. by more than 88%), because 
of a strong reduction in biogenic sources of both chemical species 


during the local winter. Conversely, the sea-salt particle mass 
concentration was estimated by Weller and Lampert (2008) to 
increase at Neumayer by more than 40%, passing from an average 
summer value of 597 ± 830 ng m“ 3 to an average winter value 
of 844 ± 1 100 ng m“ 3 . Despite the plausible strong decrease in 
sea-salt particle concentration caused by the ice coverage of the 
ocean areas near the Antarctic continent, the observed austral 
winter increase in the sea-salt concentration was probably due 
to the more intense transport of maritime particles from the 
lower latitude oceanic ice-free areas characterised by higher 
wind speed conditions at the sea-surface. 

On the basis of these remarks, it was decided to use a bimodal model 
to represent the austral winter aerosol at coastal sites, consisting of: (i) a 
water-soluble nuclei mode, and (ii) an accumulation particle mode 
composed mainly of sea-salt particles transported from remote ocean 
areas. Each mode of the three Antarctic bimodal size-distribution curves 
was assumed to have an analytical form given by Eq. (3) for the values of 
shape-parameters r c and a reported in Table 5, and to consist of the mass 
percentages defined in Table 5 for the four 6S (Vermote et al., 1997b) 
basic components. 

Because of the predominant mass fractions of nss sulphate in the nu- 
clei/fine particle mode and of sea-salt in the accumulation or coarse 
mode, the three bimodal size-distributions of Antarctic particles consid- 
ered above provided the unimodal average values of refractive index 
parts n(0.55 pm) and k( 0.55 pm) given in Table 5. For these aerosol ex- 
tinction models, values of <n(0.55 pm) ranging from -0.80 (for the 
coarse mode of coastal aerosols) to 0.99 (for the accumulation mode 
of austral winter aerosol) were obtained (see Table 5). These estimates 
of <u(0.55 pm) agree very well with the in-situ evaluations made by: 
(i) Virkkula et al. (2006a, b) at Aboa, who obtained values of 
<n o (0.55 pm) > 0.95 for about 93% of in-situ measurements, and <0.90 
for only 3.5% of cases, (ii) Weller and Lampert (2008) at Neumayer, 
who estimated that co o (0.55 pm) varied from 0.97 to 1.00 in more 
than 95% of samples, and (iii) Tomasi et al. (2012), who found values 
of <i> o (0.55 pm) mainly ranging from 0.95 to 0.98 for the particulate 
chemical composition measurements made by Hillamo et al. (1998) 
and Fattori et al. (2005) at Mario Zucchelli, and those conducted by 
Weller and Wagenbach (2007); Weller et al. (2008) and Weller and 
Lampert (2008) at Neumayer. These findings are in good agreement 
with the in-situ evaluations made by: (i) Bodhaine (1995), who obtain- 
ed monthly mean values of w o (0.55 pm) at South Pole, varying from 
0.942 to 0.972 during the austral summer, and an annual mean value 
of <» o (0.55 pm) = 0.97, (ii) Heintzenberg et al. (1997), who determined 
an annual average value of co 0 = 0.965 in the visible at South Pole; and 
(iii) Tuncel et al. (1989) and Arimoto et al. (2004) at South Pole, and 
Piel et al. (2006) at Kohnen, who obtained values of <u o (0.55 pm) ranging 
from 0.95 to 0.98. 

The daily mean values of t( 0.55 pm) retrieved by Campbell et al. 
(2012) from the 2007 MODIS and MISR satellite observations recorded 
over Antarctica were all lower than 0.10 throughout the whole year. 
Reliable measurements of AOT were not retrieved by Winker et al. 
(2010) from the CALIOP lidar observations during the Antarctic “polar 
night”, since they were lower than the lidar detection limit. This was 
confirmed by Winker et al. (2013), who showed that the CALIPSO 
mean cloud-free day-time and night-time values of r(0.532 pm) were 
all lower than 0.05 over the Antarctic continent and the peri-Antarctic 
regions during the long austral winter, from March to November. 
Therefore, the AOT contribution exhibited by nuclei is in general 
<0.01 in the austral winter months, while that of sea-salt accumulation 
particles increases appreciably yielding daily values of t( 0.50 pm) vary- 
ing from 0.02 to 0.05 during the coldest season, and values of a mainly 
ranging from 0.50 to 0.70 over the peri-Antarctic regions (Winker et al., 
2013). 

The angular diagrams of phase function P(0) obtained for the 
unimodal size-distribution curves defined in Table 5 for the three 
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Antarctic aerosol types are presented on the right-hand side of Fig. 26 as 
a function of scattering angle 0. The most pronounced forward scatter- 
ing lobe is given by the coarse particle mode of the austral summer 
coastal aerosol, while the austral summer coarse particle mode of the 
Antarctic Plateau size-distribution yields a slightly lower forward scat- 
tering lobe and the most marked backward and lateral (at 0 = 90°) 
scattering intensities. As a result, four case studies of Antarctic aerosol 
were analysed to determine the corresponding best-fit linear combina- 
tions of a fine particle mode and an accumulation or coarse particle 
mode, for the following sets of t( 0.50 pm) and a: (i) the austral summer 
coastal aerosol case (e) obtained at Mario Zucchelli for the mean values 
of 7(0.50 pm) = 0.030 and a = 0.90, (ii) the austral summer coastal 
aerosol case (f) determined at Neumayer for the mean values of 
r(0.50 pm) = 0.043 and a = 0.68, (iii) the austral summer Antarctic 
Plateau aerosol case (g) defined at South Pole for the mean values of 
r(0.50 pm) = 0.018 and a = 1.49, and (iv) the austral winter aerosol 
case (h) defined on the basis of in-situ aerosol measurements conduct- 
ed during the austral winter months at Neumayer and evaluated to give 
the average seasonal values of r(0.50 pm) = 0.035 ± 0.015 and a = 
0.65 ± 0.10. The bimodal size-distribution curves obtained by applying 
the best-fit procedure for the four values of a given above are shown 
on the right-hand side of Fig. 27, where the columnar number contents 
Nf and N a/C are also given. The size-distribution curves obtained in cases 
(e), (f) and (g) by assuming the presence of a coarse mode are 
characterised by robust right-hand wings over the super-micron radius 
range. In fact, in the two coastal cases (e) and (f), the columnar number 
content N a/C is lower than Nf by more than 6 and 5 orders of magnitude, 
respectively. Fig. 27 also shows that the relative contribution of sea-salt 
coarse mode particles is appreciably higher at Neumayer than at Mario 
Zucchelli, according to the results shown in Figs. 9 and 11, respectively. 
Considerably different values of Nf and N a/C by at least nine orders 
of magnitude were found in case (g) for the Antarctic Plateau aerosol, 
giving a measure of the large differences characterising the atmospheric 
concentrations of nss sulphate nuclei mode and sea-salt coarse mode. 
By contrast, the accumulation mode used in case (h) to represent the 
sea-salt particles during the austral winter seems to contribute more 
weakly (by only 4 orders of magnitude) to enhance N a/C , with respect 
to Nf. 

6.3. Evaluations of direct aerosol-induced radiative forcing effects in polar 
regions 

Atmospheric aerosols are known to affect the radiation balance 
of the surface-atmosphere system (i) directly through interaction 
with solar (short-wave) radiation (Haywood and Boucher, 2000) and 
terrestrial (long-wave) radiation (Lubin et al., 2002), and (ii) indirectly, 
by acting as cloud condensation nuclei and modifying cloud albedo 
characteristics (Schwartz and Andreae, 1996). They can considerably 
modify the horizontal and vertical distribution of radiant energy passing 
through the atmosphere (Stone et al., 2008), generally causing more 
intense changes in the outgoing flux of solar radiation at TOA-level 
than those affecting the long-wave radiation, as clearly indicated by 
Mie (1908) theory. For this reason, the direct aerosol-induced radiative 
forcing (hereinafter referred to as DARF) is commonly evaluated by only 
considering variations induced by airborne aerosols on the short-wave 
radiation flux and by neglecting those affecting long-wave radiation. 
Instantaneous DARF effects can vary on a diurnal basis as a function of 
solar zenith angle 0 O . and closely depend on: (a) aerosol scattering 
and absorption characteristics (represented in terms of parameters 
r(\), a, n(A), k(\), P(Q). g(A), and co(\) at solar wavelengths), and 
(b) the spectral and geometrical features of surface reflectance. Three 
different instantaneous DARF terms are usually considered in these 
calculations: (i) at the TOA-level, as the difference AF T0A (t) estimated 
at a certain time t between the upward solar radiation fluxes emerging 
from the real atmosphere with aerosols and from the pristine atmo- 
sphere without aerosols (Hanel et al„ 1999), (ii) at the bottom-of- 


atmosphere (BOA) level (i.e. at the surface), as the difference AF B oA(t ) 
estimated at a certain time t between the net short-wave fluxes 
determined at surface-level in the atmosphere with aerosols and 
in the same atmosphere assumed without aerosols (Satheesh and 
Ramanathan, 2000) and (iii) within the atmosphere, as the difference 
A F AT M(t) between terms A F TOA (t) and AFboaM (Ramanathan et al., 
2001 ). To give an average measure of the daily DARF effects, the instan- 
taneous terms AF TO a(f). AF B0A (t) and AF ATM (t) are generally evaluated 
at pre-fixed hours of the day, and then integrated over the entire period 
from sunrise to sunset, and divided by the 24-hour period. The corre- 
sponding diurnally averaged DARF terms A F T0A , A F B oa and AF A tm were 
calculated by following such a procedure, based on the Bush and 
Valero (2003) criteria. Negative values of A F TO a indicate that aerosols 
cause an increase in the radiation budget of the surface-atmosphere 
system at TOA-level, producing direct cooling effects on the climate 
system, while positive values of A F T0A indicate the occurrence of a 
warming effect during the day. Term A F B0A at the surface gives a mea- 
sure of the perturbation induced by airborne aerosols in the net flux 
reaching the surface, which may be positive (warming) or negative 
(cooling), while term AF atm = A F T oa — AF boa defines the radiation sur- 
plus (or deficit) induced by aerosols within the atmosphere. In practice, 
AF atm is the amount of latent heat released in the atmosphere by 
aerosols and internally redistributed with the effects of modifying 
both temperature gradients and atmospheric circulation. 

Calculations of A F TOA , A F boa and A F ATM have been made by Tomasi 
et al. (2014) using the radiative transfer 6S code of Vermote et al. 
(1997a) for the aerosol loads and composition, and optical characteris- 
tics defined in Table 6, pertaining to 1 1 Arctic aerosol types and 4 
Antarctic background austral summer aerosol types. The calculation 
procedure consisted of the following seven steps: 

(i) Determination of the columnar aerosol extinction parameters 
r(A) and a from the field sun-photometer measurements col- 
lected at Barrow, Ny-Alesund, Summit, Sodankyla, Tiksi, Mario 
Zucchelli, Neumayer, Dome Concordia and South Pole, from 
which the average values of r(0.50 pm) and a( 0.40-0.87 pm) 
given in Table 6 were obtained. 

(ii) Calculation of the spectral values of parameters n( A) and k(\) for 
the chemical composition data based on the OPAC mass percent- 
ages of fine and coarse particles defined by Hess et al. (1998) for 
the six OPAC basic components reported in Table 6 at 50% rela- 
tive humidity of air. 

(iii) Definition of the multimodal linear combinations of the size- 
distribution curves for the 15 aerosol types considered in 
Table 6 to fit the values of a(0.40-0.87 pm) given in Table 6. 

(iv) Calculation of the aerosol single scattering albedo co( 0.55 pm) 
reported in Table 6 for the 15 multimodal size-distribution 
curves defined at the previous step. 

(v) Selection of the BRDF surface reflectance models of Tomasi et al. 
(2013) to represent the surface reflectance characteristics of the 
land surfaces described in terms of MODIS Level 3.0 MCD43C3 
products around the principal sun-photometer stations, and 
those of Arctic and Antarctic oceanic surfaces. The analysis of 
such observational data suggested that the following models 
were suitable in the present calculations: (a) the BRDF ocean sur- 
face reflectance model OS1, defined for surface wind velocity 
V w = 2 m s^ 1 and giving a value of broadband albedo A b (9 0 = 
60°) = 0.193, (b) the BRDF surface reflectance model VS1, yield- 
ing a value ofA b (0o = 60°) = 0.155, and used to represent a land 
area covered by tundra and conifer forests in the surroundings of 
Barrow, Sodankyla and Tiksi, (c) the BRDF land surface reflec- 
tance model PS1 for land surfaces covered by snow fields and 
glaciers, giving a value of A b (6 0 = 60°) = 0.854, and used to 
represent the surface reflectance characteristics observed around 
Summit in the Arctic and Neumayer, Dome Concordia and South 
Pole in Antarctica for fresh snow conditions, (d) the BRDF surface 


Table 6 

Multimodal aerosol extinction models based on the OPAC (Hess et al., 1 998) components used to represent 4 Arctic and 4 Antarctic columnar contents of polar aerosol, with (i) the mass percentages of the OPAC components of fine and accumulation/ 
coarse particle modes in the atmospheric column; (ii) the values of real part n and imaginary part k of complex refractive index at the 0.55 pm wavelength, (iii) the values of aerosol optical thickness 7(0.50 pm) and Angstrom's exponent a calculated 
over the 0.40-0.87 pm wavelength range, and (iv) the values of columnar aerosol single scattering albedo co(0.55 pm). The last three columns provide the values of diurnal average aerosol radiative forcing terms at the TOA-level (A F T oa), at the surface 
(A F boa ) and within the atmosphere (A F ATM ), estimated for the BRDF surface reflectance models OS1 (oceanic), VS1 (vegetation-covered) and PS1, PS2, PS3 and PS4 (polar ice- and snow-covered) determined by Tomasi et al. (2014). The acronym BG 
stands for “background”. The OPAC aerosol components are indicated by the acronyms WASO (water-soluble), SSAM (sea-salt accumulation mode), SSCM (sea-salt coarse mode), 1NSO (insoluble), SOOT (soot), and MITR (mineral-transported). 

Aerosol type Mass percentages of the basic OPAC components calculated Particulate matter Aerosol optical Exponent co BRDF surface Diurnal average aerosol 

for relative humidity = 50% to give form to the fine and refractive index parts thickness a: (0.55 pm) reflectance radiative forcing (W m -2 ) 

accumulation/coarse particle modes r(0.50 pm) model 


Fine particle modes Accumulation/coarse particle modes 



WASO 

SOOT 

SSAM 

SSCM 

INSO 

MITR 

n(0.50 pm) 

k(0.50 pm) 





AF^oa 

AFbo/i 

A Fatm 

Winter-spring Arctic haze (Barrow) 

95.8 

4.2 

22.9 

37.7 

39.4 

- 

1.399 

3.3 10“ 3 

0.116 

1.28 

0.840 

OS1 

-5.5 

+6.5 

-12.0 













PS1 

+ 9.0 

+ 0.0 

+ 8.9 













PS2 

+ 8.8 

+ 0.1 

+8.7 

Winter-spring Arctic haze (Ny-Alesund) 

98.0 

2.0 

- 

90.9 

9.1 

- 

1.424 

7.3 1(T 3 

0.080 

1.29 

0.949 

OS1 

-5.3 

-0.6 

-4.7 













PS3 

+1.0 

-0.4 

+1.4 

Winter-spring Arctic haze (Sodankyla) 

95.3 

4.7 

94.8 

5.2 

- 

- 

1.399 

3.9 1(T 3 

0.066 

1.25 

0.840 

OS1 

-7.7 

-1.1 

-6.6 













VS1 

-3.5 

-5.3 

+1.8 













PS2 

+ 2.2 

-0.3 

+2.5 

Background summer aerosol (Barrow) 

96.8 

3.2 

11.0 

43.6 

45.4 

- 

1.444 

9.2 10“ 3 

0.078 

1.40 

0.978 

OS1 

-9.2 

-1.9 

-7.3 













VS1 

-4.2 

-6.5 

+ 2.3 

Background summer aerosol (Ny-Alesund) 

98.0 

2.0 

- 

50.0 

50.0 

- 

1,461 

7.9 10“ 3 

0.041 

1.20 

0.966 

OS1 

-7.2 

+1.5 

-8.6 













PS4 

+1.0 

-2.3 

+3.4 

Arctic dense summer aerosol (Ny-Alesund) 

97.6 

2.4 

6.3 

59.7 

34.0 

- 

1.437 

3.7 10“ 3 

0,120 

1.00 

0.852 

OS1 

-10.4 

+4.2 

-14.6 













PS4 

+ 13.0 

-5.2 

+ 18.2 

Background summer aerosol (Summit) 

98.0 

2.0 

- 

50.0 

50.0 

- 

1.449 

8.8 10“ 3 

0.039 

1.48 

0.969 

PS1 

+ 2.2 

-0.2 

+2.5 

Background summer aerosol (Sodankyla) 

98.0 

2.0 

- 

50.0 

50.0 

- 

1.452 

8.6 10“ 3 

0.060 

1.42 

0.965 

OS1 

-6.5 

-1.6 

-4.9 













VS1 

-3.0 

-5.5 

+2.4 













PS2 

+ 3.0 

-0.4 

+3.4 

Background summer aerosol (Tiksi) 

98.0 

2.0 

- 

50.0 

50.0 

- 

1.444 

9.2 10“ 3 

0.085 

1.60 

0.977 

OS1 

-10.3 

-1.6 

-8.7 













VS1 

-5.6 

-9.0 

+ 3.3 

Asian dust (Barrow) 

100.0 

- 

5.8 

2.4 

56.5 

35.3 

1.527 

6.3 10“ 3 

0.200 

0.80 

0.858 

OS1 

-9.3 

+ 16.0 

-25.3 













PS2 

+ 18.5 

+ 0.3 

+ 18.2 

Boreal forest fire smoke (Barrow) 

95.5 

4.5 

21.7 

26.1 

52.2 

- 

1.469 

2.5 10“ 3 

0.300 

1.20 

0.758 

OS1 

-32.9 

-6.7 

-26.2 













VS1 

-20.2 

-28.6 

+ 8.3 

Background austral summer aerosol 

99.5 

0.5 

- 

73.9 

26.1 

- 

1.468 

2.3 10“ 3 

0.030 

0.90 

0.964 

OS1 

-6.2 

+ 0.9 

-7.1 

(Mario Zucchelli) 












PS2 

+ 2.0 

-0.1 

+2.1 

Background austral summer aerosol 

99.7 

0.3 

- 

85.6 

14.4 

- 

1.457 

1.3 10“ 3 

0.043 

0.68 

0.975 

OS1 

-8.2 

+1.4 

-9.6 

(Neumayer) 












PS1 

+1.9 

-0.2 

+2.1 













PS2 

+2.1 

+ 0.0 

+2.1 

Background austral summer aerosol 

100.0 

- 

4.5 

- 

- 

95.5 

1.441 

2.0 10“ 3 

0.019 

1.77 

0.999 

PS1 

+ 0.4 

+ 0.2 

+ 0.2 

(Dome Concordia) 












PS2 

+ 0.6 

+ 0.4 

+ 0.2 

Background austral summer aerosol 

99.6 

0.4 

- 

83.6 

16.4 

- 

1.445 

3.2 10“ 3 

0.018 

1.49 

0.988 

PS1 

+ 0.6 

-0.1 

+0.7 

(South Pole) 












PS2 

+0.7 

-0.1 

+ 0.8 
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albedo model PS2, yielding a value of A b (6 0 = 60°) = 0.775, 
which was adopted to represent the surface reflectance charac- 
teristics for land areas covered by slightly polluted snow (at 
Sodankyla) and not entirely snow-covered areas (Mario Zucchelli, 
Neumayer, and other coastal sites), (e) the BRDF model PS3, 
giving A b (8 0 = 60°) =0.564, which was chosen to represent the 
surface albedo characteristics observed over Spitsbergen in late 
winter and early spring, and (f) the BRDF model PS4, yielding 
A b (0 o = 60°) = 0.329, which realistically represents the mixed 
ice-covered and ice-free land areas around Ny-Alesund in late 
spring and summer. 

( vi) Calculations of the daily time-patterns of instantaneous and diur- 
nally averaged DARFs were made using the 6S (Vermote et ah, 
1997a) radiative transfer code for: (a) the values of r(0.50 pm) 
and ct(0.40-0.87 pm) defined in Table 6 for the various aerosol 
types and assumed to be stable from sunrise to sunset, (b) the 
values of angle 6 0 calculated at the various hours of the day, for 
the above stations and seasonal periods, (c) the multimodal 
size-distribution curves determined at step (iii), (d) the spectral 
values of <u(A) for each of the multimodal size-distribution curves 
determined at step (iv), and (e) the surface albedo models chosen 
at step (v) for the various Arctic and Antarctic sites. The time- 
patterns of AF TOA (t), AF B0A (t) and AF ATM (t ) were then integrated 
to calculate the daily values of A F TOA , AF boa and A F ATM given in 
Table 6 for the various aerosol types and BRDF surface albedo 
models. 

Examining the results given in Table 6, it can be seen that AF toa 
varies appreciably at Barrow as one passes from the OS1 to VS1 model 
in the Asian dust (AD) case, and from the OS1 to PS2 model in the 
BFF smoke case, while only limited variations of AF toa from about 
— 10W m -2 to no more than + 10 W m -2 were obtained for the back- 
ground summer aerosol and Arctic haze cases considered at the 5 
above-chosen Arctic sites. Similar features were also noted for A F B0A , 
illustrating larger variations in the AD and BFF aerosol cases monitored 
at Barrow, and mainly ranging from — 5 to + 5 W m“ 2 in all the back- 
ground summer aerosol and Arctic haze cases. Correspondingly, rather 
large variations in A F ATM were obtained not only for the Barrow AD 
and BFF cases but also for the Arctic haze cases observed at Barrow 
and Ny-Alesund for t(0.50 pm) > 0.10. In particular, we have estimated 
an Asian dust value of A F T0A « 9 W m~ 2 for r(0.50 pm) = 0.20 mea- 
sured at Barrow. 

The diurnally averaged DARF terms A F T0A , A F B0A and A F ATM given in 
Table 6 for four Antarctic background austral summer aerosol cases in- 
dicate that larger variations have been estimated at Mario Zucchelli 
and Neumayer, passing from the low surface reflectance oceanic 
model OS1 to the land surface albedo model PS2 characterised by highly 
reflecting and snow-covered surfaces. Conversely, very limited varia- 
tions in the three DARF terms have been obtained at the four Antarctic 
sites, when passing from one PS land surface albedo model to another, 
since they all pertain to land surfaces covered by glaciers and snow 
fields that do not exhibit marked differences in their reflectance charac- 
teristics. Actually, it is important to mention that the above DARF 
evaluations may be affected by large uncertainties, associated with the 
numerous approximations made in the above calculations when: 

(a) assuming constant time-patterns of AOT and stable aerosol optical 
parameters throughout the day, (b) neglecting the occurrence of varia- 
tions in the scattering, absorption and extinction properties of aerosol 
with height, and assuming that the vertical profiles of such optical 
parameters did not exhibit multi-layered features, (c) using simplified 
radiative transfer codes in a cloud-free atmosphere, which can lead to 
significant errors in estimating the DARF effects, as pointed out by 
Bush and Valero (2003), and (d) choosing inappropriate surface reflec- 
tance models to represent the surface albedo characteristics, as would 
especially be the case for the VS1 and OS1 models, which may differ 


considerably in some cases from the real surface conditions, as pointed 
out by Stone et al. (2008). With regard to this, calculations of the DARF 
effects made by Tomasi et al. (2014) using the sets of four oceanic sur- 
face (OS) and four vegetation-covered (VS) BRDF models indicate that 
changes varying from 1 0% to 50% may affect the calculations of A F TOA 
and AF BO a, and consequently even greater differences can be found in 
evaluating AF atm . Therefore, it is fair to highlight that the present eval- 
uations of AFtoa, A F BOA and A F ATM obtained using the OS1 calm-wind 
ocean model could differ appreciably from those occurring for stronger 
winds causing higher surface reflectance conditions in sea areas. 

Calculations of the DARF efficiency parameters E T0A , E BOA and E ATM 
were finally made by dividing the values of A F TOA , A F BOA and A F ATM 
given in Table 6 for the various polar aerosol types by the corresponding 
mean values of t(0.50 pm). The evaluations of E TOA , E BOA and E ATM are 
shown in Fig. 28 as a function of the broadband albedo A b (8 0 = 60°) 
determined for the various BRDF surface reflectance models chosen 
above. The values of E ToA , E BoA and E Atm obtained at the five Arctic sites 
for various aerosol types indicate that DARF efficiencies are subject to 
large variations when passing from the low surface albedo characteris- 
tics of the VS1 and OS1 models to the high reflectance surfaces of the PS 
models. The results shown in Fig. 28 indicate that: 

(a) £ TO a increases from less than — 100 W m“ 2 over the VS1 surface 
(for BG summer aerosol over sea, at Ny-Alesund) to more than 
+ 100 W m -2 over the PS4 surface (for Arctic dense summer 
aerosol at Ny-Alesund), and assumes slightly positive values as 
the broadband albedo increases from 0.56 to 0.85, ranging from 
about +30 W m“ 2 (over the PS2 surface for Arctic haze at 
Sodankyla and for BG summer aerosol at Dome Concordia) to 
no more than +92 W m~ 2 (over the PS2 surface, for Asian 
dust at Barrow). These estimates had relative standard deviations 
varying on average from 10% over the PS4 surface to 50% over 
the VS1 and OS1 surfaces (i.e. from 6 to 26 W m~ 2 ). 

(b) E BOA varies from less than — 1 00 W m~ 2 over the VS1 surface (for 
BG summer aerosol at Tiksi) to around + 80 W m“ 2 over the OS1 
surface (for Arctic haze at Barrow). The estimates of E BOA obtain- 
ed for the high-reflectance PS surface reflectance models did not 
vary by much with the aerosol type, ranging from — 57 W m~ 2 
(for Arctic BG aerosol at Ny-Alesund) to + 20 W irrf 2 (for Antarctic 
BG summer aerosol at Dome Concordia). 

(c) Fa™ varies from about — 240 W m 2 over the OS1 surface (for 
Antarctic BG summer aerosol at Mario Zucchelli) to more than 
+ 150 W m -2 over the PS4 surface (for Arctic dense summer 
aerosol at Ny-Alesund). Fig. 28 shows that E ATM exhibits a 
large range over the VS1 and OS1 surfaces, and then assumes 
stable and positive values over the PS surfaces, varying from 
+ 10 W m~ 2 (for Antarctic BG summer aerosol at Dome 
Concordia) to + 151 W m -2 over the PS4 surface (for dense 
Arctic haze at Ny-Alesund). 

The present DARF efficiency evaluations are affected by considerable 
uncertainties, since they were determined with relative standard devia- 
tions vaiying from about ±10% (over the PS4 surface) to around ± 50% 
(over the VS1 and OS1 surfaces). 

The results shown in Fig. 28 have been obtained for various Arctic 
and Antarctic aerosol types with values of to(0.55 pm) ranging from 
about 0.76 (for BFF smoke at Barrow) to nearly 1.00 (for Antarctic back- 
ground aerosol at Dome Concordia). The value of E TOA for the BFF smoke 
load observed at Barrow was found to be of — 68 ± 27 W m“ 2 over the 
tundra (VS1) surface, while Stone et al. (2008) estimated an average 
efficiency value of — 20 W m“ 2 in summer 2004. The corresponding 
value of £ B oa was estimated by us at — 95 ± 36 W m -2 , while Stone 
et al. (2008) estimated a lower average efficiency of about — 40 W m“ 2 . 
Thus, we obtained a value of Fa™ equal to +27 ± 63 W m -2 , against 
an estimate of +20 W m“ 2 by Stone et al. (2008), resulting in an atmo- 
spheric warming rate of about 1 K per day. 
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Fig. 28. Scatter plots of the daily mean values of DARF efficiencies E T oa at the TOA-level (upper part). E B oa at the BOA-level (middle part), and E A m in the atmosphere (lower part) shown 
versus the broadband albedo calculated by Tomasi et al. (2014) at solar zenith angle 6 0 = 60°, for the BRDF oceanic surface (OS1 ), vegetation-covered surface (VS1 ) and snow-covered 
polar surface (PS1, PS2, PS3 and PS4) models, and for the 15 polar aerosol types defined in Table 6 and represented using different symbols. 


7. Conclusions 

Monthly mean values of r(0.50 |jm) and Angstrom exponent a 
determined from multi-year sets of ground-based sun-photometer 
measurements at 21 polar sites, yielded a comprehensive picture of 
atmospheric turbidity conditions due to polar aerosols. The results 
are presented as a function of the winter-spring and summer periods 
over the whole Arctic region and include observations conducted in 
Alaska, Northern Canada, Greenland, Svalbard, Northern Scandinavia 
and Northern-Central Siberia. The analysis of the measurements ac- 
quired at the 12 Arctic sites listed in Table 1 highlighted the seasonality 
of columnar aerosol extinction parameters, showing that a ranges 
mainly from 1.0 to 1.6 throughout the year at all stations, and 
r(0.50 pm) values were mainly lower than 0.08 in summer (for 
background aerosol) and appreciably higher during the winter-spring 
period. The latter phenomenon is associated with the occurrence of 
dense Arctic haze episodes that are often observed at remote Arctic 
sites as the result of the transport of anthropogenic polluted aerosols 
from the industrialised and most densely populated areas of mid- 
latitude North America, Europe and Asia. 

Sun-photometer measurements collected during the Antarctic 
summer were analysed for six coastal sites (Marambio, Neumayer, 
Novolazarevskaya, Mirny, Syowa, Mario Zucchelli), one mid-altitude 
station (Troll) and two high-altitude bases on the Antarctic Plateau 
(Dome Concordia and South Pole). Monthly mean r(0.50 pm) values 
usually ranged from less than 0.02 on the Antarctic Plateau to no more 
than 0.06 at the coastal stations, with a decreasing respectively from 
extremes of 1.8 (at South Pole and Dome Concordia, for aerosols domi- 
nated by nss sulphates), to around 0.6 at Mario Zucchelli, Neumayer and 
the other coastal sites (where aerosol light extinction is predominantly 
influenced by sea-salt particles). 

Ship-borne sun-photometer measurements were conducted during 
14 spring and summer AERONET/MAN cruises across the Greenland Sea 
and Norwegian Sea, the West Siberian Sea and the North American Arctic 
Ocean. The results yielded monthly mean estimates of r(0.50 pm) and a 
that were found to be comparable with results obtained from ground- 
based sun-photometer measurements taken at the stations located 
on the coasts of the same oceanic sectors during the March to September 


period, thus confirming the covariance features of t( 0.50 pm) and a 
determined at the ground-based sun-photometer sites. The analysis of 
ship-borne sun-photometer measurements performed during the 18 
AERONET/MAN cruises conducted in the coastal and off-shore Antarctic 
areas of the Indian, Pacific and Atlantic Oceans and along the Antarctic 
Peninsula coasts showed that the monthly mean values of t( 0.50 pm) 
varied between the extremes of 0.02 and 0.06 from December to April, 
while decreasing from 0.12 in October to 0.07 in November for measure- 
ments acquired around the Antarctic Peninsula. The large magnitude 
AOTs in the latter case were likely the result of stronger winds from the 
northern and western quadrants producing more sea-salt particles, 
which are transported towards the continent, while the rapidly decreas- 
ing trend results from a rapid decrease in the amplitude of these winds 
(a contention that is further supported by the corresponding increase in 
a values, between October and January, while r(0.50 pm) systematically 
decreases). In more general terms we would argue that scattergrams of 
a versus t( 0.50 pm) for the ensemble of both land-based and ship- 
borne based sun-photometry showed a roughly inverse relationship 
of high a, small r( 0.50 pm ) to low a, high t( 0.50 pin ) that was represen- 
tative of a progressive transformation from nuclei dominated nss 
sulphates above the Antarctic Plateau to sea-salt dominated aerosols 
for coastal sites. 

The vertical distribution of the aerosol extinction coefficient was 
investigated with ground-based lidars at Arctic and Antarctic sites. In 
particular, the monthly mean vertical profiles of backscatter coefficient, 
along with lidar ratio trends and depolarisation ratio covariance with 
backscattering coefficient were illustrated for tropospheric altitudes 
using data acquired at Ny-Alesund by the AWI KARL lidar from 1 
November, 2012, to 31 October, 2013. As part of this investigation, we 
illustrated the seasonal variations of integrated backscatter coefficient 
and lidar ratio, as influenced by Arctic haze episodes in late winter 
and spring, and the background aerosol in summer. 

Evaluations of AOT in the visible have been obtained using the 
operational algorithms applied to MODIS (Levy et al., 2013; Glantz 
et al., 2014), MISR (Diner et al., 2005b; Kahn et al., 2010) and AATSR 
(Holzer-Popp et al., 2013) satellite observations, which were found to 
substantially agree with the values of AOT obtained from the ground- 
based and ship-borne sun-photometer measurements conducted in 
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the Arctic and Antarctic open-water regions. Reliable evaluations 
of t( 0.55 pm) have also been retrieved from the AATSR measurements 
recorded over the Arctic oceanic and land regions covered by sea ice, 
land ice and snow using observations taken in different spectral chan- 
nels (Istomina et al., 2011), and by using the DVMS method (Mei 
et al., 2013b). These retrievals allow the aerosol spatial distribution 
(e.g., pollution events) to be analysed on short term scales, but may 
contain an AOT offset due to many challenges connected to virtually 
unknown surface types and aerosol optical properties. The global appli- 
cation of these procedures is at present a challenging problem because 
of the very high sensitivity of the retrievals to unscreened thin clouds. 
In fact, these data are already affected by large uncertainties and need 
to be further validated, to overcome the errors that may arise from 
cloud contamination and the complex picture of polar conditions (low 
values of AOT, various aerosol/cloud types, atmospheric inversions, 
variable features of snow- and ice-covered surface reflectance ). 

The main characteristics of multimodal particle size-distribution 
curves have been defined by: (i) taking into account the in-situ aerosol 
sampling measurements conducted at various Arctic sites (Barrow, Ny- 
Alesund, Sodankyla) and Antarctic sites (Mario Zucchelli, Neumayer, 
Aboa, South Pole, Dome Concordia), and (ii) assuming that the overall 
size-distribution curves determined at the various sites consist in 
general of a nuclei mode and an accumulation or coarse particle mode, 
linearly combined to give a bimodal aerosol columnar load for each 
value of a obtained from the field measurements. Eight log-normal 
curves were defined for the Arctic stations on the basis of in-situ particle 
samples collected for Arctic haze, summer background aerosol, Asian 
dust and BFF smoke. They can be used for each spectral set of sun- 
photometer t{\) measurements to fit the derived value of tx The radia- 
tive parameters of the eight unimodal curves were defined by taking 
into account the complex refractive index and single scattering albedo 
derived from in-situ measurements conducted with nephelometers 
and aethalometers and/or using the simulations of the aerosol optical 
properties made with the aerosol extinction models given by the 6S 
code. Six log-normal curves were defined to represent the unimodal 
size-distribution curves of Antarctic nuclei, accumulation and coarse 
particles, using field data obtained from the aerosol collected at coastal 
and high-altitude sites. The optical properties of particulate matter were 
defined according to the data provided by the in-situ active remote 
sensing measurements, and used to determine the optical characteris- 
tics of six Antarctic log-normal size-distribution curves shown in 
Table 5. 

Assumptions about the aerosol optical properties are totally based on 
the particulate chemical composition measured at the ground and, there- 
fore, may differ appreciably from airborne particles suspended within the 
upper layers of the atmosphere in all cases of aerosol transport from 
remote areas (like Arctic haze, Asian dust, and BFF smoke), as shown 
for instance by the PAMARCMiP airborne measurements conducted in 
April 2009 (Stone et al., 2010). However, the best-fit procedure combin- 
ing the nuclei mode with the accumulation or coarse particle mode to fit 
the value of a will have probably allowed us to evaluate the average 
spectral features of the column loading consisting of larger particles. 

The bimodal size-distribution curves of Arctic and Antarctic aerosol 
obtained as linear combinations of the nuclei or fine particle modes 
and accumulation or coarse particle modes have been defined by vary- 
ing their columnar number contents until a best-fit value of a. was ob- 
tained for each spectral series of r(A). All these bimodal models were 
assumed to exhibit the well-defined optical properties given in 
Table 5 and to represent the most relevant optical characteristics of Arc- 
tic and Antarctic columnar aerosol. Additional realistic multimodal ex- 
tinction models of nuclei, accumulation and coarse particles have been 
defined as linear combinations of the OPAC components proposed by 
Hess et al. (1998), which fit the mean values of a determined for differ- 
ent columnar aerosol types. Using these multimodal aerosol extinction 
models, instantaneous and diurnally averaged DARF effects induced by 
atmospheric aerosol at the TOA- and BOA-levels, and within the 


atmosphere have been calculated over sea and land surfaces by using 
the sets of BRDF surface reflectance models determined by Tomasi 
et al. (2013) for vegetation-covered, oceanic, and polar snow-covered 
surfaces. The results indicate that the diurnally averaged DARF terms 
vary strongly as a function of t( 0.50 pm), aerosol single scattering albe- 
do and surface albedo. The corresponding efficiency at the TOA-level 
was found to increase appreciably as surface albedo over sea and land 
increases until higher and more stable values over the snow- and ice- 
covered surfaces are attained. Less variable values and a more limited 
dependence of the BOA-level DARF efficiency on surface albedo were 
found, so that the atmospheric DARF efficiency resulted in: (1) a stron- 
ger increase when passing from vegetation-covered and oceanic sur- 
faces to surfaces covered only in part by snow, and (2) more stable 
positive values over the polar surfaces covered by fresh snow and 
clean glaciers in the Arctic and Antarctic regions. 
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